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Foreword 



It is a great personal pleasure to be introducing the 100th Volume of Structure 
and Bonding. The series was launched when I was an undergraduate and it 
proved to be an important resource for my initial exploration of inorganic 
chemistry and subsequently as a medium for the publication of some of my 
reviews when I developed a research interest in valence problems. Therefore, 
the series has served as a historical backdrop for my association with the 
subject. Initially the series was set up with the following intention; “A valuable 
service is performed by bringing together up-to-date authoritative reviews 
from the different fields of modern inorganic chemistry, chemical physics and 
biochemistry, where the general subject of chemical bonding involves a metal 
and a small number of associated atoms.... We are particularly interested in 
the role of the “complex metal-ligand” moiety and wish to direct attention 
towards borderline subjects.... We hope that this series may help bridge the 
gaps between some of these different fields and perhaps provide in the process 
some stimulation and scientific profit to the reader”. 

Of course in the last 30 years the subject has progressed and the initial 
reviews on ligand field theory, hard and soft acid base phenomena and 
bonding issues in cluster compounds have been progressively replaced by 
reviews which summarise and interpret interdisciplinary research at the 
borders of chemistry, biology and physics and reflect the breaking down of the 
barriers between inorganic, organic and physical chemistry. Therefore, in 1998 
the Editorial Board refined the aims of the series in the following terms: “We 
expect the scope of Structure and Bonding series to span the entire periodic 
table and address structure and bonding issues wherever they may be relevant. 
Therefore, it is anticipated that there will reviews dealing not only with the 
traditional areas of chemical bonding based on valence problems and 
dynamics, but also nano-structures, molecular electronics, supramolecular 
structure, surfaces and clusters. These represent new and emerging areas of 
chemistry”. 

With these aims in mind it is noteworthy that the current volume effectively 
reinforces and illustrates these ideals and is titled n-Electron Magnetism from 
Molecules to Magnetic Materials. Professor Veciana is to be congratulated on 
bringing together a series of reviews on organic molecular magnets. This area 
has considerable potential applications, but also presents considerable 
challenges for the synthetic, solid state and theoretical chemists. An 
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Foreword 



underlying theme is how systematic variations in the molecular structure of 
unsaturated organic molecules and the packing modes of the molecules may 
control their collective solid state magnetic properties. It is interesting how 
magneto-chemistry which historically was an area associated vdth transition 
metal ions in 1950’s has evolved into a topic of current interest to organic 
chemists. An understanding of the relationships between the packing modes of 
the molecules and the intermolecular magnetic interactions provides an 
equivalent challenge to the interpretation of the magnetic properties of 
transition metal complexes using crystal field theory half a century ago. 

May I take this opportunity of thanking Springer-Verlag and their editorial 
staff for all those behind the scenes activities which have enabled us to make 
this 100th volume a celebration and the members of the Editorial Board who 
have had the imagination to regularly propose new and exciting topics for 
individual volumes. It takes considerable enthusiasm to approach and badger 
colleagues to produce the manuscripts reasonably close to the submission 
dates. They may have lost friends in the process, but the scientific community 
generally has gained from the presence on library shelves of such a fine series. 

Oxford, February 2001 Michael Mingos 





Preface 



The notion of organic molecular materials with magnetic properties started 
three decades ago as a mere dream of some members of the chemical 
community. The goal was to create a collection of molecules, composed only of 
light elements (C, H, N, O, S), with unpaired electrons located in orbitals with a 
n symmetry, in which the whole material possesses the macroscopic magnetic 
properties exhibited by some classical metals or oxides - ferromagnetism, 
ferrimagnetism, etc. Such magnetic properties are unattainable from the 
isolated molecular building-blocks and took a couple of decades to be realized 
after its first proposal because of the subtle and complex structural and 
electronic interactions required by these materials. 

There is an awareness among scientists and technologists that organic 
magnetic materials offer considerable potential for applications in many 
diverse areas owing to the intrinsic characteristics of organic molecules 
(processability, transparency, flexibility, etc.). The skills of the organic chemist 
in allowing subtle variation in molecular structure of building blocks allow 
solid state characteristics to be fine-tuned. It is this rich potential for extreme 
diversity that gives organic molecular magnetic solids an important advantage. 
An underlying theme is, therefore, in the systematic variation of molecular 
structure and its effect on molecular packing and hence in the intermolecular 
electronic interactions that finally control the solid state magnetic property. 
An understanding of the relationship existing between packing arrangements 
in the solid and the intermolecular magnetic interactions are essential and 
many of the present developments in this area are linked to this aspect. This 
book contains a series of contributions chosen to highlight developments and 
methodologies in this direction of research. It is not intended to cover either 
all the topics of research in organic magnetic materials or all details of each of 
these topics. Thus, high-spin macromolecules have been deliberately excluded 
from this volume because they probably would require an entire volume of this 
book series. 

All topics contained in this volume deal with organic molecule-based 
magnetic materials in which the unpaired spins of the building-block units are 
located on 7i-orbitals conferring special characteristics that differentiate them 
from materials based on transition metal complexes. The first two chapters 
discuss the most representative examples of purely organic 3-D magnets as 
well as the through-space ferromagnetic interactions responsible of the bulk 
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An example of purely organic ferromagnets was first found nine years ago in the /!-phase 
crystal of the p-nitrophenyl nitronyl nitroxide radical (abbreviated as p-NPNN), which 
consists only of light elements, H, C, N, and O. This finding was prompted by extensive 
studies on the mechanism governing intermolecular ferromagnetic interaction in the 
galvinoxyl radical crystal. In the mid-1980s, only a few organic radicals were known to 
exhibit intermolecular ferromagnetic interaction. Galvinoxyl had been a typical example of 
such radicals. Strategies learned from galvinoxyl were first applied to p-NPNN, and the 
ferromagnetic transition was observed at 0.6 K. This article describes the progress of the 
research in this regards and the recent developments since then. The antiferromagnetism of 
the y-phase crystal of p-NPNN is also briefly described. 

Keywords: Organic ferromagnet, Intermolecular ferromagnetic interaction, Nitronyl 
nitroxide radicals, Spin polarization, p-NPNN 
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List of Abbreviations 



a.c. 


alternating current 


AFM 


antiferromagnetic 


CT 


charge-transfer 


FM 


ferromagnetic 


FOMO 


fully occupied molecular orbital 


galvinoxyl 


4-[[3,5-bis(l,l-dimetylethyl)-4-oxo-2,5-cyclohexadien-l- 

ylidene]methyl]-2,6-bis(l,l-dimethylethyl)phenoxy 


hydrogalvinoxyl 


4-[[3,5-bis(l,l-dimethylethyl)-4-hydroxyphenyl]methyl- 

ene]-2,6-bis(l,l-dimethylethyl)-2,5-cyclohexadien-l-one 


INDO 


incomplete neglect of differential overlap 


M 


magnetization 


MO 


molecular orbital 


NHOMO 


next highest occupied molecular orbital 


NLUMO 


next lowest unoccupied molecular orbital 


p-NPNN 


p-nitrophenyl nitronyl nitroxide (2-(4-nitrophenyl)-4,4,5, 
5-tetramethyl-4,5-dihydro- 1 H-imidazol- 1 -oxyl-3JM-oxide) 


SOMO 


singly occupied molecular orbital 


Tc 


critical temperature 


UHF 


unrestricted Hartree-Fock 


ZF-pSR 


zero-field muon spin rotation 


Z 


magnetic susceptibility 



1 

Introduction 

The study of magnetic properties of organic radical crystals has attracted 
much attention in the last few years. In particular, a movement to search out 
bulk ferromagnetism in organic materials composing exclusively of light 
elements was quite active in the 1980s. After the achievement of organic bulk 
ferromagnetism in 1991 [1], the magnetic properties of a variety of stable 
organic radicals such as phenoxy, nitroxide, and verdazyl radicals have been 
extensively studied. Among them, the derivatives of nitronyl nitroxide are very 
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stable and suited for the study of their solid state properties. The prominent 
feature of organic radicals is twofold. One is the highly isotropic nature of 
electron spins, which results from weak spin-orhit coupling in a molecule 
composed only of light elements such as H, C, N, and 0. The other is the 
distribution of spin densities over a molecule of Ti-electron system. These 
result in little anisotropy of the magnetic interactions. Therefore, the 
magnetism of organic crystals is, in most cases, analyzed on the basis of 
ideal Heisenberg spins. The distribution of spin densities also causes a 
complicated interaction scheme between adjacent radicals in a crystal, and 
gives rise to a variety of magnetic behavior. This is another feature of 
magnetism of organic crystals in comparison to that of inorganic crystals, in 
which localized character of electron spins is dominant. 

2 

Magnetism and Chemical Bonding 

Magnetism and the concept of chemical bonding are closely related to each 
other apart from the energy scale involved therein. In chemistry, a single bar 
indicates a single covalent bond, which corresponds to a pair of electrons 
coupled antiferromagnetically, that is with their magnetic moments in 
antiparallel fashion in accordance with the Pauli exclusion principle. In other 
words, a pair of electrons gets together in a small space between the atoms and 
becomes an adhesive agent for positively charged nuclei. Such a bond is quite 
stable in most cases, although these electrons repel each other in a small space. 
The stability comes from the overlap (or non-orthogonality) of the atomic 
orbitals associated with the bonded atoms. In the Heitler-London model for a 
hydrogen molecule (thus in terms of valence bond theory), it is the finite 
overlap integral which brings the stability to the spin singlet ground state. 

A similar situation occurs when two organic radicals are brought together 
closely, resulting in bond formation. An example is the triarylmethyl radical in 
solution, which yields hexaarylethane. In this case, a true covalent bond is 
formed between the methyl carbon atoms and the unpaired electrons are 
strongly coupled antiferromagnetically. From these and other examples, bond 
formation can be said to be equivalent to antiferromagnetic (AFM) spin 
coupling. 

There is the possibility of another kind of bond formation in the case of 
organic radicals. Some 7i-radicals, such as diphenylpicrylhydrazyl, nitroxides, 
phenoxyl derivatives, and many ion radical salts, are very stable even in the 
crystalline state. In these crystals, weak AFM couplings are usually recognized 
by magnetic measurements (i.e., the magnetic susceptibilities tend to become 
smaller at low temperature than those expected from the Curie law). From the 
analogy to the equivalence of chemical bond and AFM spin coupling, a weak 
chemical bond can be said to form between adjacent radicals in the crystal, 
where any atom in the radical species may not be specified to be involved in 
the bond formation. The bond is formed between radical molecules them- 
selves, i.e., the radical molecule as a whole behaves like an atom. These 
situations are illustrated in Fig. 1. 




4 



M. Kinoshita 




Fig. la, b. Schematic illustration of bond formation between radicals: a real bond formation; 
b weak intermolecular bond formation 



On the other hand, when the overlap between the orbitals of the unpaired 
electrons is small or zero, the triplet state or the state of parallel spin alignment 
is stable. This is the case with an oxygen molecule or carbene type radical. The 
physical basis of the spin parallel arrangement, in these molecules, is the Hund 
rule (or orthogonality). In the case of ferromagnetism, the spin magnetic 
moments on radical species align in a parallel fashion throughout the magnetic 
domain. Then, from the chemical point of view, no bond forms between 
adjacent radicals in the crystal. Therefore, to design an organic radical 
ferromagnet is, in a sense, a kind of chemistry that does not make a chemical 
bond between unpaired electrons on adjacent radicals. 

3 

Galvinoxyl 



3.1 

Magnetic Properties of the Neat Crystal 

Until the mid-1980s, only two organic radicals were known to exhibit 
intermolecular ferromagnetic (FM) interaction. One of them is the galvinoxyl 
(4-[[3,5-bis(l,l-dimetylethyl)-4-oxo-2,5-cyclohexadien-l-ylidene]methyl]-2,6- 
bis(l,l-dimethylethyl)phenoxy) radical, in which the ferromagnetic interac- 
tion was first observed by Mukai in 1969 [2]. In the inset of Fig. 2, the 
temperature dependence of the paramagnetic susceptibility of the galvinoxyl 
crystal is shown [3]. The susceptibility above 85 K follows the Curie-Weiss law 
with a ferromagnetic Weiss constant of 19 K. The crystal, however, undergoes 
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hydrogaMnoxyl 



Fig. 2a, b. Paramagnetic susceptibilities of the neat galvinoxyl crystal (inset) and of 6:1 
mixed crystal of galvinoxyl and hydrogalvinoxyl (main frame) 

a phase transition at 85 K and becomes almost diamagnetic below 85 K. This 
had caused difficulty in studying FM interaction in the temperature region of 
the Weiss constant. 

3.2 

Magnetic Properties of Mixed Crystals 

It was found that the phase transition is suppressed in making mixed crystals 
with a small amount of hydrogalvinoxyl (4-[[3,5-bis(l,l-dimethylethyl)-4-hy- 
droxyphenyl]methylene]2,6-bis(l,l-dimethylethyl)-2,5-cyclohexadien-l-one) 
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[3-5]. Hydrogalvinoxyl is a precursory closed shell compound for which the 
crystal structure is known to be isomorphous to that of galvinoxyl, but phase 
change does not occur [6]. The temperature dependence of the susceptibility 
of the 6:1 mixed crystal is shown in the main frame of Fig. 2. As is shown in 
Fig. 3, the reciprocal susceptibility of the mixed crystal crosses the temper- 
ature axis in the positive region, thereby confirming that the FM interaction is 
maintained even in the mixed crystal down to the cryogenic temperature. 

To comprehend the above matter in detail, the magnetic field dependence of 
the magnetization at 2 K was measured by making three types of the crystals of 
4:1, 6:1, and 9:1 mixing ratios [4]. The result is shown in Fig. 4. Depicted with 
the dotted curves are the theoretical curves of the Brillouin function for the 




0 20 40 60 80 100 120 140 

TIK 



Fig. 3. Reciprocal susceptibilities of the neat galvinoxyl crystal and 6:1 mixed crystal 




Fig. 4. Magnetization isotherms of the 4:1, 6:1 and 9:1 mixed crystals. The dotted curves 
represent the theoretical ones for S = 1, 2, 3, 4, and 5 
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spin quantum numbers S = 1, 2, 3, 4, and 5. When compared with the 
theoretical curves, it is seen that the experimental result for the crystals of n:l 
mixing ratio almost corresponds to the theoretical curve for S = nil. Therefore 
it has turned out that n pieces of the galvinoxyl radicals of S = 1/2 are amassed 
as an average, and their magnetic moments are aligned in parallel at low 
temperatures. 

As shown in Fig. 5, the crystal of galvinoxyl belongs to the monoclinic 
system (C2/c) and the almost planar radical molecules are arranged with the 
plane facing with each other along the c-axis at an equal interval [7]. If it is 
assumed that the one-dimensional interaction along the c-axis is effective, the 
above results are easy to understand. That is, the one-dimensional chain of 
radicals is divided into a number of segments comprised, on average, of n 
radicals separated by the closed shell molecules having similar molecular 
structure, and the magnetic moments of the radicals are aligned in parallel 
within the individual segments [5]. 

The magnitude of the exchange interaction was obtained by using 
paramagnetic resonance absorption [8]. In this case, we used a dilute mixed 
crystal containing about 4% of the radical and observed the temperature 
dependence of the absorption intensity of the Am = ±2 transitions, which are 
typical of the radical pairs existing statistically. The result is shown in Fig. 6. 
The ordinate, the product of the intensity and temperature, is the quantity 
proportional to the occupation number in a magnetic state (the triplet state in 
this case). The occupation number increases as the temperature decreases. 
Evidently the ground state is a magnetic state. The energy difference between 
the magnetic and non-magnetic states can be obtained as 1.5 meV by fitting 
to the theoretical equation for the susceptibility of radical pair. The value 




Fig. 5. Crystal structure of galvinoxyl. The tert-butyl groups are omitted 
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1.2 



S; 1 

a 



0.8 

0 20 40 60 80 100 

TIK 

Fig. 6. The temperature dependence of paramagnetic absorption intensity of the Am = ±2 
transitions observed in the dilute (4%) mixed crystal of galvinoxyl in hydrogalvinoxyl 

corresponds well to the Weiss constant obtained from the susceptibility 
measurements. 




3.3 

Mechanism of Ferromagnetic Interaction 



3.3.1 

SOMO-SOMO Overlap 

The mechanism governing the FM interaction is, of course, related to the 
electronic structure of the galvinoxyl radical and its crystal. We first examined 
the molecular orbitals (MO) of galvinoxyl by the INDO UHF method. 

As was mentioned above, the AFM interaction most frequently encountered 
among a number of organic radicals is explained by means of an analogy with 
the hydrogen molecules because two neighboring radicals have a tendency to 
form a weak covalent bond by making a pair of the unpaired electrons. For the 
formation of a covalent bond, overlap of the orbitals of the unpaired electrons 
plays an important role. With galvinoxyl where directions of the spins are 
prone to be in parallel, how is the overlap of the orbitals? This point was 
examined first. 

The intermolecular overlap integral was calculated with the singly occupied 
molecular orbitals (SOMO) using the molecular location of the crystal 
structure. The value of the SOMO-SOMO overlap became very small. When 
calculation was made by changing the relative location of the molecules step by 
step, it turned out that the overlap is close to the minimum when the 
molecular location is in the crystal structure. 

Furthermore, when the calculation process was examined it was found that 
relatively large positive and negative contributions cancel out each other, 
resulting in a small value. Since the 7i-orbitals are spread out over the 
molecule, the partial overlaps become positive and negative here and there and 
cancel out each other. As a result, the SOMOs on the adjacent radicals are 
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nearly orthogonal and degenerate. This is a situation quite similar to the case 
of the d-orbitals of an isolated transition metal ion. In the case of d-orbitals on 
one center, the partial overlap becomes positive and negative, but cancels out 
to null and the d-orbitals are orthogonal and degenerate. The exchange 
integral becomes positive in such a case and the high spin state is stabilized in 
accordance with the Hund rule. With galvinoxyl, even though it is a multi- 
center molecule and the intermolecular interaction is now considered, the 
Hund rule may be applied in a modified manner and the spin state must be 
stable when the spins are maintained in parallel. This is the so-called potential 
exchange. 



3.3.2 

Charge-Transfer Interaction 

The average distance between the neighboring galvinoxyl radicals along the 
c-axis is as much as 4.05 A, and the positive and negative overlaps themselves 
are not very large. Accordingly it appears that the interaction large enough 
to be 1.5 meV cannot be explained only by the static exchange interaction 
mentioned above. Some other mechanism, namely charge-transfer interaction, 
should also be taken into account. 

The 7 i-MO energy levels of galvinoxyl near the SOMO level are depicted in 
Fig. 7a [9-11], where the a-spin orbitals and the corresponding /1-spin orbitals 
are connected with dotted lines. The most conspicuous feature of this figure is 
that there exists the next highest occupied MO level of /1-spin (NHOMO-/1) 



a p 



A 

-f i". 



i" 



4 






-f- 

L-h i"J 



NLUMO 

SOMO 



-10 



- ^ 



NHOMO 



-15 



-20 






(a) 




Fig. 7. a Molecular orbital energy levels of galvinoxyl calculated by the INDO UHF method. 
Only the levels near SOMO are shown, h The electronic configurations in the radical pair 
coupled by CT interaction 
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situated higher than the SOMO of a-spin (SOMO-a). The exchange interaction 
within the molecule is great enough to stabilize the SOMO-a. In other words, 
the spin correlation causes a large spin polarization effect in galvinoxyl. 

In Fig. 7b, some low energy charge-transfer (CT) configurations are shown 
for the neighboring pair of radicals. NS and NT are the no-bond structures of 
singlet and triplet multiplicity and S, and T, (i = 0, 1, 2) are the excited CT 
state configurations, respectively. Among the excited CT configurations. So is 
the lowest and the resonance between So and NS usually stabilizes NS, resulting 
in an AFM interaction. This is the reason why most organic radicals exhibit 
AFM intermolecular coupling. With galvinoxyl, however, the SOMO-SOMO 
overlap is quite small as mentioned above. Therefore, the stabilization of NS 
by admixture of So is expected to be minimized. On the other hand, Ti and T 2 
must be lower in energy than Si and S 2 , respectively. In particular, Ti and T 2 
are much stabilized with respect to Si and S 2 in galvinoxyl because of the large 
spin polarization effect. Thus stabilization of NT by admixture of Ti and T 2 is 
expected to overweigh that of NS and the ground state becomes magnetic. It is 
in this way that the FM interaction is brought about in galvinoxyl. 

For the sake of assurance, the intermolecular overlap integrals, which are 
thought to be proportional to transfer integral, are calculated for the orbitals 
relevant to the above configurations with the molecular arrangement in the 
crystal. The results are given in Table 1 [9-11]. As is seen, the overlap integrals 
for the Ti and T 2 configurations are, respectively, larger than those for Si and 
S 2 . Thus it is concluded that the off-diagonal interaction is also favorable for 
the NT stabilization. 

From these observations, it is concluded that the cooperative effect of the 
spin polarization caused by a large exchange interaction within a molecule and 
the charge transfer interaction among molecules is essential for the FM 
interaction of the galvinoxyl crystal and probably of organic radical crystals in 
general [9]. The arguments described in this section may be summarized 
schematically as those shown in Fig. 8, where the states, S 2 and T 2 , are omitted 
for the sake of simplicity. In the case of FM coupling, the electrons in SOMO 
and NHOMO play an important role to make intermolecular bonds (Fig. 8a) 
because the overlap between them is large enough. In other words, the 
molecular design giving an opportunity of participating in a bond formation to 
the electrons in NHOMO, or more generally those in FOMOs (fully occupied 
MOs) is of crucial importance for FM interaction. The case of AFM interaction 
is shown in Fig. 8b, where the overlap between SOMOs are large enough to 
make a weak intermolecular covalent bond. 



Table 1. The calculated overlap integrals for the various CT configurations 



So 


(SOMO-a|SOMO-i?) 


0.72 X 10“^ 




(NHOMO-iS|SOMO-/?) 


1.60 X 10“^ 


Si 


(NHOMO-o£|SOMO-/1) 


0.87 X 10“^ 


T 2 


(SOMO-a|NLUMO-a) 


2.73 X 10"^ 


S 2 


(SOMO-o(|nLUMO-j3) 


1.33 X 10“^ 
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Dimer CT model for FM Interaction 
(large spin polarization) 



CT excited 
configuration 



No-bond 

configuration 



Non-magnetic (S=0) 



Magnetic (S=l) 




less stabilized 
ift-<SOMO|SOMO>-4) 



FM bonds in a chain (even in a solid) 



more stabilized 
because U' is small 




Dimer CT model for AFM Interaction 



Non-magnetic (S=0) Magnetic (S=l) 



CT excited 


r / -1 




SOMO 






configuration 






NHOMO 







resonance resonance 

t^/U (large) t^/U' (small) 



No-bond 


1 


— ^ 


SOMO 




^ — 


configuration 






NHOMO 







more stabilized (U small) less stabilized 



AFM bonds in a chain (solid) 



V 



(b) 







Fig. 8a, b. Schematic illustrations of; a ferromagnetic; b antiferromagnetic intermolecular 
bond formation. The bonds are shown by the elliptic circles. Note that the spins on NHOMO 
are aligned parallel in FM configuration, because the bonds are formed between the 
electrons in SOMO and NHOMO, while the intermolecular bonds are formed between the 
electrons only in SOMOs in the case of AFM configuration 
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4 

Conditions for Ferromagnetic Interaction in Organic Crystals 

In the preceding section we have examined the FM interaction in the 
galvinoxyl crystal. Now we extract the conditions for FM interaction applicable 
to other organic radicals from the study on galvinoxyl. 

As the property of radicals, it is first of all necessary for the spin 
polarization to be large enough. In other words, the exchange interaction 
within a molecule should be large enough. In order to obtain large exchange 
interactions in a molecule consisting of a number of atoms, it is 
recommended to use radicals having heteroatoms, such as nitrogen and 
oxygen. The unpaired electron resides largely on such electronegative 
heteroatoms, and interacts with the nonbonding electrons in them. This 
interaction is expected to be very large because of its one-center character 
and virtually determines the magnitude of the exchange interaction of the 
molecule, even though the unpaired electron is delocalized over the 
molecule. In galvinoxyl, the oxygen atoms at both terminals play this role. 
When the exchange interaction within a molecule is very large, the CT 
configuration of Fig. 7b is greatly stabilized and the energy of Ti is close 
to that of So- This means that the exchange interaction almost overwhelms 
the energy interval between SOMO and NHOMO; namely SOMO and 
NHOMO behave as if they were degenerate. In order for this to be realized, a 
molecule having an extended Ti-electron system is preferred because the 
interval between SOMO and NHOMO is related to the size of Ti-electron 
system. The extension of this argument would lead us to the situation where 
Ti is below So energetically; this is exactly the situation proposed by 
McConnell for FM coupling in 1967 [12]. 

Next is the problem of crystal structure. As mentioned above, the SOMO- 
SOMO overlap should be minimized. This presents a problem of where to 
locate the radicals in a crystal. Although it is not easy to control the crystal 
structure, it has been attempted by changing the position or size of chemical 
substituents and by introducing a charge or hydrogen bonding. However, the 
condition favoring ferromagnetism is herein pursued from the viewpoint of 
the electronic property of radicals. 

The galvinoxyl radical is known as a radical with large spin polarization in 
the field of chemistry, and has been studied as an example of a radical which 
possesses the so-called negative spin density. The negative spin density is 
usually observed in a radical of an odd-alternate Ti-electron system. The 
feature of an odd-alternate system lies in SOMO; SOMO has nodes 
alternately on the bonded atoms and lobes with opposite polarities on each 
of the neighboring atoms. Therefore, when two planar odd-alternate radicals 
are located close to each other in a face-to-face manner and are relatively 
deviated by one bond or an odd number of bonds, the local overlap is 
expected to become small through effacing with each other. In Fig. 9, a 
schematic diagram of SOMO of galvinoxyl and the relative location of the 
neighboring radicals are illustrated. It is easily seen that the local overlap is 
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Fig. 9. a Schematic drawing of SOMO and NHOMO of galvinoxyl (the fert-butyl groups are 
omitted). Filled and open circles show the difference in polarity of the p-orbitals. b The 
arrangement of adjacent molecules viewed along the normal to the best-fit molecular 
plane 



well effaced, resulting in a minimal SOMC-SOMO overlap throughout the 
whole of the molecules. Thus there is a high possibility of achieving such 
relative positioning of small SOMO-SOMO overlap by the use of radical of 
the odd-alternate system [9]. 

This situation may be evaluated from another point of view. Even though 
the charge density is zero, a finite spin density is induced at the nodes by the 
electronic correlation (or exchange interaction) effect. The induced spin 
density is negative in the sense that the direction of the induced spin is 
opposite to that of the unpaired electron. Since there are positive spin 
densities on the other atoms located at the lobes, the relative face-to-face 
location of two radicals deviating by one or an odd number of bonds 
corresponds to the situation where the atom of negative spin density in one 
molecule overlaps the atom of positive spin density in another molecule. 
This situation was also proposed for FM interaction by McConnell in 1963 
[13]. 

Another advantage of the use of the odd-alternate system is that there is a 
high probability of the overlap of SOMO with NHOMO or NLUMO being 
enhanced once such a relative location is established, because the charge 
density is usually distributed over the molecule in the tt-MOs other than 
SOMO. 

The arguments given in this section are summarized as follows. The radicals 
of the odd-alternate system having heteroatoms and relatively developed 
7i-conjugation provide a high possibility of realizing the intermolecular FM 
interaction in a crystal [9]. 
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5 

Ferromagnetism of p-Nitrophenyl Nitronyl Nitroxide 



5.1 

Electronic Structure of p-NPNN 

The first radical employed fulfilling the conditions in the preceding section is 
p-nitrophenyl nitronyl nitroxide (2-(4-nitrophenyl)-4,4>5,5-tetramethyl-4,5- 
dihydro-lH-imidazol-l-oxyl-3N-oxide, abbreviated as p-NPNN hereafter). Its 
large spin polarization effect has long been an object of research in the field 
of chemistry. When the molecular orbitals are actually calculated, the spin 
polarization effect manifests itself in a more obvious manner than in 
galvinoxyl; i.e., the fourth occupied MO for the ;S-spin maintains higher 
energy than SOMO-a, as shown in Fig. 10. 

Also, the charge density in SOMO is mostly concentrated on the two NO 
moieties and only a little is distributed over the other parts of the molecule. 
Accordingly, it is expected that the SOMO-SOMO overlap is minimized 
provided that the NO moieties of the neighboring radicals do not approach 
each other in the crystal (this actually holds in the /i-phase crystal, see below). 
On the other hand, as the charge distribution in other orbitals ranges over the 
whole molecule, the overlap with SOMO becomes large. These are quite 
favorable for FM interaction, as discussed above. 



5.0 



p-NPNN 



; SOMO 

i NLUMO 









4-th 

HOMO 



Fig. 10. Molecular orbital energy levels of p-NPNN calculated by the INDO LJHF method, 
showing large spin polarization effect. Only the levels near SOMO are shown 
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5.2 

Crystal Structure of p-NPNN 

There are four polymorphic forms, a-, P-, y-, and d-phases, known in p-NPNN 
[14-17]. Each of them can he separately prepared by properly adjusting the 
conditions for depositing crystals from solutions [18]. The crystallographic 
constants of these phases are given in Table 2 and the structure of the ^-phase 
is shown in Fig. 11. The orthorhombic /i-phase is the most stable form, and the 
other forms are subject to change to the /i-phase when they are kept at room 
temperature or below room temperature. Hereafter, details of the magnetism 
of the jS-phase are described. 



5.3 

Magnetic Properties of p-NPNN Crystal 

The magnetic properties of the /i-phase crystal of p-NPNN was first reported in 
1989 [19]. The temperature dependence of the paramagnetic susceptibility 
follows the Curie-Weiss law and gives a ferromagnetic Weiss constant of about 
1 K. In Fig. 12, the result of measurements with Hlla is shown [1, 18, 20]. There 
is no sizable anisotropy in the susceptibility, and nearly the same results are 
obtained, within the range of errors, in the other directions of the applied field. 

Since the Weiss constant 6 is very small, the FM interaction should also be 
checked by measuring the field dependence of the magnetization at low 
temperatures [1, 18, 20]. Figure 13 shows the results. The magnetization 
curves at several temperatures are unified info the single curve of the Brillouin 
function for S = 1/2 by the molecular field correction. In this case, the best fit 
is obtained with the coefficient of A = 2.8 Oe mole emu“^ (Heff - H + AM), 
yielding an FM Weiss constant off/ = 1.1 K in agreement with the result of 
susceptibility measurements. In addition, this experiment assures that the 
sample is not contaminated with an FM impurity. 



Table 2. Crystal structures of p-nitrophenyl nitronyl nitroxide 





a-phase 


/i-phase 




y-phase 


d-phase 


T/K 


=300 


=300 


6 


=300 


=300 


System 


monoclinic 


orthorhombic 




triclinic 


monoclinic 


Space group 


Plilc 


F2dd 




Pi 


P2i/c 


a/A 


7.302 


12.347 


12.16 


9.193 


8.963 


b/A 


7.617 


19.350 


19.01 


12.105 


23.804 


elk 


24.677 


10.960 


10.71 


6.471 


6.728 


a/[(7t/180) rad] 








97.35 




p/[(7i/180) rad] 


93.62 






104.44 


104.25 


>'/[(7i/ 180) rad] 








82.22 




Z 


4 


8 


8 


2 


4 


ViP 


1369.7 


2618.5 


2475.7 


687.6 


1391.3 


Density, p/g cm“^ 


1.354 


1.416 


1.498 


1.349 


1.333 
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Fig. 11a, b. Structure of the orthorhombic /J-phase crystal of p-NPNN: a molecules on the 
flc-plane; b view along the n-axis 
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r/K 



Fig. 12. The temperature dependence of paramagnetic susceptibility of the /!-phase crystal of 
p-NPNN 




(^+AA/;7’"VkOe K"^ 



Fig. 13. The field dependence of magnetization of the /?-phase crystal of p-NPNN at several 
temperatures. The curves are analyzed in terms of the molecular field approximation. 
The data almost coalesce into the curve of the Brillouin function for S = 1/2 with 
A = 2.8 Oe mol emu“^ 

5.4 

Transition to a Ferromagnetic Ordered State 

In Fig. 14, the temperature dependencies of the heat capacity and the a.c. 
susceptibility are shown. The heat capacity has a sharp peak at Tq = 0.6 K, and 
reveals the existence of a transition. The corresponding entropy amounts to 
85% of i?ln 2 in the range up to 2 K. Thus it is concluded that the transition is 
magnetic and bulk-like in nature. Since the a.c. susceptibility (shown in the 
inset) diverges around Tq, the ordered state is, without a doubt, a 
ferromagnetic state. 
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Fig. 14. The temperature dependence of heat capacity (main frame) and a.c. susceptibility 
(inset) of the /?-phase crystal of p-NPNN 



The magnetization curves measured at temperatures above and below the 
transition point are illustrated in Fig. 15 [1, 18, 20]. Although the curves have 
slight gradients at 1.22 K and 0.81 K in the paramagnetic region, the curve at 
0.44 K clearly traces a hysteretic loop characteristic of ferromagnetism. The 
magnetization is almost saturated at about 50 Oe, and the coercive force is 
small. The reason for small coercive force is the small anisotropy in the 
g-factor and the dipole interaction. The g-factors observed in the paramagnetic 
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Fig. 15. The field dependence of magnetization of the /!-phase crystal ofp-NPNN above and 
below the transition temperature (0.6 K). The field has been corrected for the demagne- 
tization 
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resonance are ga - 2.0070, gh — 2.0030, and g^ - 2.0106 [18]. The linewidth 
also shows an almost isotropic angular dependence at room temperature. 

5.5 

Other Evidence of Ferromagnetism 
5.5.1 

Field Effect on Heat Capacity 

Illustrated in Fig. 16 are the temperature dependencies of the heat capacity at 
various magnetic field strengths [18]. The sharp peak in the zero field is 
eventually slightly rounded and is shifted to the higher-temperature side as the 
magnetic field is increased. This is a feature of ferromagnetic materials. In the 
low magnetic field region, the demagnetizing field compensates the applied 
field up to a certain temperature, the substance behaves as if it were in the zero 
applied field, and the sharp peak remains as shown in Fig. 16a. 




Fig. 16a, b. The temperature dependence of the heat capacity of the /?-phase crystal of 
p-NPNN in: a the low; b the high magnetic field region 
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As the temperature increases above a certain temperature, however, the 
demagnetizing field becomes small and no longer cancels out the applied field. 
For FM substances, the critical temperature cannot be defined in a finite 
magnetic field. When there is an FM interaction among the spins, they have a 
tendency to align themselves in parallel along the magnetic field at very low 
temperature, and the spin system is ordered by a weak field even above the FM 
transition temperature. Thus the sharp peak of the heat capacity shifts and 
becomes rounded. This is shown in Fig. 16b. In the case of antiferromagnetic 
order, the peak remains up to a certain field strength determined by the AFM 
interaction (see Fig. 24). Therefore, this experiment ensures the ferromagne- 
tism of the /1-phase crystal below 0.6 K. 



5.5.2 

Zero-Field Muon Spin Rotation 

Another piece of evidence for ferromagnetism was obtained by the measure- 
ments of the zero-field muon spin rotation (ZF-juSR). Figure 17 shows some 
of the results of ZF-/zSR experiments conducted by Uemura and coworkers 
[21, 22]. The oscillating signals observed at 640 mK and 20 mK are due to the 
precession of the muons implanted into the crystal. Since there is no applied 
field, it is obvious the precession is caused by the internal field from the 
spontaneous magnetization. The long-lasting oscillation indicates that the 
muons experience a rather homogeneous local field which requires that the FM 
spin network is commensurate with the crystallographic structure. 

This is the result for the muon spin polarization perpendicular to the h-axis. 
When the muon spin is polarized along the h-axis, the amplitude of the 
oscillating signal becomes very small (about 20% of that of the perpendicular 
orientation). This suggests that the spin orientation in different domains is not 
aligned randomly and is most likely along the h-axis. Recent FM resonance 
experiments by Oshima et al. [23, 24] and neutron diffraction measurements 
by Schweizer’s group [25] also show that the magnetic easy axis is along the 
h-axis. 

In Fig. 18, the frequency of the oscillation, which is proportional to the 
spontaneous magnetization, is plotted against the temperature. The dotted 
curve shown is the one calculated using the random phase approximation 
under the assumption of a three-dimensional (3-D) isotropic Heisenberg 
model with a specific choice of 2 J/k^ = 470 mK for the interactions with the 
eight neighboring radicals. These results of the ZF-pSR experiments clearly 
demonstrate the appearance of spontaneous magnetic order in the /1-phase 
crystal of p-NPNN. 

Although the 2 J/k^ value of 470 mK is cited above only as an average, the 
exchange interaction is supposed not to be of only one kind. The crystal 
structure of the ;S-phase is schematically shown in Fig. 19, where each ball 
denotes the p-NPNN radical. The lattice may be divided into two face-centered 
orthorhombic sublattices deviating by a/4, b/4, and c/4 to each other, similarly 
to a diamond or, more precisely, to a zinc blende structure. From the crystal 
structure, we expect that at least two kinds of interactions dominate the 
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p-nitrophenyl nitronyl nitroxide 




Fig. 17. The ZF-fiSR spectra observed in the /i-phase crystals ofp-NPNN with initial muon 
spin polarized perpendicularly to the h-axis 

ferromagnetism. One is the interaction within the ac-plane and the other is the 
interaction between the neighboring ac-planes. The former interactions, 
shown by the dotted lines in Fig. 19, form an approximate 2-D square lattice 
[15], while the latter interactions, shown by the solid lines, result in a distorted 
tetrahedral coordination [20, 26]. According to the calculation by Okumura 
et al., the former interaction is estimated to be 2 / 12 /kB = 0.48 K, and the latter 
to be 2 = 0-22 K [27]. It has been suggested that the magnetic order in 

p-NPNN is due to a magnetic dipole interaction instead of the electronic 
exchange interaction [28]. However, our calculation of dipole-dipole interac- 
tion, D, based on the spin density data from the neutron diffraction 
experiments results in the values of DJk^ - -0.016 K, - -0.029 K, 

and DJk^ - 0.045 K. These values indicate that the dipole interaction in the fi- 
phase crystal stabilizes the spin system when the ferromagnetic spin alignment 
is assumed along the h-axis, but that it is smaller by one order of magnitude to 
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Fig. 18. The temperature dependence of muon spin precession frequency in the /!-phase 
crystals. The frequency is proportional to the spontaneous magnetization. The solid curve is 
a fit to M(T) = [l-(T/rc)“]^* with a = 1.86 and /! = 0.32 and the dotted curve represents 
calculation for isotropic 3-D Heisenberg model 




Fig. 19. Schematic drawing of the crystal structure of the /?-phase crystal of p-NPNN. Each 
ellipsoid represents the p-NPNN radical molecule. The three kinds of exchange interactions, 
/i 2 > /i 3 > and /i 4 , are identified 



account for the experimental Tq value. Therefore, we are of the opinion that 
the Curie temperature, To is essentially governed by the electronic exchange 
interaction, but that the direction of the magnetic easy axis is determined by 
the dipole interaction. 



5.6 

Pressure Induced Ferro- to Antiferromagnetic Transition 

When the pressure is applied to the crystal, the a.c. susceptibility exhibits 
remarkable changes. Figure 20 shows the pressure dependence of the a.c. 
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Fig. 20. The temperature dependence of a.c. susceptibility of the polycrystalline sample of 
the /1-phase p-NPNN as a function of applied pressure. The crystals undergo a ferro- to 
antiferromagnetic transition at around 650 MPa as shown in the inset 



susceptibility of the polycrystalline sample of //-phase p-NPNN up to 
p - 1.04 GPa measured under the a.c. field Hac(v) < 1 Oe (=10^/471 A m“^) 
with the frequency v = 15.9 Hz [29, 30]. It is found that the critical 
temperature, Tq, defined by the crossing point of the extrapolated lines from 
above and below Tq, approximately agrees with that determined from the heat 
capacity peak. As shown in Fig. 20, Tq shifts towards the lower temperature 
side with the initial gradient d[Tc{p)/Tc{po)]/^p = -0.48 GPa“^ and the 
magnitude of the susceptibility decreases gradually as the pressure increases 
from pQ (=0 MPa). In the low pressure region below p ~ 650 MPa, however, 
the ferromagnetic behavior is still preserved below Tdp), as characterized in 
the shape of Xac- 

In the high pressure region above 650 MPa, in contrast, the magnitude of 
Xac becomes quite small, the susceptibility in the ordered state decreases 
sensitively with the pressure increase, and the shoulder-like curve of Xac 
around Tdp) changes into a cusp as shown in the inset of Fig. 20. 
Furthermore, Tdp) increases as the pressure increases with the gradient of 
<^[Tdp)/Tdpc)]/‘^P - +0.04 GPa“\ where p^ - 650 ± 50 MPa is the critical 
pressure. These results suggest that the magnetic order below Tq is that of an 
antiferromagnet under high pressure. 
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The antiferromagnetic behavior is also recognized in the external field 
dependence of Xac at constant pressure, p = 690 MPa as shown in Fig. 21 
[29, 30]. Tcip) shifts to lower temperature as the field increases, in contrast to 
the case of a ferromagnet. Thus, we can conclude that pressurization induces 
a ferromagnetic to antiferromagnetic transition in the ;S-phase crystal. 



5.7 

Charge-Transfer Mechanism 

These experimental results could be explained in terms of charge transfer 
mechanism mentioned in Sect. 3.3 by competition between the ferromagnetic 
and antiferromagnetic interactions. The effective exchange interaction be- 
tween A and B molecules contains the kinetic (/fg) and the potential (/^g) 
terms as 

7ab = /ab + 7ab (1) 

Both terms depend on the overlap of molecular orbitals (MOs) on A and B 
molecules. The essential point of the charge transfer mechanism is that the 
kinetic exchange interaction is described by a sum of terms contributing to 
antiferro- and ferromagnetic interactions: 

p p . 

7^g = — ^ ^ 7™ + (terms related to other paths) (2) 




Temperature, T/K 



Fig. 21. The temperature dependence of a.c. susceptibility of the polycrystalline sample of 
the /?-phase p-NPNN in the antiferromagnetic region under 690 MPa as a function of applied 
field. The fields corresponding to the susceptibility peaks are plotted against the critical 
temperature in the inset 





An Organic Radical Crystal Showing Spontaneous Ferromagnetic Order 



25 



where fss stands for the transfer integral between SOMOs of A and B 
molecules, for that between SOMO and other FOMOs, U is the on-site 
Coulomb repulsion, and /“ is the intramolecular exchange integrals. Then 
interplay or frustration among these contributions would result in /ab ~ 0 in a 
certain condition, giving Tdp) = 0 K. In the case of the /1-phase of p-NPNN, 
we must take at least twelve interacting molecules adjacent to a central 
molecule in its zinc blende-like structure of Fig. 19. The corresponding 
exchange integrals are classified into three types / 12 , / 13 , and J 14 from the 
symmetry of the lattice. As mentioned before, J 12 and /13 are, from theoretical 
calculation, known to be ferromagnetic and /14 is only weakly antiferromag- 
netic. On the other hand, reduction of the dimensionality from the three- to 
two-dimensional ferromagnetic system has been pointed out from the 
appearance of short-range order effect by the heat capacity measurements 
under pressure [31]. From the crystal symmetry, it is obvious that only the 
exchange integral J 12 is responsible for the two-dimensional ferromagnetic 
interaction. The transition temperature, T 3 D, in such a reduced system, can be 
written in terms of the mean field theory as 

^bTsd OC S^<f2D(fl2, T3 d){|/i3 +/ 14 I}, (3) 

where ^ 2 d is the spin correlation length in the nc-plane in which / 12 /k ~ 0.8 K 
is estimated from the heat capacity curve [30, 31]. Therefore, the antiferro- 
magnetic behavior zX p > p^- 650 ± 50 MPa can be ascribed to a change in 
the sign of / 13 . This means that the relative importance of the first and second 
terms in Eq. (2) relieve each other for /13 under high pressure. 



5.8 

Lattice Constants Under High Pressure 

The change of the magnetic interactions discussed above would be closely 
related to changes in molecular shape and packing in the crystal. The lattice 
constants under various pressures have been determined by the Riedvelt 
method from the X-ray powder patterns obtained on the polycrystalline 
sample of ;S-phase p-NPNN by the use of imaging plates. The powder patterns 
observed are similar to one another, indicating that the crystal symmetry is 
maintained under pressure up to 1.26 GPa [30, 31]. The peak positions, of 
course, shift toward the direction of wider angles or to the direction of lattice 
contractions as the pressure increases. The lattice constants are plotted against 
pressure in Fig. 22. The crystal changes in two steps. The linear and volume 
compressibilities are summarized in Table 3. The biggest contraction ('^4.5%) 
is found again along the c-axis as in the case of thermal contraction. The 
crystal density increases up to as high as 1.58 g cm“^ at 1.26 GPa. 

The two-step contraction observed here could be understood in the 
following way. The molecular shape at 6 K is compared with that at room 
temperature in Fig. 23 [32]. The molecules are rotated by ±3.3 (7i/180) rad 
from the orientation at room temperature, the nitro groups being further 
rotated by ±1 (7r/180) rad. The lattice constants at 6 K (see Table 2) [25] nearly 
corresponds to those at about 400 MPa. From this, we could expect that the 
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Fig. 22. The pressure dependence of the lattice constants of the /i-phase crystal of p-NPNN 



Table 3. Compressibility of the /?-phase p-nitrophenyl nitronyl nitroxide crystal in units 
of Pa“^ 





K — — 


Ku — — 

~ bdP 


K — — 

“ cdP 


hr - -1^ 
~ V dP 


P < 550 Pa 


34.4 X 10^“ 


56.6 X 10“" 


55.1 X 10“" 


143 X 10“" 


P > 550 Pa 


11.5 X 10““ 


18.4 X 10“" 


22.3 X 10“" 


49.0 X 10“" 



molecules are librationally rotated about the a-axis to some further extent 
under pressures up to about 550 MPa. However, the plane of the nitro group 
has already been at the upright orientation at 6 K (or equivalently at about 
400 MPa) with respect to the ac-plane, and further pressurization beyond 
550 MPa would not increase the tilt angle any more for the nitro group. Then, 
the other parts of the molecule would rotate internally about the long 
molecular axis under higher pressures. Internal rotations of the five- 
membered ring having the bulky tetramethylethylene group would be the 
most probable candidates to take place in order for the molecule to become 
more and more planar and for the crystal to become more compact. 

Such molecular deformation would cause a change in the electronic 
structure of the molecule and in the intermolecular magnetic interactions. At 
this moment we cannot conclude definitely, because the crystal structure at 
very low temperature and under high pressure is not available. However, it is 
natural to expect that the molecules take a more planar form and a more 
parallel arrangement in the crystal on compression at low temperature. The 
intermolecular charge transfer interaction between the SOMOs on adjacent 
molecules would become more efficient, resulting in interchange of the 
importance of ferromagnetic and antiferromagnetic interactions, as discussed 
in the preceding section. 
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Fig. 23. The molecular shape of p-NPNN in the /!-phase crystal at room temperature and at 
6 K 



6 

Antiferromagnetism of y-Phase Crystal of p-NPNN 

It was initially suggested that the y-phase crystal of p-NPNN can become a 
ferromagnet because the heat capacity exhibited a sharp peak and the a.c. 
susceptibility became quite large at 0.65 K [33]. However, it was later found 
that the heat capacity peak remains even in the magnetic field up to 1800 Oe, 
as is seen in Fig. 24 [34]. This is in contrast to the case of the /1-phase and is 
characteristic of an antiferromagnet. Despite that, it was not established why 
the a.c. susceptibility became so large. Various efforts have made it clear that 
the phenomenon is induced by contamination of the ferromagnetic ;S-phase, to 
which the y-phase crystal transforms during the measurements at low 
temperatures [34]. 

After various trials to determine the proper experimental conditions, it 
became possible to measure the a.c. susceptibility of the y-phase itself with a 
powder sample. Figure 25 shows the results of a.c. susceptibility measure- 
ments with the powder sample. The curve at the zero field shows the 
characteristic feature of antiferromagnets, namely peaking at the transition 
temperature and decreasing to about 2/3 of the maximum in the low 
temperature side. Thus it is concluded that the y-phase crystal of p-NPNN 
undergoes a transition to an AFM order at = 0.65 K. 
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Fig. 24. The temperature dependence of heat capacity of the y-phase crystal of p-NPNN as a 
function of applied field. Note the heat capacity peak persists up to the field as high as 
1800 Oe 
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It should be noted that the magnetic susceptibility and the heat capacity of 
the y-phase are well explained by a 1-D isotropic FM Heisenberg model with 
2 J/k^ = 4.3 K in the temperature range above 4 K [35]. This FM interaction is 
believed to operate along the [Oil] direction. For the other two directions, FM 
coupling with 2 J'/k^ = 0.22 K and AFM coupling with 2 ]"!k^ - —0.18 K are 
suggested, from the analysis of the field dependence of the heat capacity. 
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Fig. 26. Some organic ferromagnets reported 
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7 

Other Ferromagnetic Organic Crystals 

Following the finding of the ferromagnetism in p-NPNN, a dozen organic 
compounds have been found to become a ferromagnet [36-44]. These are 
shown in Fig. 26. Most substances undergo the ferromagnetic transition below 
1 K. Diazaadamantane dioxide, which has the triplet ground state, exhibits the 
highest transition temperature of 1.48 K [38]. The transition temperature Tc is 
proportional to S(S + 1) in the mean field approximation, and this means that 
the magnitude of the ferromagnetic interaction is nearly the same as that of the 
other compounds. 

In order to increase the ferromagnetic transition temperature, it is 
indispensable to introduce other mechanisms other than those described 
here. Introduction of conductivity is one of the possibilities and various trials 
are being attempted in several laboratories. 
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The mechanism of the intermolecular magnetic (or through-space) interactions found in 
purely organic molecular crystals is reviewed, using results from the ab initio studies of 
model systems, and from the statistical analysis of the packing of nitronyl nitroxide crystals 
presenting dominant ferro or antiferromagnetic interactions. First of all, the foundations of 
the McConnell-I mechanism (more properly called a proposal) are reviewed from a rigorous 
theoretical point of view. It is shown that this proposal lacks a rigorous foundation and 
works in some prototype systems due to error compensations associated to the high 
symmetry of the model systems employed to evaluate such a mechanism. It will be shown 
how the McConnell-I mechanism fails in rationalizing the magnetic properties of well 
known nitronyl nitroxide crystals. We will also show the existence of pitfalls in many of the 
magneto-structural correlations employed today. One example of an erroneous correlation 
is that which associates the magnetic character of a nitronyl nitroxide dimer with the relative 
orientation of the ONCNO groups in these two molecules. Consequently, new magneto- 
structural relationships are needed based on unbiased assumptions. For such a purpose, we 
need to have solid data on the properties of the through-space magnetic interactions at the 
microscopic level. We will review in the last sections the current knowledge about these 
properties, obtained from ab initio computations and crystal packing studies. 

Keywords: Molecular magnetism, Nitronyl nitroxides, Though-space magnetic interactions, 
Magnetic interaction mechanism, McConnell-I mechanism, Quantum chemical ab initio 
studies. Statistical analysis 
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1 

Introduction 

Molecular magnetism is a macroscopic property found in solids having long- 
living, spin-containing units (that is, stable molecular radicals, in the case of 
purely organic molecular magnets) oriented in such a way that they allow the 
propagation of the magnetic interactions throughout the solid [1]. In the case 
of purely organic molecular magnets, the magnetic interaction occurs mainly 
by direct overlap of the orbitals of the nearby units, in what are called through- 
space interactions, to distinguish them from the through-bond interactions 
found in organometallic solids, where the spin-containing metallic centers 
interact by means of the orbitals of the ligands placed between these centers 
(this is one example of the superexchange mechanism) [la, b, 2]. 

The properties of the mechanism of the through-space magnetic interac- 
tions are not well understood at the present moment, although some basics 
facts are well established. Thus, for instance, the dependence of the magnetic 
interaction on the relative orientation of the radicals within the crystal is well 
known. Such a property is clearly demonstrated by the different magnetic 
properties found in different polymorphic forms of the same radical [3] 
(which, by definition, can only differ in the relative orientation of the same 
radicals). However, for a proper rationalization of the current molecular 
magnets or for a rational design of new, purely organic materials showing 
magnetic properties, one needs to go one step further and be able to identify 
which relative orientations of the radicals give rise to ferromagnetic 
interactions in the crystal, and which ones generate antiferromagnetic 
interactions. One should mention here that having good magneto-structural 
correlations is only half of the problem in designing molecular magnets. In 
addition, one has to learn how to control the packing of the radicals within 
the crystal [4, 5] to guarantee the presence of the desired magnetic 
interactions in the crystal. This is an example of what is nowadays called 
crystal engineering [6]. We will not discuss the complexities of crystal 
engineering in detail here. Instead, our aim is to review the current 
knowledge on the mechanism of the through-space magnetic interactions, at 
a microscopic level. 

At the present time, there are some empirical magneto-structural 
correlations published in the literature, generated after looking for common 
factors in the relative orientations of the radicals whose crystals show 
dominant magnetic properties. Thus, for instance, it is commonly accepted 
that short N — O- • 0 — N contacts in the nitronyl nitroxide crystals are 
indicative of dominant antiferromagnetic interactions between the dimers in 
which such contacts are found [7]. Besides these empirical correlations, the 
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most popular magneto-structural correlation is the so called McConnell-I 
proposal (or mechanism) [8], which is based on theoretical considerations. It 
receives this name because it was the first of two proposals presented by 
McConnell to describe the mechanism of the through-space magnetism [9]. 
Originally introduced to rationalize the through-space magnetism in crystals 
whose packing motif are piles presenting n-n short contacts, it is nowadays 
used on any kind of molecular crystal. It is a very popular mechanism 
because it is easy to use: the rationalization of the crystal magnetic 
properties is done by looking at the products of the atomic spin populations 
of the atoms making the shortest intermolecular contacts between pairs of 
adjacent radicals [1]. 

For a long time, the McConnell-I mechanism was believed to work [1], in 
particular after its predictions had been successfully matched against the 
experimental results obtained on the [2.2]paracyclophane isomers [10]. 
However, there are also cases in which the McConnell-I predictions are in 
conflict with the experimentally measured magnetic properties. Here, we will 
illustrate some of these conflicting cases. Detailed theoretical and numerical 
investigations of the McConnell-I proposal showed that these failures are just a 
consequence of the deficiencies of the proposal: the basic expression in which 
the McConnell-I proposal is based is not theoretically sound [11]. Moreover, it 
was shown that the success of this proposal in the [2.2]paracyclophane 
isomers came from a fortuitous cancellation of errors, due to the high 
symmetry of these isomers [11]. Therefore, such a cancellation is not likely to 
take place in the general case. 

Given this state of affairs, one needs to find new magneto-structural 
correlations based on sound principles. This can be done by investigating the 
microscopic properties of the through-space magnetic interactions along two 
lines of research: 

(a) searching for general trends on the relative orientations of neighboring 
radicals in the known molecular magnets, by detailed statistical analysis of 
their crystal packings, and 

(b) by quantum chemical computations on model systems, aimed at repro- 
ducing the type of magnetic interactions present in the purely organic 
molecular crystals. 

Here we will present the main findings from investigations in which the two 
lines of research are used. 



2 

The McConnell-l Mechanism of Through-Space Magnetism: 

A Theoretical Analysis 

In order to explain the intermolecular magnetism between aromatic radicals, 
in 1963 McConnell suggested [8] that the magnetic interaction between two 
aromatic radicals A and B could be approximated by the following Heisenberg 
spin Hamiltonian (Eq. 1) 
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Jfsfsf ( 1 ) 

ieAJeB 

in which are two-center exchange integrals, whose form is (Eq. 2) 

Jij = [Wj\ + Mi\j){i\h\j) (2) 

and SfS^ is the product of the spin operators of atoms i and j placed on 
fragments A and B, respectively. One should note that this expression is an 
approximation of the general Heisenberg Hamiltonian (Eq. 3) [12] 

= + (3) 

in which the intrafragment terms have been neglected. Such an approximation 
is only valid if the expected values of the intrafragment components are the 
same for the states for which the general Hamiltonian acts. 

McConnell [8] also proposed that the Hamiltonian of Eq. 1 could be 
replaced by a simplified form (Eq. 4) 

Jfpfpf (4) 

In this expression the first two operators are the total spin operators for 
fragments A and B, are same two-center exchange integrals defined above, 
and pfpj is the product of the atomic spin density on atoms z and j of 
fragments A and B, respectively. This Hamiltonian is purely phenomenolo- 
gical, as there is no systematic set of approximations that derives it from Eq. 1 
in a rigorous way [11]. 

One can gain a better insight into the properties of the McConnell-I 
Hamiltonian (Eq. 4) by computing energy differences between two states. If 
only one unpaired electron is present in each fragment, these electrons can be 
paired into a singlet (S) and a triplet (T) state. Using the following expression 
(Eq. 5) 

= i[S(S+ 1) - Sa{Sa + 1) - Sb{Sb + 1)] (5) 

it is possible to demonstrate that the energy difference between the singlet and 
triplet states obtained using the Hamiltonian of Eq. 4 is equal to Eq. 6 [11] 

E^-E'^ = jfpfPj = J (6) 

where J is an effective coupling constant between the two states of interest. 

According to this equation, it can be concluded that the McConnell-I 
mechanism predicts that a triplet ground state is obtained in the interaction 
between two doublet fragments when the atoms making the shortest contacts 
present atomic spin populations of opposite sign (as is normally negative). 
Such an interpretation rests on the assumption that the values of the /d® 
integrals are determined mainly by the intermolecular atom- atom distance and 
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are all similar. However, given the complexities of the molecular orbitals in 
which the interacting electrons are located, this is not necessarily the case and 
these integrals could be affected by other factors like the symmetry of the 
orbitals, or the relative position of the atoms in their fragments, for instance. 
Therefore, it is possible that the shortest contacts are not those presenting the 
dominant values. As a consequence, the association is not always valid 
[ 11 ]. 

Another weak point of the “triplet-opposite signs” interpretation comes 
from the validity of Eqs. 4 and 6. As mentioned above, Eq. 4 is a simplification 
of the general case, whose validity has been never demonstrated. In fact, the 
derived expression for the energy difference obtained using the 

general Heisenberg Hamiltonian (Eq. 3) is different to that obtained with 
the simplified form (Eq. 7) [11] 

E^-E^ = J2lij{Plj-Pt) ( 7 ) 

‘j 

where P| and Pjj are the singlet and triplet exchange density matrices obtained 
from the singlet and triplet wave functions, using the following Eq. 8 

Pij=(^-(^2SrSj+^-Iii^'j ( 8 ) 



Therefore, it is clear that the McConnell-I relationship is valid only if one 
can make the following association (Eq. 9): 



Pfpf ^Pii-Pn 



( 9 ) 



However, there is no reason why the product of atomic spin densities should 
be related to the difference of exchange density matrices. In the best case, this 
association has to be tested numerically before accepting it. 

In order to test the applicability of Eq. 9, we carried out numerical 
computations on simplified models of the ortho-, meta-, and para-isomers of 
the bis(phenylmethylenyl)[2.2]paracyclophanes [10] (see Eig. 1), which rep- 
resent their main structural and electronic properties. These computations 
were carried out by solving the exact Heisenberg Hamiltonian in a Valence 
Bond configurational space [11], using the so called MMVB method [13]. This 
method uses the MM2 potential to describe the inert cr-bonded framework, 
and a Heisenberg Hamiltonian to represent all the electrons involved in the n- 
conjugation or in the formation of new c-bonds. The Heisenberg Hamiltonian 
has been parametrized to reproduce the results from accurate CASSCF 
computations on various model systems [13]. 

The ground states for the three model paracyclophane isomers agree well 
with the experimental results and the McConnell-I predictions. An analysis 
of the energetic components in these computations indicates that in the three 
isomers the dominant contribution to the singlet-triplet energy difference is 
that coming from the atoms making the shortest contacts, as one would 
expect if the McConnell-I proposal would work. The reason for this 
domination is the small numerical value of the integrals associated to the 
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Fig. 1. Structure of the model employed to evaluate the properties of the [2.2]paracyclo- 
phanes in our MMVB computations. The structure corresponds to the ortho isomer 



Other terms present in Eq. 7, due to the high symmetry of the system. 
Therefore, on systems where such a high local symmetry environment is 
destroyed and the atoms are not well aligned, one should not expect that this 
cancellation will hold. Test computations on distorted paracyclophanes 
support such a conclusion. Consequently, there is no reason to expect that 
the McConnell-I proposal will hold in general. Our computations also showed 
that the /|^ integrals have strong orientational and directional properties. 
Thus, even if the proposal would hold, a precise determination of the singlet- 
triplet stability cannot be done without a computation of the /d® integrals, 
in contrast to common practice which makes the predictions solely on the 
basis of the atomic spin populations. 

3 

Failure of the McConnell-l Mechanism in Nitronyl Nitroxide Crystals 



3.1 

Evidence from Individual Crystal Studies 

Even if not all the interactions present in the crystals of purely molecular 
radicals are of the n-n type, the McConnell-I mechanism has been the main 
tool used in the analysis of the magnetic interactions present in molecular 
crystals. However, there is now enough evidence indicating that the McCon- 
nell-I predictions are not always valid. This is the case for many magnetic 
crystals of the nitronyl nitroxide family, whose general formula is indicated in 
Fig. 2 [ 1 ] . Here we will describe some of the evidence indicating the failure of 
the McConnell-I mechanism. 

A first evidence against the validity can be obtained when analyzing the 
magnetic properties of many nitronyl nitroxides with well defined structures. 
We can illustrate the failure with one of these examples, the phenylnitronyl 
nitroxide crystal (PhNN, see Fig. 3) [14]. The X-rays structure of this crystal 




The Mechanism of the Through-Space Magnetic Interactions 



39 



CH, 



/. 

.N 

\ 

/ 



Fig. 2. Chemical structure of a general nitronyl nitroxide radical (R is the substituent which 
changes from one radical to another in this family) 




Fig. 3. a Structure of the A plane of the PhNN crystal, b Structure of the B planes of the 
PhNN radical. Both planes extend along the crystallographic b-c plane, and are visualized 
along the a axis 
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has been determined with precision [14], and so was the spin density 
distribution using polarized neutron diffraction [15]. The crystal belongs to 
the P2i/c symmetry group, with cell parameters a = 21.14 k,b- 10.14 A, c = 
12.122 A, and P = 108.1°, and eight molecules per unit cell. The radical packs 
building planes on the b-c directions, which pile up in the a direction. There 
are two types of planes, the A and B planes shown in Fig. 3. Four planes pile up 
within each unit cell in an ABBA sequence. The experimental spin distribution 
agrees well (see Fig. 4) with that computed using the BLYP non-local exchange 
and correlation density functional [16, 17] and the 6-31G(2d,2p) basis set [18]. 
These results show that 90% of the spin density is located on the ONCNO 
atoms of the five-membered ring of the PhNN molecule (Fig. 4), being positive 
on the N and 0 atoms (0.29 and 0.35 e~, respectively), and negative on the 
central C atom (-0.22 e~). The remaining 10% of the spin density is 
distributed on the other atoms, with atomic spin populations in each atom 
smaller than 0.03 e“ in absolute values. There is sign alternation between 
adjacent C atoms of the aromatic ring, and also between the N and C atoms of 
the five-membered ring. The amount of atomic spin population on the 
hydrogen atoms is always smaller than 0.003 e~ in absolute values. It is 
important to mention here that the atomic spin population is not always 
proportional to the spin density at the nucleus, as obtained from the NMR or 
EPR hyperfine coupling constants (hfccs) [19]. This fact is particularly 
important on some of the hydrogens of the PhNN radical, and is inherent to 
the fact that the atomic spin population is an integration of the spin density on 
the nucleus and the rest of the space associated to that atom [20]. 

The PhNN crystal presents dominant antiferromagnetic properties [14]. 
However, according to the McConnell-I mechanism one should expect this 
crystal to be a paramagnet. This conclusion is reached after analyzing all the 
direct X- • Y short contacts between the atoms of the ONCNO group. Note 
that when extending the McConnell-I mechanism to a interactions the short 
H- • O or H- • C contacts are disregarded because the atomic spin populations 




Fig. 4. Representation of the atomic spin population found in the atoms of the PhNN radical. 
Shaded atoms represent positive values, and unshaded negative ones. The size of the circles 
reflect the amount of atomic spin population (larger circles mean more population). The 
population on the H atoms is too small to he drawn 
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located in the H atoms are two orders of magnitude lower than the values in 
the ONCNO atoms. The shortest contacts between ONCNO groups within the 
A and B planes are nearly collinear O- • 0 contacts located at 5.02 A and 
6.65 A. Thus, given the large distance of these contacts, paramagnetism 
within the A and B planes should be expected if McConnell-I works properly. 
The situation for the interplane contacts is similar: the shortest 

ONCNO- • ONCNO contacts found in between adjacent A-A planes is 
4.51 A, which becomes 4.05 and 3.69 A when looking between adjacent A-B 
and B-B planes, respectively. The A-A and A-B distances are too large to allow 
an effective interplane magnetic interaction to be propagated all over the 
crystal. Therefore, as already mentioned, the McConnell-I predicts that the 
PhNN crystal should behave as a paramagnet, if one looks at the direct 
overlap of the ONCNO atomic spins. 

The failure of the McConnell-I mechanism to predict magnetic properties of 
the PhNN crystal in good agreement with the experimental data, is only one of 
the many cases found in the nitronyl nitroxide crystals. Among them, there are 
various hydrogen bonded crystals whose short contacts are always of the 
A — H- • B hydrogen bonded type. In such crystals, usually one finds that the 
ONCNO- - -ONCNO distances are very large and the McConnell-I mechanism 
systematically predicts a paramagnetic behavior. Therefore, to be able to 
explain the existence of the ferromagnetic or antiferromagnetic behaviors 
detected experimentally, the existence of magnetic interactions throughout 
the hydrogen bond is postulated [21]. This is a clear transgression of the 
McConnell-I mechanism. However, the existence of through-hydrogen bond 
magnetic interactions on hydrogen-bonded radical dimers is supported by 
theoretical computations on H 2 NO dimers [22]. 

When one analyzes the magnetism of the PhNN crystal after allowing for the 
existence of magnetism through the C — H- - -ON contacts, the crystal is no 
longer a paramagnet: each PhNN molecule makes many C — H- - - ON contacts 
with other adjacent radicals in all directions of the space, due to the fact that 
the C — H- - -ON contacts are the strongest attractive contacts in the PhNN 
crystal in conjunction with the C — H- - -n contacts [5b]. Therefore, if one 
accepts that the McConnell- 1 mechanism can be extended to the C — H- - -ON 
contacts, and leaves aside the small value of the atomic spin population found 
in all hydrogens (this is equivalent to saying that the Jy integrals for these 
interactions are stronger), it is possible to know if the contact is ferromagnetic 
or antiferromagnetic looking at the sign of the atomic spin population located 
on the hydrogen of the C — H group. A simple qualitative approach to define 
the sign of the atomic spin population on such hydrogens is to assign to 
that atom the opposite sign to that found in the carbon of the C — H group 
(see Fig. 4 for the atomic spin populations). Using this approach, all 
C(sp^) — H- - -ON contacts should be ferromagnetic, while the magnetic charac- 
ter of the C(sp^) — H- - -ON contacts depends on the aromatic C atom: H atoms 
attached to the Cmeta atoms give rise to ferromagnetic C(sp^) — H- - -ON con- 
tacts, while these contacts are antiferromagnetic in the Cortho and Cpara case. 
The overall nature of the interaction is difficult to define, as there are many 
contacts of opposite character within hydrogen-bonded crystals. Further 
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studies are under way to determine the importance of these interactions in real 
crystals. 



3.2 

Evidence from ab initio Computations on Small Model Radical Dimers 

The results of the previous section show the existence of inconsistencies 
between the experimental magnetic character and the predictions of the 
McConnell-I mechanism. However, these studies do not give information 
on the cause of the failure. One form of obtaining such information is 
by comparing the McConnell-I predictions on simple radical dimers 
with the results obtained from accurate ab initio computations on the same 
dimers. 

We started our investigation on model dimers, by comparing the ab initio 
and McConnell- 1 predictions for the ground state of an H 2 NO dimer oriented 
in such a way that the shortest contacts between the fragments are established 
between the oxygen atom of each fragment (Fig. 5). We made a potential 
energy surface scan in which the relative position of one fragment respect 
to the other was modified while keeping the 0- • O distance fixed at 3.0 A, also 
forcing the N — 0 axis of one of the H 2 NO fragment to be collinear with the 
oxygen atom of the other fragment. This scan was done while keeping all 
atoms in the plane of the fixed fragment (the a-b plane of Fig. 5). Alternatively, 
the NO group of the second molecule can be moved along the a-c plane of 
Fig. 5, while the hydrogens are in symmetrical position above and below that 
plane. In the two sets of geometries, the McConnell-I mechanism predicts that 
the interaction between these radicals must be antiferromagnetic, as the 
density in the two O atoms has the same sign (in this case, positive, like the N 
atom). 

The singlet-triplet energy difference of the H 2 NO dimers computed at the 
CASSCF(6,4)/6-31G(d) level for angles in the 90-180° range are shown in Fig. 5 
[23]. By CASSCF(6,4)/6-31G(d) level one indicates computations done using 
the CASSCF method with a large (6,4) complete active space (that is, built by 
putting four electrons into six orbitals in all possible ways) and the 6-31G(d) 
basis set. We tested that the singlet-triplet energy separations did not change 
when the size of the active space or the basis set was increased. Therefore, 
there are regions of the dimer in which the singlet is the ground state and others 
in which it has a triplet ground state. This dual behavior depending on 
the relative orientation, goes against the predictions obtained using the 
McConnell-I model in the usual form, that is, by looking only at the atomic 
spin populations of the overlapping atoms. 

The disagreement between the McConnell-I predictions and the ab initio 
results found in the H 2 NO dimer can be extended to other systems, and, 
in particular, to the methyl- • allyl radical dimer. Here, while the McConnell-I 
predictions are correct when the plane of the allyl radical is parallel to the 
plane of the methyl radical [23], it fails to predict the rapid change from 
ferromagnetic to antiferromagnetic when the methyl plane is tilted off the 
parallel arrangement. 
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(b) <0...0-N / degrees 

Fig. 5. a Diagram representing the orientation of the H 2 NO- • ONH 2 dimers in the scans 
along the a-b or a-c planes (see text for details), b Energetic difference between the singlet 
and triplet states as a function of the O- • O — N angle 

3.3 

Evidence from Statistical Magneto-Structural Studies in Nitronyl Nitroxide Crystals 

The previous two examples (Sects. 3.1 and 3.2) could be taken as exceptional 
cases, which, therefore, do not reflect the overall trend found in the 
experimental crystals. Therefore, it seems appropriate to investigate if it is 
possible to correlate, in a general statistical form, the magnetic properties of 
the crystal with the relative geometrical orientation of the adjacent ONCNO 
groups present in the nitronyl nitroxide crystals. According to the 
McConnell-I mechanism, crystals presenting dominant ferromagnetic inter- 
actions should present a different geometrical arrangement of the ONCNO 
groups than those found in crystals presenting dominant antiferromagnetic 
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interactions. We thus selected those crystals presenting dominant magnetic 
properties, either ferro- or antiferromagnetic, because in this way we can 
assume that all their ONCNO- • ONCNO short contacts exhibit the same 
magnetic character. The study was done on the nitronyl nitroxide crystals 
because they include many examples of crystals having dominant magnetic 
interactions, besides the fact that these crystal are the most interesting ones 
to analyze, because they include most of the purely organic magnets found 
up to now. 

For our statistical study of the geometry of the ONCNO- • ONCNO short 
contacts [24], we created a subset of nitronyl nitroxide crystals presenting 
dominant ferromagnetic interactions (the FM subset), and another separate 
subset with those crystals presenting dominant antiferromagnetic interactions 
(the AFM subset) [25]. These subsets were made by combining the data on 
nitronyl nitroxide crystals deposited in the Cambridge Crystallographic 
Database, with others provided by various authors. Of the initial set of 
crystals, we discarded those having R factors larger than 0.10 or having large 
molecular distortions. We also discarded co-crystals or crystals containing 
transition metals, and those crystals not showing dominant magnetic 
interactions in the 2-300 K range. We ended up with a database of 47 
crystals, 23 of them in the dominant FM subset, and the remaining 24 in the 
AFM subset. The structures of the molecules included in each subset, with 
their refcodes (when available, otherwise we have assigned to them an internal 
refcode which begins by a 0) are shown in Figs. 6 and 7. 

We then analyzed the relative orientation of the ONCNO- - -ONCNO short 
contacts using as parameters the six coordinates depicted in Fig. 8 . Other 
studies [26] have shown that the internal geometry of these groups is nearly 
invariant among the crystals of this family. Therefore, under these circum- 
stances, the six geometrical parameters of Fig. 8 are sufficient to fully describe 
the relative orientation of two ONCNO groups. 

The results from this statistical analysis [24] show that many commonly 
accepted magneto-structural correlations are not valid. The first one is that 
short NO- - -ON contacts are indicative of dominant antiferromagnetic 
interactions in the crystal, a behavior consistent with the McConnell-I 
predictions. As shown in Table 1, short NO- - -ON contacts are found in the 
ferro- and antiferromagnetic subsets and in a similar number. In particular, in 
the 3-4 A range of O- - -0 distances, the number of contacts in each subset, 
42% and 58%, is close to the 50% proportion that one should expect if the 
packing of the subsets of crystals would be identical. 

Another correlation that does not hold is that associating the magnetic 
properties with the relative disposition of the ONCNO groups. If this asso- 
ciation would work, there should be specific orientations of the ONCNO groups 
only present in the ferromagnetic subset, and similarly for the antiferromag- 
netic subset. In order to see such a failure, one has to include in the analysis 
the angular parameters (that is, Aj, A 2 , Ti, T 2 , and T 3 ). Figure 9 plots the 
position of the O 21 atom relative to the O 12 — N 12 — Ci — Nn — On atoms (which 
depends only on the D, Ai, and T 2 parameters) for the ONCNO- - -ONCNO 
contacts having an O- - -0 distance shorter than 5 A. It is clear from these plots 
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Fig. 6. Chemical structures of the R substituent for all the nitronyl nitroxide radicals whose 
crystals are included in the ferromagnetic subset 




that the distribution of the O 21 atom does not present regions in the FM subset 
excluded in the AFM subset, and vice versa. Instead, the distribution spreads 
over the whole range of values plotted, in a similar form. In fact, a cluster 
analysis of the sets of positions of the FM plus AFM subsets showed that the 
relative orientations of the FM subset are interlocked with those from the AFM 
subset, thus making impossible to distinguish them by the geometry of their 
contacts, that is, there are no specific orientations of the ONCNO groups 
associated to a type of magnetic interaction. The same conclusion was obtained 
when the range of 0- • 0 contacts included in the analysis was increased up 
to 10 A. 
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Fig. 7. Chemical structures of the R substituent for all the nitronyl nitroxide radicals whose 
crystals are included in the antiferromagnetic subset 



When a similar analysis packing analysis is carried out for the shortest 
C — H- • 0 — N contacts, the same conclusions are reached [24]: there are 
no statistically significant differences in the relative dispositions of the 
shortest C—H- ■ 0—N contacts in the FM and AFM subsets of crystals. 
Therefore, by combining the results for the O — N — CN — O- • O — N — C — N — O 
and the C — H- • O — N contacts, one can safely conclude that it is not 
possible to determine the nature of the dominant magnetic interaction in 
a crystal by looking at the geometry of only one of these two types of 
contacts. 

The impossibility to generate a structure- correlation relationship based on 
the geometrical properties of one class of contacts is indicative of one of the 
following options: 




The Mechanism of the Through-Space Magnetic interactions 



47 




D = O21 Oil 

= ^ 21 "' ^11 -^11 
^ 2=^21 - ^21 ■" ^11 







T, = f^n-On -0,,-N,, 

T2 - O21 ••• Oil ■ ^11 ■ 



T3-C2- N21 - 021 • 



o 



11 



Fig. 8. Definition of the coordinates used to define the relative position of two ONCNO 
fragments when the internal geometry of these fragments is frozen. D is a distance, Ai and 
A2 are angles, while Tj, T2 and T3 are dihedral angles 



(a) the McConnell-I is not working or is an oversimplification of the real 
mechanism, or 

(b) the magnetic nature of the interaction between adjacent dimers is 
determined by the combined effect of more than one class of contacts 
and, consequently, it is not possible to obtain a reasonable magneto- 
structural correlation by looking only at one of the components. 

The results presented in the previous sections indicate that the first option is 
certainly a valid one, but we will also explore if the second option is 
possible. 

There is a final point to be stressed about this statistical analysis: the 
similarity between the geometry of the ONCNO- • ONCNO contacts found in 
the FM and AFM subsets, also found for the C — H- • ONCNO contacts. This is 
clearly manifested in the results of Tables 2 and 3, in which we collect the 
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Table 1. List of ONCNO- • ONCNO contacts for crystals of the FM and AFM subsets within 
the range of distances indicated. Percentages of cases with intermolecular ferro- and 
antiferromagnetic interactions are also given 



FM subset 








AFM subset 




Distance 
range (A) 


Total number 
of contacts 


Number of 
contacts 


% 


Number of 
contacts 


% 


[0,3] 


0 


0 


0 


0 


0 


[0,4] 


24 


10 


42 


14 


58 


[0,5] 


92 


36 


39 


56 


61 


[0,6] 


204 


90 


44 


114 


56 


[0,7] 


378 


167 


44 


211 


56 


[0,8] 


608 


274 


45 


334 


55 


[0,9] 


901 


416 


46 


485 


54 


[0,10] 


1312 


611 


47 


701 


53 



FM 



AFM 





Fig. 9. Scatterplot of the values of the D, Ai and T 2 parameters for all ONCNO- • ONCNO 
contacts found in the FM and AFM subsets whose contacts lie between 3 and 5 A. The 
scattergram in the left plots the FM subset contacts, and that in the right the AFM subset 
contacts 



shortest contacts of each type found in the FM and AFM subsets. Clearly, there 
is a striking similarity in the distribution of distances and angles. This 
similarity is also found when one compares two crystals having the same range 
of distances. Thus, for instance, when we compare the two crystals of the FM 
and AFM subsets in which the shortest NO - -ON contacts are found, the 
LICMIT [27] and WILVIWIO [28] crystals, the geometry of the shortest 
ONCNO- - -ONCNO contacts is very similar. This is illustrated in Fig. 10, 
where we have plotted the position of the dimers in which the shortest 
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Table 2. Values of the geometrical parameters defining the ONCNO- • ONCNO contact for 
the five crystals of the FM and AFM subsets showing the shortest D distances. The distances 
are given in A and the angles in degrees. The refcodes of crystals in which these interactions 
are found are given in Figs. 6 and 7 



Refcode 


D 


-^1 


^2 


Ti 


T2 


T, 


FM subset 


LICMIT 


3.158 


146.8 


146.8 


180.0 


-0.1 


0.1 


ZORHIX 


3.168 


126.9 


71.5 


119.0 


46.5 


-76.2 


MMEPYC 


3.429 


153.2 


68.0 


173.3 


-40.1 


70.1 


OOOMPY 


3.499 


127.5 


76.4 


-162.1 


104.2 


-92.1 


PEYPUA 


3.719 


114.6 


59.0 


172.9 


-114.0 


-66.5 


AFM subset 


WILVIWIO 


3.159 


117.2 


134.1 


139.0 


134.3 


77.9 


5CL20H 


3.369 


82.6 


82.6 


180.0 


75.6 


-75.6 


OOOOAH 


3.431 


120.0 


75.8 


-38.9 


75.1 


62.7 


SUKBOP 


3.522 


77.7 


77.7 


180.0 


-66.7 


66.7 


ZIPTAT 


3.589 


79.2 


79.2 


180.0 


-82.7 


82.7 



Table 3. Values of the geometrical parameters defining the C(sp^) — H- • O — N contacts for 
the five crystals of the FM and AFM subsets showing the shortest D distances. The distances 
are given in A and the angles in degrees. The refcodes of crystals in which these contacts are 
found are also indicated 



Refcode 


D 


-^1 


-^2 


Ti 


FM subset 


HAFXOB 


2.339 


137.9 


167.2 


-145.3 


YOMYII 


2.382 


111.9 


162.2 


94.9 


OOOMPY 


2.419 


118.0 


168.6 


160.5 


ZORHIX 


2.582 


126.8 


148.1 


169.2 


OOGPNP 


2.600 


124.1 


132.1 


-150.9 


AFM subset 


OOOOAH 


2.295 


88.9 


158.0 


54.5 


YISCIM 


2.343 


162.6 


161.5 


-109.3 


SUKBIJ 


2.442 


142.1 


161.2 


-148.0 


OOOOBR 


2.447 


126.3 


149.6 


112.3 


YOXMED 


2.489 


131.3 


157.7 


-157.8 



ONCNO- • ONCNO contacts are found. Notice that there are two dimers in 
Fig. 10 for the WILVIWIO crystal, the reason being that there is a second short 
contact at an 0- • 0 distance of 3.384 A, besides the shortest 3.158 A contact. 
All the other short ONCNO- - -ONCNO contacts in the LICMIT and WILVIWIO 
crystals are located at distances larger than 4.3 A. The analysis of the dimers of 
Fig. 10 shows that the nearby ONCNO groups are oriented in a similar form in 
the two crystals, with these groups in a parallel arrangement: one of them 
displaced relative to the other in such a way that the shortest interfragment 
atom-atom distance between these ONCNO groups is a lateral 0- - -0 contact 
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3.384 A 



Fig. 10. Geometry of the dimers in which the shortest ONCNO- • ONCNO contacts are found 
in the LICMIT (leftmost dimer) and WILVIWIO (two dimers in the right) crystals. The value 
of the O- • O distance is also indicated 



(see Fig. 11 for tridimensional views of the positions of the ONCNO groups). 
The fact that in the LICMIT crystal the radicals of the dimer are in an up-down 
disposition while in the WILVIWIO crystal they are in an up-up disposition, 
is totally irrelevant from the point of view of the McConnell-I mechanism. 
Furthermore, not all the crystals of the AFM subset presenting short NO- • ON 
contacts have dimers packed in an up-up disposition. For example, the 
5CL20H AFM crystal presents a shortest NO- - -ON contact at 3.369 A (the 
second shortest in the AFM subsets) and the dimers follow an up-down 
disposition. According to the McConnell-I mechanism all these contacts 
should present magnetic interactions of the same character (the atoms making 
the shortest contacts have always the same sign for the atomic spin 
population). Therefore, such a similarity creates an interesting problem: 
how two crystals presenting similar ONCNO- - -ONCNO contacts can have 
different magnetic interactions? We will see in the next section how this puzzle 
can be solved. 
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4 

Ab initio Study of the Nature of the Through-Space Magnetic 
Interactions Present in Nitronyl Nitroxide Crystals: 

The LICMIT and WILVIW10 Case 

The previous sections have shown the weak points of many of the most 
common magneto-structural correlations used nowadays. Therefore, new and 
more sound magneto-structural correlations are required to explain the nature 
of the through-space interactions found in purely organic molecular crystals. 
These correlations should also serve as guiding principle in the design of new, 
purely organic molecular crystals. 

In this line of thought, we decided to investigate the nature of the magnetic 
interaction in cases where there are strong conflicts between the McConnell-I 
predictions and the experimental findings. As described in the previous 
section, one such conflicting case is found in the LICMIT [27] and WILVIWIO 
[28] crystals, that is, the a-phase of the 2,5-dihydrophenylnitronyl nitroxide 
crystal and the p-methylpyridylnitronyl nitroxide crystal. The first crystal is 
the member of the FM subset presenting the shortest N — O- • 0 — N contact 
(3.158 A), while the second presents the same property within the AFM 
subset, being the shortest N — 0- • O — N distance 3.159 A. We can assume that 
the shortest contact is sharing the same magnetic character as the crystal. 
Therefore, given the fact that the relative disposition of the ONCNO groups in 
these contacts is similar (Fig. 11), the obvious question is: why two similar 
ONCNO- • ONCNO contacts can have different magnetic nature? 

The LICMIT crystal (Fig. 12) orders ferromagnetically below 0.5 K, and 
shows an effective J/k^ value of 5-0.93 K [27]. The crystal belongs to the Pljn 
spatial group, with cell parameters a = 15.142 k, b - 12.320 A, c = 7.196 A, 
and = 99.18°. It has four molecules per unit cell, packed forming parallel 
strips along the c direction, in which the molecules are ordered in an up-up 
disposition (primary structure) [4d]. The strips form planes (secondary 
structure) [4d], ordered in a T-shaped disposition among themselves (tertiary 
structure) [4d]. The shortest NO- - -ON contacts (Figs. 10 and 11) are found 
between molecules of parallel strips, connecting radicals in one strip with these 
in the nearby parallel strip (remember that these molecules present an up- 
down disposition, previous section). 

The WILVIWIO crystal (Fig. 12) shows dominant antiferromagnetic 
interactions, its effective }/k^ constant being -74 K [28]. It belongs to the 
P-1 space group, with cell parameters a - 11.843 A, b = 12.695 A, 
c - 9.532 A, (X = 95.53°, ^ = 90.55°, and y - 146.89°. This crystal has two 
molecules per unit cell, packed forming planes along the b-c directions. Each 
plane is made by replicating a strip of molecules, ordered within the strip in an 
up-up disposition. The strips are the primary structure, and the plane is the 
secondary structure [4d]. These planes then pile up to form the crystal 
(tertiary structure) [4d]. The shortest NO- - -ON contacts (Fig. 10) are found 
within the planes, between adjacent molecules of the same strip. 

Given the exponential decay of the strength of the through-space magnetic 
NO- - -ON interactions, [29, 30] the dominant magnetic interactions in a crystal 
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Fig. 11. Frontal, top and side views of the geometry of the shortest ONCNO- • ONCNO 
contacts found in the LICMIT (upper) and WILVIWIO (middle and bottom) crystals. The 
value of the O- • O distance is also indicated 

can be expected to be those with the shortest ONCNO- • ONCNO contacts, in 
the case of the LICMIT and WILVIWIO crystals, the contacts shown in Figs. 10 
and 11. All the other ONCNO- - -ONCNO contacts in these two crystals lie at 
distances higher than 4.3 A and can be ignored in a first approach. 

If one applies the McConnell-I mechanism to the three contacts of Fig. 10, 
one would expect the three of them to be antiferromagnetic in nature, as the 
shortest contact is that made by the O atoms of nearby groups (in the 3.384 A 
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Fig. 12. a Structure of the LICMIT crystal, visualized along the c axis, b Structure of the 
WILVIWIO crystal, visualized along the a axis 



contact found in WILVIWIO one can think of an overlap involving also the N 
atoms, but this does not change the magnetic character of the interaction). 
Consequently, these contacts involve atoms having the same sign for the 
atomic spin population. Therefore, the McConnell-I mechanism will predict 
that these contacts are all antiferromagnetic and, because they are the 
dominant contacts, will predict these two crystals as presenting dominant 
antiferromagnetic interactions. 

Given the previous incorrect prediction, one has to conclude that the 
McConnell-I mechanism is oversimplifying the problem along one of the 
following two possibilities: 
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(a) the character of the magnetic interaction depends on the small geometrical 
changes found between the ONCNO- • ONCNO contacts of these three 
isomers (see Fig. 11), or 

(b) there are other contacts, besides the ONCNO- • ONCNO ones, that are 
important in defining the magnetic character of these dimers. 

We will explore these two possibilities with the help of ab initio computations 
on model systems. 

To discard the first possibility, we performed ab initio computations on 
(ONH — C — HN0)2 and (H 2 NO )2 model dimers, selected to reproduce the 
interaction of the isolated ONCNO groups of the Figs. 10 and 11. In the 
(H2N0)2 dimer we explored the magnetic character of the NO- - -ON contacts 
when the N — O group is isolated, while in the (ONH — C — HNO )2 dimer we 
explored the influence on the magnetic character of the NO- - - ON contacts 
of the conjugation with the other N — 0 group of the ONCNO groups. The 
geometry of the (ONH — C — HNO )2 dimers is that found in the LICMIT and 
WILVIWIO crystals for the shortest ONCNO- - -ONCNO groups. It was 
obtained by deleting all the atoms in Fig. 10 dimers except those in the 
two ONCNO groups, and then attaching an H atom to each NO group at the 
same distance it has in the H 2 NO molecule. The (H 2 NO )2 model dimer was 
obtained by also deleting the C and NO atoms not belonging to the NO group 
making the shortest contact. The magnetic character of the interaction in these 
model dimers can be obtained by computing the singlet-triplet energy 
difference in the dimer = E{S) — F(T)]. A positive value indicates a 

triplet ground state, that is, a ferromagnetic interaction, while a negative value 
is obtained when the singlet is the ground state and the interaction is 
antiferromagnetic. 

Previous computations have shown that one can obtain qualitative values of 
the singlet-triplet energy differences at the CASSCF level, [29, 31] although a 
precise computation of these values requires methods which include the 
dynamical electron correlation, like the CASPT2 methods. Alternatively, one 
can obtain quantitative the singlet-triplet energy differences in these dimers 
using dedicated energy difference methods [3 lb]. Table 4 collects the values of 
the AE^ ~ ^ values obtained for the two model dimers of the three dimers of 
Fig. 10. The CASSCF computations on the (H 2 NO )2 model dimer were done 
using a (2,2) complete active space, which produces the same qualitative 
singlet-triplet energy differences as the larger spaces used in Sect. 3.2. 
However, a larger (6,6) active space was used for the (ONH — C — HNO )2 model 
dimer, in order to describe the polarization effects found in each monomer. 
These effects are responsible for the presence of negative atomic spins in the 
central C atom of these monomers. The computations on the smaller dimer 
were done with three basis sets, the 6-31G, 6-31G(d), and 6-31 + G(d) basis, to 
explore the basis set truncation effects. On the larger dimer, we used only the 
smaller two basis sets for the same evaluation, since the addition of diffuse 
functions did not strongly affect the AE^ ~ ^ results in the (H 2 NO )2 dimer. 

A look at the results in Table 4 shows that, independently of the basis set 
and model dimer used, the singlet is the ground state in the two dimers, that is. 
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Table 4. Values of the AE^ ^ energy difference computed for the (ONH — C — HNO )2 and 
(H2N0)2 model dimers using the CASSCF and B3LYP broken symmetry methods (see text) 



Model system 


Crystal 


Basis set 


AE^ ~ ^ CASSCF cm“‘ 


AE^ ■ ^ B3LYP cm“ 


(H2N0)2 


LICMIT 


6-31G 


-3.4 


-4.6 






6-31G(d) 


-4.2 


-4.4 






6-31G -1- (d) 


-4.2 


-8.0 




WILVIWIO 


6-31G 










6-31G(d) 










6-31 -1- G(d) 






(ONH— C— HN0)2 


LICMIT 


6-31G 


-0.8 


-2.0 






6-31G(d) 


-1.0 


-1.8 




WILVIWIO 


6-31G 


- 2 . 2 Y-I 6 . 4 '’ 


-5.4Y-14.0'’ 






6-31G(d) 


-1.8Y-17.6'’ 


- 5 . 2 Y-I 5 . 0 '’ 



“ Dimer with a NO- ■ ON contact of 3.159 A in Fig. 7; 
Dimer with a NO- - -ON contact of 3.383 A in Fig. 7. 



the magnetic interaction in these dimers is always antiferromagnetic. The 
numbers also show a non-negligible influence of the second group in 
the singlet-triplet stability, which becomes a lot less antiferromagnetic in the 
larger dimer. A CASPT2 calculation on the CASSCF(6,6) wave function of 
the (ONH — C — HN0)2 model dimer confirmed the higher stability of the 
singlet state (by 2.28 cm“^ in the LICMIT geometry, and by 0.77 and 9.68 cm“^ 
in the WILVIWIO geometries). The CASPT2 values can be taken as a reference 
for other less accurate methods than the CASSCF (in this case, one finds that 
the CASSCF singlet-triplet separations have the right sign but are smaller than 
the CASPT2 numbers, although of the same order of magnitude). In 
conclusion, the results from these highly accurate ab initio methods confirm 
the McConnell-I prediction and indicate that looking only at the 
ONCNO- • -ONCNO groups it is not possible to reproduce the magnetic 
character of the dimers found in the LICMIT crystal. This fact justifies the lack 
of correlation between the magnetic character of the interactions and the 
relative geometry of the ONCNO groups reported in the previous statistical 
analysis of the nitronyl nitroxide crystals. 

Given the impossibility to reproduce the experimental character of the 
magnetic interaction by looking only at the ONCNO groups, one can suspect 
that the nature of the magnetic interaction is strongly affected by atoms whose 
atomic spin density is known to be small. To test at the ab initio level the 
validity of such a proposition, we computed the singlet-triplet energy 
difference (A£^ “ ^) using the whole geometry of the dimers of Fig. 10. Given 
the size of these molecules, the computations were done using a broken- 
symmetry formulation [32] of the unrestricted B3LYP determinant (hereafter 
identified as BS-UB3LYP), where by B3LYP we mean the usual three 
parameters non-local exchange and correlation functional [33]. Such a 
formulation has been shown to provide a very precise description of the 
singlet-triplet gaps in organometallic complexes [34]. We have tested 
the quality of the BS-UB3LYP computations on our dimers by recomputing 
the singlet-triplet gap of Table 4 for the (ONH — C — HNO )2 and (H 2 NO )2 
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dimers (last column in Table 4). The results show that the BS-UB3LYP singlet- 
triplet splitting has always the same sign and order of magnitude as the 
splitting computed at the CASSCF level, although the BS-UB3LYP values are 
always a bit larger (in fact, close to the CASPT2 values). Therefore, we can feel 
confident about the ability of the BS-UB3LYP method to compute the nature of 
the through-space magnetic interactions. 

The BS-UB3LYP computations on the whole LICMIT dimer gave a AE^ ~ ^ 
value of 1.3 cm“^ [34] (to be compared with an experimental value of 
1.29 cm“^), while for the two WILVIWIO dimers the values were -9.1 cm“^ for 
the shortest one, and -28.2 cm“\ for the largest one (to be compared with an 
experimental averaged value of -102.9 cm“^) [34]. Therefore, when the whole 
geometry of the dimer is included, the ab initio methods can reproduce the 
experimental character of the magnetic interaction in the LICMIT and 
WILVIWIO crystals. This is a numerical proof that the atoms of the radical 
which hold in total only 10% of the atomic spin density play an important role 
in determining the nature of the magnetic interaction. An inspection of the 
geometry of the LICMIT dimer and the WILVIWIO dimers suggests that 
the NO- • aromatic ring contacts could be the key interaction in establishing 
the ferromagnetism in the LICMIT dimer; these are left out when only the 
ONCNO groups are considered. However, more detailed studies are needed 
before this fact can be stated without discussion. In any case, the ab initio 
computations on the whole dimer have shown that the magnetic character of 
the interaction is not a consequence of the relative geometry of only one 
functional group, but depends on the relative positions of all the groups. This 
relative position of many functional groups receives the name pattern. 
Therefore, we can talk about magnetic patterns, that is, the patterns found in 
magnetic crystals, each of them having a magnetic net interaction, not 
necessarily of the same nature in all the patterns of the same crystal. The 
nature of the pattern is not solely defined by the relative geometry of the 
ONCNO groups. Consequently, although the orientations of these groups in 
various patterns can be similar, the magnetic character of patterns sharing the 
same ONCNO relative orientations is not necessarily the same. This explains 
the apparent inconsistency found in the statistical analysis of the nitronyl 
nitroxide crystals presenting dominant magnetic interactions, in particular, 
the lack of dependence of the magnetic character with the relative orientation 
of the ONCNO groups. 

To complete our evaluation of the character of the magnetic interactions 
present in the LICMIT crystal, the only point that remains to be tested is our 
hypotheses that the dominant magnetic character of the crystal was associated 
to the shortest contact. Therefore, we computed the singlet-triplet energy 
difference for all the nearest neighbors that one of the radicals makes with the 
surrounding ones, irrespective of the type of contacts found between these 
radicals and their shortest distance. We found [35] five different types of 
nearest neighbor dimers in the LICMIT crystal, plotted in Fig. 13. They present 
NO- - -ON distances of 3.158, 4.594, 5.525, 5.856, and 6.294 A, the shortest 
one being that for the dimer of Fig. 10. The singlet-triplet energy difference in 
each of these dimers are 1.3, 0.15, 0.04, -0.15, and 0.11 cm“\ respectively [35]. 
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Therefore, the strongest magnetic interaction is hy far the shortest one (by one 
order of magnitude) and will dominate over the rest of radical- • radical 
interactions. It is interesting to note the ferromagnetic character of all but one 
of the nearest neighbors radical- - -radical interactions in the LICMIT crystal. 
The same applies to the WILVIWIO crystal [35]. Consequently, the dominant 
magnetic assumption made in the analysis of the LICMIT and WILVIWIO 
crystals seems to be a valid one. 

5 

Concluding Remarks 

From the previous studies we can extract some useful conclusions for a better 
understanding of the through-space magnetism found in purely organic 
molecular crystals. First of all, we have seen that the McConnell-I mechanism, 
the most popular used nowadays, lacks of a consistent methodological basis 
and part of its success in the [2,2]paracyclophanes relies in a fortuitous error 
cancellation due to the high symmetry of these biradicals. The McConnell-I 
predictions are shown to fail in the H 2 NO- - -ONH 2 dimers, and also in some 
orientations of the methyl- - -allyl dimer. We have also presented one example 
of a nitronyl nitroxide crystal in which McConnell-I fails, and discussed the 
problems found when using this mechanism on nitronyl nitroxides presenting 
packings driven by hydrogen bonds. 




Fig. 13 . Different types of nearest neighbor dimers found around a radical in the LICMIT 
crystal. The O- - -O distance is given for each dimer 
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A detailed statistical analysis of the packing of many nitronyl nitroxide 
crystals presenting dominant magnetic interactions was also done. No 
preliminary assumptions were made on which geometrical parameters were 
going to play a determinant role in defining the magnetic character. This 
analysis shows that the magnetic character of the crystals does not present any 
relationship with the relative orientation of the individual ONCNO- • ONCNO 
or C — H- • ONCNO contacts. Therefore, it suggested that some kind of 
cooperative effect among all the functional groups present in the molecule is at 
work in these crystals. Ab initio studies on the LICMIT and WILVIWIO 
crystals provides a numerical proof of existence on these types of cooperative 
effect: while one cannot explain the ferromagnetic character of the LICMIT 
crystal looking only at the relative orientations of the ONCNO groups, the 
ferromagnetism is found when the six-membered rings of the radicals are 
included in the computation. This proves that the character of these magnetic 
interactions in the dimer is determined by all groups present in the dimer, 
regardless of the atomic spin population of their constituting atoms. That is, 
we have to look at the magnetic patterns (combination of short contacts 
defining a energetically stable supramolecular aggregate of the magnetic 
crystal). Each presents a dominant magnetic interaction, resulting from the 
combination of the magnetic interactions generated by the individual contacts, 
which depend on all the functional groups present in the pattern. In the 
LICMIT dimer, it seems that the NO- • aromatic ring contact works in 
conjunction with the ONCNO- - -ONCNO contacts in establishing the magnetic 
character of the dimer, but further numerical tests are needed to fully prove 
this point. 

As a consequence of that cooperativity, one cannot expect to carry out 
reliable magneto-structural analysis of purely organic molecular crystals 
using oversimplified empirical magneto-structural correlations which as- 
sume that only one group is determining the magnetic character of the 
dimer. This fact discards the oversimplified analysis performed using the 
McConnell-I mechanism in which the position of the ONCNO groups are 
analyzed. We have also seen that McConnell-I also fails when other groups 
are considered, in part because of the simplistic form in which it is 
generally used (the exchange integrals are not considered). Further studies 
are under way to be able to define reliable magneto-structural correlations 
based on the correct assumptions. At the present moment, these studies 
[35] seem to indicate that one should use a more elaborate mechanism, like 
some sort of improved version of the McConnell-II mechanism, capable of 
solving the problems detected by Kahn and others in some well known 
magnets, but along the avenues described by Kinoshita in the first chapter 
of this book [36]. 
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The nitroxide radical is one example of a stable organic radical, and the oxygen atom of the 
nitroxide group exhibits weak basicity; therefore, it demonstrates coordinating ability with 
respect to transition metal ions. There are many ^-conjugated oligo-nitroxides which have 
high-spin ground states. When the oligo-nitroxides are used as a bridging ligand for 
paramagnetic metal ions, we can control both of the structure of the complex and the spin 
structure. This chapter discusses one-dimensional, two-dimensional, and three-dimensional 
magnetic materials derived from coordination complexes consisting of magnetic metal ions 
and nitroxide radicals. 

Keywords: Organic radicals. Magnetism, Transition metal complexes 



1 Introduction 62 

2 Nitroxide Radicals 62 

2.1 Magnetic Interaction Between the Transition Metal Ions 

and Nitroxide Radicals Through Coordination Bonds 64 

2.2 Design Strategy for Nitroxide Radical-Metal 

Magnetic Complexes 68 

2.3 Preparation of Metal-Poly(Nitroxide) Radical Complexes 69 

3 One-Dimensional Nitroxide-Metal Systems 73 

3.1 Structure and Magnetic Properties of Ferrimagnetic 1-D 

Chain Formed by Manganese(II) and Nitroxide Radicals 75 

3.1.1 Manganese(II) Complexes with Nitronyl Nitroxides 75 

3.1.2 Manganese(II) Complexes with Triplet Bis-Nitroxide Radicals . . 75 

3.2 Conclusion 77 

4 Two-Dimensional Metal-Nitroxide Systems 78 

4.1 Structure and Magnetic Properties of Ferrimagnetic 2-D 

Sheet Formed by Manganese(II) and Nitroxide Radicals 80 

4.1.1 2-D Manganese(II) Complexes with Nitronyl Nitroxides 80 

4.1.2 Manganese(II) Complexes with Triplet Bis-Nitroxide Radicals . . 81 



Dedicated to the memory of Prof. Olivier Kahn 



Structure and Bonding, Vol. 100 
© Springer-Verlag Berlin Heidelberg 2001 




62 



K. Inoue 



4.2 Conclusion 84 

5 Three-Dimensional Metal-Nitroxide Systems 86 

5.1 Crystal and Molecular Structures 86 

5.2 Magnetic Properties of 3-D Systems 86 

5.3 Conclusion 88 

6 Conclusions 89 

7 References 89 



1 

Introduction 

The construction of molecular-based magnets that have a well-defined 
magnetic structure is a scientific topic of increasing interest. The idea is to 
establish unprecedented macroscopic spins of long-range order in molecular 
systems. There are mainly three steps to generate molecule-based magnets: 

i) align the magnetic moments in the molecule; 

ii) arrange the molecules in the crystals; and 

iii) order all the spins in a mesoscopic scale. 

In approaches employing purely organic crystals there are many difficulties 
in steps ii) and iii). A network structure may not be a problem if an 
appropriate crystal design is made. However, since the exchange coupling 
between neighboring molecules through van der Waals forces, hydrogen 
bonding, hydrophobic interaction, etc., is not necessarily strong, it is difficult 
to expect strong intermolecular magnetic coupling. There are principally two 
conditions to obtain high-Tc magnets: 

1) make strong exchange interactions between spins and 

2) all the spins are ordered in two- or three-dimensions. 

In order to satisty the above conditions, it is necessary to create polymeric 
networks among spins. It is, therefore, a powerful design strategy to make 
extended polymeric structures by assembling free radicals by means of magnetic 
metal ions. Such polymeric complexes satisfy the above-mentioned conditions 
(i-iii). In this chapter magnetic materials made up of coordination complexes of 
magnetic metal ions and nitroxide radicals are discussed. 

2 

Nitroxide Radicals 

In order to use free radicals as ligands for a magnetic metal, it is necessary for 
the radicals to have sufficient coordinating ability [4-9]. The oxygen atom of 
the nitroxide group has a weak basicity; and thus has coordinating ability to 
transition metal ions. 
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There are some Ti-conjugated oligoradicals that have a high-spin ground 
state. Various his-nitroxide radicals with triplet ground states and tris- 
nitroxide radicals with quartet ground states have been designed on the basis 
of the m-phenylene topology and prepared by oxidation of the corresponding 
hydroxyamines. The magnetic interaction between nitroxides depends on the 
distances. High-spin oligo-nitroxide radicals that have been used as ligands for 
magnetic metal ions are summarized in Scheme 1 [10-18]. 
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2.1 

Magnetic Interaction Between Transition Metal Ions 
and Nitroxide Radicals Through Coordination Bonds 

Syntheses of metal-nitroxide complexes were first reported in the early 1970s 
as spin labels for metal complexes or catalysts for living polymerizations. In 
principle, when the overlap between the magnetic orbitals is large, a normal 
covalent bond will be formed and the spins are paired. In practice, however, 
the bond is only partial as judged from the observed pairing energy of the 
electrons. On the other hand, when the overlap is small or zero, direct 
exchange will be operative, giving either weak antiferro- or ferromagnetic 
coupling. The coupling modes are classified in terms of the geometry of the 
two magnetic orbitals involved. 

Structures of metal-nitroxide complexes are classified in the first approxi- 
mation into two limiting arrangements, the nitroxide oxygen occupying an axial 
or an equatorial coordination site. In both limits, the magnetic orbital of 
copper(II) may be loosely described as xy, lying in the equatorial plane of the 
pyramid with the lobes pointing towards the ligands. When the nitroxide ligand 
occupies an axial position, its magnetic orbital is essentially orthogonal to the 
magnetic orbital of copper(II) and the coupling is expected to be ferromagnetic. 
Indeed, this has been found to be the case for a number of complexes in which 
magnetic susceptibility data showed the existence of ferromagnetic coupling in 
the range 30-100 K [19-23]. On the other hand, when the nitroxide oxygen atom 
is in the equatorial (or basal) position in copper(II) complexes, the exchange 
interaction tends to be antiferromagnetic [19-21]. 

In the case of high spin manganese(II), at least one of five magnetic 
orbitals must have the correct symmetry for overlapping with those of the 
nitroxides. As a result, the interaction is antiferromagnetic and a ground 
S = 3/2 (=5/2 - 1/2 - 1/2) state will result in complexes of the general formula 
[Mn(hfac) 2 (nitroxide) 2 ]. For TEMPO and PROXY ligands, J/k^ values are -227 
and -302 K, respectively [19-21, 24, 25]. 

3-Imidazoline nitroxides and Ullman’s nitronyl nitroxides display have 
basic centers and, therefore, serve as good bridging ligands for extended 
structures with coordinatively doubly unsaturated metal ions. These species 
will be discussed in Sect. 3 and 4. At the same time, since Cu(II) has the ability 
to form complexes with coordination numbers of both 5 and 6, trinuclear 
complexes with a 3:2 stoichiometry [3Cu(hfac)2 • 2nitroxide] are obtained 
[26, 27]. The magnetic coupling in these complexes is ferromagnetic except for 
one with the 2-ethyl derivative of nitronyl nitroxide [26]. In general, exchange 
coupling of Mn(II), Ni(II), and Co(II) with nitronyl nitroxide ligand is 
antiferromagnetic in character. 

There are cyclic transition metal complexes with nitroxide radicals or 
nitronyl nitroxide radicals. The number of spins is finite for these complexes; 
therefore, the systems can afford a rigorous solution of the spin Hamiltonian 
to give the sign and magnitude of the exchange coupling. These values serve 
as good measures for designing and analyzing of new extended systems 
containing two or more types of exchange coupling parameters. 
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The cyclic dimer complex of bis(hexafluoroacetylacetonato)manganese(II) 
with 5-tert-butyl-l,3-phenylenebis(N-ferf-butylnitroxide) [28, 29] is such an 
example. The systems have been analyzed in order to obtain two sets of two 
intramolecular exchange coupling parameters Ji and J 2 (Scheme 2). The crystal 
structure (Fig. 1) of the first complex, [Mn(hfac )2 • BNO(_bu] 2 > is found to be a 
cyclic analog of 1-D complex polymers made of similar components lacking a 
ring system, a tert-butyl moiety, or with halogen substituents in its place (see 
Sect. 3) [30-33]. The molecule consists of two Mn(II) ions (S — 5/2) and four 
nitroxide radicals (S = 1/2) arranged in a cyclic array as shown in Scheme 2. 
The molar paramagnetic susceptibility of the complex was investigated in 
the temperature range 2-300 K. The Xmoi^’ values exhibit a continuous increase 
from room temperature to ca 50 K, reaching a plateau at 5.71 K cm^ mol“\ 










Scheme 2. 




Fig. 1. Cluster structure of [Mn(hfac )2 ■ BNO,.bu] 2 - Fluorine atoms, N-ferf-butyl groups and 
hydrogen atoms are omitted. Thermal ellipsoids are drawn at the 50% probability level 
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This agrees well with the theoretical spin-only value for an S = 3 spin ground 
state, which can he explained by ferromagnetic (F) interaction between two 
S = 3/2 pseudo-spins (5/2 - 2 x 1/2), obtained by assuming a predominant 
antiferromagnetic (AF) coupling {J 2 ) between the manganese(II) ions and the 
two coordinated radicals. This high-spin ground state is stabilized with respect 
to upper energy levels by several hundred wave numbers, as deduced from the 
decrease in Xmoi^’ with increasing temperature. A decrease in Xmoi^’ at the 
lowest temperatures may be analyzed in terms of intermolecular AF coupling. 

The XmoiT vs. T plots were analyzed in two ways: one is a rigorous solution 
of the spin Hamiltonian for the six-spin system [34-36], and the other by a 
model of the coupling of a cyclic array of effective spins S = 3/2. The first 
method together with the molecular field contribution taken into account by 
the Curie-Weiss law gives h/k^ = -1-1005 K (±11 K), / 2 /kB = -539 K (±19 K) 
and Weiss constant A = -0.10 K (±0.01 K), with g = 1.953 [37], provides a 
very good description of the experimental data over the entire temperature 
range. The second method gave J/k^ = ±68.6 K (±0.6 K) for repeating S = 3/2 
units and the Lande factor g = 2.002. Simulation of the observed data is quite 
good except that, the high-temperature data are poorly reproduced by the 
second model. This is due to the lack of precision in the approximation 
method in locating correctly the higher energy levels where spins are more 
populated at higher temperature. 

A zero-dimensional complex of bis(hexafluoroacetylacetonato)manga- 
nese(II) and 4-N-fert-butylnitroxide-4'-(l-oxyl-3-oxide-4,4,5,5-tetramethylim- 
idazolin-2-yl)biphenyl (BIPNNNO) [38] is also such an example. The complex 
consists of discrete molecules of [Mn(hfac) 2 ] • (BIPNNNO )2 and has a triclinic 
structure (space group PI) where a - 11.267(3) A, b - 12.605(2) A, c = 
10.661(3) A, a = 99.16(1)°, P = 96.97(3)°, y = 82.60(1)°, V = 1471.5(6) A^ 
The [Mn(hfac) 2 ] • (BIPNNNO )2 molecule is shown in Fig. 2. The Mn(II) ion 
lies on a crystallographic inversion center and is coordinated by a pair of 
BIPNNNO moieties and a pair of bidentate hfac anions, which afford an MnOe 
octahedral environment. The two BIPNNNO molecules coordinate trans to one 
another. The Mn — O bond distances are 2.128(5) A for the BIPNNNO radical 
and 2.153(5) and 2.130(5) A for the hfac ligands. The formation of a 1:2 
complex, [Mn(hfac) 2 ] • (BIPNNNO) 2 , suggests the higher coordination ability 
for N-terf-butyl nitroxide than for nitronyl nitroxide to the Lewis acidic 
manganese ion. The temperature dependence of Xmo\T of the complex is 
illustrated in Fig. 3. The low temperature value of Xmoi^ tends toward 
0.375 emu/mol, a limit for the S = 1/2 ground state. This value is expected for 
the negative Mn — NO and positive intraradical NO — NN exchange parameters. 
The complex remains paramagnetic down to 1.8 K. As demonstrated in Fig. 3, 
the magnetization curve at this temperature can be well described by a 
Brillouin function for S = 1/2, Bi/ 2 {H/T), neglecting the intercluster exchange 
interaction. 

The high temperature limit expected for separately-acting four spins 1/2 
and one spin 5/2, 5.875 emu/mol, is far from being reached. The character of 
XmoiT temperature dependence indicates that both of the exchange interactions 
are comparable and that new levels are successively excited with increasing 
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Fig. 2. View of the Mn(hfac)2(BIPNNNO)2 molecule 




Fig. 3. Temperature dependence of XmoiT of the complex Mn(hfac)2(BIPNNNO)2. Experi- 
mental data are denoted by symbols. Solid line corresponds to the theoretical fit (see text) 
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Fig. 4. Magnetic structure of Mn(hfac)2(BIPNNNO)2 



temperature T up to room temperature. The temperature variation of XmoiT 
was simulated numerically in the Heisenberg approximation with use of the 
Van Vleck equation. The energy level scheme was determined by exact 
diagonalization of the following Hamiltonian (Fig. 4): 

H = — 2 /[SiS 2 + 7(6263 + 6364) + 6465] ( 1 ) 

where / and yj are the NO — NN and Mn — NO exchange integrals, respectively, 
and S 3 is the Mn^"^ spin operator. The best fit parameters are J/k^ = 72 ± 5 K 
and 7 = -1.9 (yj/k^ = -135 ± 10 K) with the purity factor P = 1.1. 



2.2 

Design Strategy for Nitroxide Radical-Metal Magnetic Complexes 

As we have previously seen, a set of two conditions makes the necessary and 
sufficient conditions for realizing metal-radical complexes having interesting 
magnetic properties. 

1) A ligand must possess at least two nitroxide coordination sites. Any 
nitroxide group may form a metal complex having this nitroxide as a ligand, 
but there is no chance of forming an extended structure without additional 
coordination sites. 

2) When high-spin nitroxide radicals are utilized as ligands for the complex, 
the interaction among nitroxides must be ferromagnetic. 

These are the reasons why nitronyl nitroxides and high-spin oligo-nitroxide 
radicals are of importance for material design. When an organic free radical 
ligand has two ligation sites, e.g., nitronyl nitroxides or bis-nitroxide radicals, 
extended complex structures that are often chain polymers or macrocycles, are 
formed with coordinatively doubly unsaturated metal ions. Extension of this 
design strategy to tris- nitroxide radicals exhibiting quartet ground states has 
led to the construction of two- (2-D) and three-dimensional (3-D) network 
structures in which both organic 2p and metallic 3d spins have been ordered in 
macroscopic scales [1, 3, 39-41]. This schematic is shown in Fig. 5. 

TT-Conjugated oligo-nitroxide radicals (Scheme 1) have been employed as 
bridging ligands in which the spins of unpaired electrons interact ferromag- 
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Fig. 5a-c. Schematic drawing of the coordination networks of the high-spin oligo-nitroxide 
radical and manganese systems 



netically {/(intraligand) > 0}. Dimensions of the complexes, as well as sign and 
magnitude of the exchange coupling between the adjacent spins, may be 
readily tuned in this methodology [1, 3, 39-41]. 

The tris-monodentate triradicals TNO, TNOPB, TNOPA, TNOB, with 
quartet ground states (S = 3/2) in which the radical centers are arranged in 
a triangular disposition, would form 3:2 complexes with a coordinatively 
doubly unsaturated 3d metal ion M. In an ideal case a 2-D hexagonal network 
structure would be generated (Fig. 5b). A T-shaped quartet triradical carrying 
two non-equivalent ligating sites, e.g., TNOP, would form a 1-D chain by using 
two terminal nitroxide groups. The middle nitroxide group might then be used 
to cross-link the chains to form a 3-D network structure (Fig. 5c) depending 
on whether the second bridging takes place between the same chains as cross- 
linked by the first bridging. The spin alignment in these systems would be very 
much stabilized and is expected to give higher- Tc magnets. 



2.3 

Preparation of Metal-Poly(Nitroxide) Radical Complexes 

The complex formation is a kinetically and/or thermodynamically controlled 
self assemblage of the reactants. Typical procedures are as follows. A 
suspension of manganese(II) bis(hexafluoroacetylacetonate) dihydrate, 
[Mn(hfac )2 • 2 H 2 O], in n-heptane is heated at reflux to remove water of 
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hydration by azeotropic distillation. To the resulting cooled solution is added 
BNOh in n -heptane ether. The mixture is then concentrated under reduced 
pressure and the concentrated solution is allowed to stand to furnish black 
needles of [Mn(hfac )2 • BNOh] from a deep brown solution. It is preferable to 
carry out the reaction in an inert atmosphere and under anhydrous conditions, 
sometimes in a refrigerator. The reaction is completed by precipitation. Some 
complexes can be recrystallized but others dissociate in solution. 

Excess of either component can generate complexes of different composi- 
tions. For example, the reaction of [Mn(hfac) 2 ] with tris(nitroxide) TNOP is 
complex while an equimolar mixture in ether containing n-hexane at -10 °C 
provides black blocks of a 1:1 complex [{Mn(hfac) 2 | • TNOP • n-C 6 Hi 4 ], a 
mixture containing a 1.7 molar excess of [Mn(hfac) 2 ] in n-heptane-ether 
generates black blocks of a 3:2 complex [{Mn(hfac) 2}3 • TNOP 2 ] in ten days at 
0 °C. The complex [{Mn(hfac) 2}3 • TNOPB 2 • n-CyHig] is obtained by dissolv- 
ing [Mn(hfac )2 • 2 H 2 O] in a mixture of diethyl ether, n-heptane, and benzene, 
followed by addition of TNOPB in benzene. Black blocks are formed from a 
deep violet solution. While [Mn(hfac) 2 l gave similar black violet 3:2 complexes 
with tris(nitroxide) TNOB, but TNO and bis(nitroxide) PNNNO did not form 
polymeric complexes, probably due to steric hindrance around the ligand 
molecule. 

Recently, a notorious side reaction has been elucidated leading to undesired 
by-products that have unique [3 -t 3] benzene-dimer structures [42]. 1-D 
ferrimagnetic complexes, [Mn(hfac )2 • BNOx]n (X = Cl or Br), are typically 
obtained by the reaction of Mn(hfac )2 with BNOx- When crystallization 
requires a period of several days, the black solutions often turn yellow in about 
one day and do not afford the expected black polymer complexes. Instead, 
yellow crystalline precipitates are obtained under these conditions (Scheme 3). 

An X-ray structure analysis revealed that BBNOx • [Mn(hfac )2 • 
H20]2 ■ CH 2 CI 2 has a [3-1-3] benzene-dimer structure [X = Cl or Br; 
BBNOx = 3, lO-dihalo-5, 8, 11, 12-tetrakis(N-tert-butylimino)-tricyclo[5, 3, 1, 
1^’®] dodeca-3, 9-diene N,N',N",N"'-tetraoxide} (Fig. 6). Two crystallographi- 
cally equivalent manganese(II) ions have an octahedral coordination and are 
coordinated with four oxygen atoms of two hfac ligands, one oxygen atom of 
water, and one oxygen atom, 01, of BENObi- Since the dimer complex is 
chiral, both enantiomers are contained within each unit cell (Scheme 4). The 
BBNOci complex is isostructural to the bromine derivative. 

The resonance structure BNObi' must be responsible for the reaction 
leading to the dimer complex (Scheme 5); either dimerization of BNObi' or 
attack of BNObi' to free and complexed BNObi- Whereas the 1-D ferrimag- 
netic complex [Mn(hfac )2 • BNOBr]« is a kinetic product and precipitates 
from solution at an early stage of reaction, the yellow crystals of 
BBNObi ■ [Mn(hfac )2 • H20]2 ' CH 2 CI 2 appear to be a thermodynamic product. 

Another interesting feature of this work is the liberation of the dimer ligand 
BBNObi) free from manganese ions, by dissolving the complex in ether. Two 
water molecules of hydration appear to be crucial for the stability of 
BBNObi ■ 2 H 2 O; it is stable in water at 100 °C but starts to decompose by 
dissociation even at -78 °C in CH 2 CI 2 when dehydrated by molecular sieves. 




Metal-Aminoxyl-Based Molecular Magnets 71 




Mn(hfac)2*H20 

d 



Mn(hfac)2*H20 

BBN0x*[Mn(hfac)2*H20]2 BBNOx 

Scheme 3. 




Mn(hfac)2*H20 Mn(hfac)2*H20 




Mn(hfac)2*H20 Mn(hfac)2«H20 



-Mn(hfac)2 0.-|^>^A|^-.0 
f-Bu f-Bu f-Bu 




Scheme 4. 
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3 

One-Dimensional Nitroxide-Metal Systems 

The first one-dimensional polymeric transition metal complex with nitroxides 
was [Cu(hfac )2 • TEMPOL] [24]. Coupling between the nitroxide radical and 
copper(II) is ferromagnetic with 2Jlk^ = 19 ± 7 K and this (S = 1) pair 
couples with adjacent pairs through quite lengthy superexchange paths of the 
cr-bonds of TEMPOL by 2J' = —78(2) mK. Since Gatteschi et al. reported a 
polymeric transition metal complex with nitronyl nitroxide (NIT) in 1986 
[19, 43-56] others with l,3-phenylenebis(nitroxide) derivatives (BNOr) [1, 15, 
30-33, 41, 57-61] have been documented (Table 1). Bis-monodentate nitronyl 
nitroxides with doublet states (S = 1/2) and bis-monodentate 1,3-phenylene- 
bis(nitroxides) with triplet (S = 1) ground states form 1:1 complexes having 
one-dimensional infinite chain structures with coordinatively doubly unsat- 
urated paramagnetic metal ions. The coupling results in the alignment of 
neighboring spins in either parallel or antiparallel fashion to one another along 
the chain thereby giving one-dimensional ferromagnets or ferrimagnets, 
unless the sizes of the spins match each another in the latter. 

The paramagnetic susceptibility values of both kinds of materials are 
expected to diverge at low temperatures, as a result of lengthening of the 
correlation of the spins along the chain. Depending on the nature of the 




Table 1. Nitroxide-metal-based one-dimensional complexes 
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Data were evaluated by quantum-classical approximation of ferromagnetic 3/2 spins. /moNO /irTr define the interactions between Mn and NO and between 
two NO— Mn— NO’s groups exhibiting S = 3/2 spins, respectively. 
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additional interchain interaction, the chain polymers become a bulk antifer- 
romagnet or a ferri/ferromagnet. In these magnets the interchain magnetic 
interactions are much weaker than the intrachain interaction and, therefore, the 
transition temperatures to three-dimensionally ordered states are relatively low. 



3.1 

Structure and Magnetic Properties of Ferrimagnetic 1-D Chain Formed 
by Manganese(ll) and Nitroxide Radicals 



3.1.1 

Manganese(ll) Complexes with Nitronyl Nitroxides 

Manganese(II) ions form crystalline chain complexes with various NITr 
(R = i-Pr, Et, w-Pr, P-alkyloxyphenyl) in which Mn(II) is hexacoordinated 
with four oxygen atoms of two hfac molecules and with two oxygen atoms of 
two different NIT radicals. The other oxygen atoms of the two NITr radicals 
are coordinated to the adjacent Mn(II) ions (Scheme 6). All the above 
complexes are ferrimagnetic and order at ca. 8 K. Quantitative analyses of the 
magnetic susceptibility data are performed utilizing a model of classical- 
quantum chains. The results are summarized in Table 1. 

[Mn(hfac) 2 ] gave a crystalline 1-D chain complex with the NITchirai radical. 
The X-ray crystal structures of the NITchirai complex solved in the chiral space 
group P2i2i2i (No. 19). The asymmetric unit of [NITchirai ' Mn(II)(hfac) 2 ]„ 
consists of one Mn(II) ion, two hfac anions, and one chiral radical NITchirai- The 
Mn(II) ion and the NITchirai form a one-dimensional structure along the crystal 
fl-axis. The saturation magnetization of this complex was 3.8 /Tb, which agrees 
well with the theoretical limit, assuming antiferromagnetic coupling between 
the manganese and radical spin. This complex is ferrimagnetic and orders at 
4.5 K [16]. 



3.1.2 

Manganese(ll) Complexes with Triplet Bis-Nitroxide Radicals 



Manganese(II) ions couple form crystalline chain complexes with triplet bis- 
nitroxide benzene derivatives [30-33]. X-ray crystal structures of the achiral 
BNOr complexes solved in the monoclinic P2i/n space group (No. 14) with 
Z - 4 reveal that the manganese(II) ions have an octahedral coordination with 
four oxygen atoms of two hfac anions and two oxygen atoms of two different 
BNOr molecules. The latter is bound to the Mn(II) ion in a cis-configuration. 



(hfac )2 

R 




R = ;-Pr, Et, n-Pr 



Scheme 6. 
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As a result, the Mn ions and biradical molecules form a helical 1-D polymeric 
chain structure along the crystal b axis. All the hexacoordinated Mn(II) ions 
have either the A- or the A-configuration along a given chain. Two tert- 
butylnitroxide groups are rotated out of the phenylene ring plane in a 
conrotatory manner, but with different angles; each BNOr molecule in the 
crystal has no symmetry element and is therefore chiral, i.e., R or S. The 1-D 
polymeric chains are therefore isotactic as all units of the same chirality reside 
on a given chain (Fig. 7a). The crystal lattice is, as a whole, achiral since an 
enantiomeric chain is present. In the case of the chiral BNOchirai ligand, the 
crystal lattice is, as a whole, chiral (Fig. 7b) [15]. 
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The observed Xmoi^’ value at 300 K varies between 2.2 and 2.4 emu K mol~S 
which is slightly larger than the theoretical limit of 1.88 emu K moP^ expected 
for S = 3/2 ferromagnetic chain compounds. As T is lowered, the XmoiT value 
increases monotonically in proportion to the increase in the correlation 
length within the network. Together with the lack of a minimum at lower 
temperature, the room temperature Xmoi^’ value also points to the operation of 
strong (more negative than -300 K) antiferromagnetic coupling between the 
Mn(II) ion and the nitroxide radical of BNOr, in which the onset of the 
intramolecular ferromagnetic coupling is meaningful. 

The saturation magnetization values, Ms = 3 of all 1-D of 

Mn(hfac) 2 (BNOR) complexes agree well with the theoretical limit, assuming 
antiferromagnetic coupling between the manganese and triplet radical spins. 

The low-field susceptibilities of the Mn(hfac) 2 (BNOH) and 
Mn(hfac) 2 (BNOp) complexes show a sharp cusp at 5.5 K and 4.5 K respectively 
(Fig. 8). This behavior indicates that these complexes possess an antiferro- 
magnetic ground state. On the other hand, the Mn(hfac) 2 (BNOci) and 
Mn(hfac) 2 (BNOBr) show demonstrate ferrimagnetic behavior (Fig. 8). 

Magnetization at 1.8 K revealed metamagnetic behavior for the 
Mn(hfac) 2 (BNOH) and Mn(hfac) 2 (BNOp) complexes (Fig. 9). Namely, although 
the response of the magnetization was not sensitive to the weak applied magnetic 
field below Hq, a behavior characteristic of an antiferromagnet, a sharp rise and 
approach to saturation of magnetization, characteristic of a ferromagnet, was 
observed at higher applied magnetic field. The complexes Mn(hfac) 2 (BNOci) and 
Mn(hfac) 2 (BNOBr) show ferrimagnetic behavior (Fig. 9). These compounds 
show narrow hysteresis with coercive force being less than 20 Oe (Inset of Fig. 9). 



3.2 

Conclusion 

Mn complexes with BNOci and BNObi have the same crystal structure as the 
complexes with BNOh and BNOp. From the observed intermolecular distance, 
the strongest interchain interaction in Mn(hfac) 2 (BNOR) appears to arise from 
the Mn — Mn and Mn — Cmiddie distances. From the dipole-dipole interaction 
term, the magnetic interaction between nearest neighbor magnetic moments is 
always antiferromagnetic. This argument is confirmed by the fact that 2J'/k^ 
does not exhibit a regular change throughout the Mn(hfac) 2 (BNOR) series, 
but changes abruptly in sign when the nearest interacting pair is altered. 
Consequently, the ferri/ferromagnet ground state is stabilized in complexes 
containing BNOci and BNObi (Fig. 10a), whereas complexes containing BNOh 
and BNOp behave as antiferromagnets (Fig. 10b). Due to a weak interchain 
exchange interaction, the antiferromagnetic complexes Mn(hfac) 2 (BNOn) and 
Mn(hfac) 2 (BNOp) display metamagnetic behavior characteristic of highly 
anisotropic magnetic materials. 



◄ 

Fig. 7. (a) View of a 1-D chain formed by a bisnitroxide BNOh radical and Mn(II)(hfac) 2 . 
(b) View of a chiral 1-D chain formed by a bisnitroxide BNOchirai radical and Mn(II)(hfac )2 
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Fig. 8. Temperature dependence of the low-field susceptibility (5 Oe) of the ferrimagnetic 
complexes (a) Mn(hfac)2 ■ (BNOh) and Mn(hfac)2 ■ (BNOp) and (b) Mn(hfac)2 ■ (BNOci) 
and Mn(hfac)2 ■ (BNOpr) 



The magnetic properties of chiral 1 -D polycrystals of [NITchirai ' 
Mn(II)(hfac)2]„ are similar to those of the ferrimagnetic complexes 
[BNOci • Mn(II)(hfac)2]„ and [BNObi ' Mn(II)(hfac)2]„. The magnetic prop- 
erties of polycrystals of chiral 1 -D magnets of [BNOchirai ' Mn(II)(hfac)2]„ are 
similar to those of the metamagnetic complexes [BNOh ' Mn(II)(hfac)2l„ and 
[BNOp • Mn(II)(hfac)2]„. 



4 

Two-Dimensional Metal-Nitroxide Systems 

In principle, there are no difficulties in the design of materials that undergo 
magnetic phase transitions at higher temperatures: the magnetic network of 
coupled metal ions and radicals should be extended from one to two or three 
dimensions, and the strong magnetic coupling between the spin centers would 
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Fig. 9. Magnetization curves of Mn(hfac )2 ■ (BNOBr) (•) and Mn(hfac )2 ■ (BNOh) (A) at 
1.7 K. Inset shows the low-field cycling 
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Fig. 10. Schematic drawing of the magnetic structure of (a) Mn(hfac )2 ■ (BNOci) and 
Mn(hfac )2 ■ (BNObi) and (b) Mn(hfac )2 ■ (BNOh) Mn(hfac )2 ■ (BNOp). Broken lines repre- 
sent the antiferromagnetic interaction between the ferro/ferrimagnetic 1-D chains. Solid 
lines represent the ferromagnetic interchain interaction 



ensure high transition temperatures [64]. One of the first two-dimensional 
structures is a series of copper(II) /f-diketonate complexes [65-71]. The 
coordination of the oxygen atoms of the nitroxide groups of neighboring 
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molecules completes the distorted octahedra around the Cu(II) ions. Tem- 
perature dependence of the magnetic susceptibility of these complexes fits a 
theory of isolated ferromagnetic coupled pairs. Enaminoketones A and B of 
3-imidazoline- 1-oxyl form layered polymeric structures [72-75] with Cu(II), Co, 
and Ni. However, while all these complexes exhibit extended structures, the 
magnetic coupling is limited to the directly-bonded metal-nitroxide couples. 

Thus, the problem is not only how to control the construction and the 
structure of extended systems in a desired fashion in order to extend the 
properties from the individual building blocks to the entire lattice, but also 
how to maintain strong magnetic coupling throughout the extended structure. 
Moreover, there is much difficulty in making a single crystal of the complex 
which has high dimensionality. Therefore, there are still only few examples of 
nitroxide-metal ion-based two-dimensional magnetic complexes (Table 2) [14, 
31, 32, 39, 52, 76, 77]. The flexibility of the ligands is important to make good 
crystals of two- or three-dimensional complexes. In metal-nitronyl nitroxide 
systems, there are limitations with regard to the topological synthesis of two- 
or three-dimensional complexes. The approach through high spin oligo- 
nitroxides and transition metal complexes is more advantageous than that 
through metal-nitronyl nitroxide systems. 



4.1 

Structure and Magnetic Properties of Ferrimagnetic 2-D Sheet Formed 
by Manganese(ll) and Nitroxide Radicals 



4.1.1 

2-D Manganese(ll) Complexes with Nitronyl Nitroxides 

Some ferrimagnetic 2-D manganese(II) -nitronyl nitroxide systems were 
synthesized by the complexation of manganese(II) bispentafluorobenzoate 



Table 2. Nitroxide-metal-based two-dimensional complexes 



Formula 


Crystal or 
powder 


//^B (/interlayer/^B) 


7c or Tn 


Ref. 


{Mn(pfbz)2l2 
(NITR), R = Me,Et 


Powder 




Tc = 24 K, 
Tc = 20.5 K 


[52] 


{Mn(pfpr)2l2 
(NITR), R = Me,Et 


Powder 


-149.6 K (0.075 K) 


Tc = 24 K, 
Tc = 20.5 K 


[76] 


{Mn(hafc)2l3 
(NITBzald )2 ■ 
O. 5 CHCI 3 


Powder 


-367 K 


Tc = 6.4 K 


[77] 


{Mn(hafc)2l3 

(TN0PB)2 ■ «-C7Hi5 


Crystal 


/MnNo/^B = -175.4 K, 
/nono/^b = 0.225 
( 0.002 ± 0.001 


Tc = 3.4 K 


[31, 32, 39, 78] 


{Mn(hafc)2l3(TNOP)2 


Crystal 


/MnNo/^B < -350 K 


Tc = 9.5 K 


[14, 79] 



Data were evaluated by quantum-classical approximation of ferromagnetic 3/2 spins, /miino and 
/nono define the interactions between Mn and NO and between two NO spins, respectively. 
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Scheme 7. 
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[Mn(pfbz)2, pfbz = pentafluorobenzoate] or manganese(II) bispentafluoro- 
propionate [Mn(pfpr)2, pfpr = pentafluoropropionate] with 2 -alkylnitronyl 
nitroxide (NITr, R = Me, Et) (Scheme 7 ). None of the complexes gave single 
crystals amenable for an X-ray crystal structure analysis. It is assumed that 
2 -D sheet structures were formed in these complexes. The complexes of 
[Mn(pfbz)2]2(NITR) (R = Me, Et) and [Mn(pfpr)2]2(NITR) (R = Me, Et) are 
ferrimagnetic and the transition into an ordered state occurs at ca. 23 K. The 
quantitative analysis of the magnetic susceptibility of these complexes was 
performed by using a model of a Heisenberg classical-quantum spin system 
(see Table 2 ). 



4.1.2 

Manganese(ll) Complexes with Triplet Bis-Nitroxide Radicals 

The X-ray crystal structure of the complex {Mn(hfac)2}3 • (TNOPB)2 • u-CyHis 
revealed an isostructure and that the manganese(II) ion has an octahe- 
dral coordination with the four oxygen atoms of two hfac anions bound to 
the metal ion in the equatorial plane while the axial positions are occupied 
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◄ 

Fig. 11. (a) View along the c axis of a layer displaying hexagons comprised of six 
trinitroxide TNOPB and six manganese complexes Mn(hfac)2. (b) View of the shortest 
contact between the layers 



by the two oxygen atoms of the two nitroxide groups. Six triradical molecules 
and six Mn ions form a hexagon from which an extended honeycomb network 
is constructed by sharing edges (Fig. 11a). A disordered n-heptane molecule is 
contained within each hexagonal cavity. The 2-D network sheets form a 
layered structure in which the adjacent layers can slide by a radius of the 
hexagon from the superimposable disposition with a mean inter-plane 
distance was 3.58 A (Fig. lib). On the basis of spin density and distance, 
the strongest inter-plane spin-spin interaction is considered to arise from the 
carbons (3.78 A apart) of the benzene rings para and meta to the nitroxide 
groups (Fig. lib). This type of interaction is expected to be ferromagnetic as 
dictated by McConnell’s theory [80, 81]. 

While satisfactory analytical data were obtained for the complex 
[Mn(hfac) 2}3 • (TNOB )2 [14], the complex did not afford a single crystal 
amenable to X-ray crystal structure analysis. It is assumed that, while a 2-D 
sheet structure is formed, the reduced symmetry of TNOB relative to TNOPB 
must be responsible for the difficulty in obtaining single crystals. Magnetic and 
analytical data support the 2-D nature of this complex. The 2-D network is 
schematically presented in Fig. 12. As in the case of TNOB, the network may 
be stacked ferromagnetically across layers. 

The magnetic susceptibility of [Mn(hfac) 2 } 3 (TNOPB )2 • n-CyHig and 
[Mn(hfac) 2 } 3 (TNOB )2 were investigated in the temperature range 1.8- 
350 K and fields up to 5 Tesla. For the complex [Mn(hfac) 2 } 3 (TNOPB )2 • 
n-CyHig with increasing temperature, this dependence is first characterized 
by a sharp maximum at 2.5 K followed by a smooth minimum in the range 
80-140 K, which rather resembles a plateau. The minimum value of 
Zmoi^ =5.71 emu K/mol within this interval is very close to the theoretical 
limit of 5.625 emu K/mol which is expected for three stable non-interacting 
S = 3/2 spins per mole. According to Fig. 11a, this spin configuration can 
be considered as being formed by one S = 5/2 spin of the Mn(II) ion 
coupled antiferromagnetically with two 1/2-spins of two different TNOPB 
triradicals. 

For the complex [Mn(hfac) 2 } 3 (TNOB) 2 , the Xmoi^’ values increased with 
decreasing temperature, displaying a maximum at 10 K. The value of 
Zmoi^’ =7.1 emu K/mol is larger than the theoretical limit of 5.625 emu 
K/mol, which is expected for three stable non-interacting S = 3/2 spins per 
mole. According to Fig. 12, this spin configuration can be considered as 
being formed by one S = 5/2 spin of the Mn(II) ion coupled antiferromag- 
netically with two 1/2-spins of two different TNOB triradicals. Together with 
the lack of a minimum at lower temperature, the room temperature 
value also points to the operation of strong (more negative than -300 K) 
antiferromagnetic coupling between the Mn(II) ion and the nitroxide radical 
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Fig. 12. A schematic drawing of the estimated 2-D network structure of 
{Mn(hfac)2}3 ■ (TNOB)2 



of TNOB, in which the onset of intramolecular ferromagnetic coupling is 
meaningful. 

When measurements were conducted in a low field, the magnetization 
values of {Mn(hfac) 2}3 • (TNOPB )2 • n-CyYii^ and {Mn(hfac) 2}3 • (TNOB )2 
complexes show a sharp rise at Tc = 3.4 K and 9.2 K, respectively (Fig. 13). 
The low-field susceptibilities below Tc are extremely large, as expected for a 
ferro/ferrimagnet. Spontaneous magnetization was observed below Tq, dem- 
onstrating the transition to a bulk magnet. Saturation magnetization values of 
these two 2-D complexes (Ms = 9 ^^/mo\.) agree closely with the theoretical 
limit assuming antiferromagnetic coupling between the manganese and triplet 
radical spins (Fig. 14). These compounds exhibit narrow hysteresis with 
coercive force of less than 20 Oe. 



4.2 

Conclusion 

It has been shown that perfect 2-D ferrimagnetic sheets exhibiting heterospin 
ferromagnetic (/i > 0)-antiferromagnetic Q 2 < 0 in Fig. 5b) networks together 
with ferromagnetic stacking of layers are realized in metal-nitroxide radical 
systems. 
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Fig. 13. Observed magnetization versus T plots for the complex {Mn(hfac) 2}3 ■ (TNOB) 2 . 
Measurements were obtained at a magnetic field of 1 Oe (O) and spontaneous magneti- 
zation (•) 




Fig. 14. Field dependence of the magnetization of {Mn(hfac) 2}3 ■ (TNOB )2 measured at 1.8 K 
(O) and 10.8 K (-t) 



While the kinetic instability of the triradical TNOB did not allow us to study 
the magnitude of the intramolecular ferromagnetic exchange coupling, it is 
expected to be considerably greater for TNOPB. In cases involving TNOB, the 
interaction forms by its symmetry with respect to an isosceles triangular 
geometry. The value for a TNOB moiety located between the two nitroxide 
radicals at positions 3 and 5 of the same benzene ring is estimated from that of 
m-phenylenebis(N-ferf-butylnitroxide) and analogues to be J/k^ - 200-500 K. 
The other interaction between the 4' and 3 (and 5) nitroxide radical centers is 
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estimated to be ca. 67 K. Altogether, the intramolecular ferromagnetic 
coupling in TNOB is estimated to be stronger in TNOPB and contributes to 
the higher Tq value of 9.2 K, in comparison to 3.4 K in the latter, where self- 
assembly with the assistance of Mn(II) ions is observed. 

5 

Three-Dimensional Metal-Nitroxide Systems 

In Sects. 3 and 4 we described the one- and two-dimensional complexes, 
respectively. For one-dimensional complexes the transition temperatures to 
ferri/ferromagnets are about 5-6 K, while for two-dimensional complexes, Tq 
rises to ca. 9 K. In order to synthesize complexes with higher transition 
temperatures, the dimensionality of the spin network must be raised. The 
flexibility of ligands is more important in the construction of suitable crystals 
in three-dimensional complexes than in two-dimensional systems. Owing 
to these difficulties, in metal-nitroxide systems only three-dimensional 
complexes have been reported. 



5.1 

Crystal and Molecular Structures 

An X-ray crystal structure analysis of an orthorhombic crystal of the 
complex revealed the formation of a parallel cross-shaped 3-D polymeric 
network (Fig. 15). The oxygen atoms of the terminal nitroxide groups of the 
triradical TNOP are ligated to two different manganese ions to form a 1-D 
chain in the b/c plane of the crystal. Since any manganese ion in an 
octahedral position is attached to the two nitroxide oxygens from two 
different triradical molecules in a trans disposition, the trinitroxide 
molecules are in a zigzag orientation along the chain. The diphenyl 
nitroxide unit is in a chiral propeller conformation and the R- and S-forms 
alternate along the chain. The middle nitroxide group of the ligand molecule 
TNOP on one chain is used to link its oxygen with that of the same chirality 
occurring in adjacent chains extended in the b/-c diagonal direction through 
a third Mn(II) ion in an octahedral position with an intersecting angle of 
54.4°, thus establishing a parallel cross-shaped 3-D polymeric network 
(Fig. 16). 



5.2 

Magnetic Properties of 3-D Systems 

The Xmo\T value of 8.64 emu K mol“^ at 300 K is larger than the theoretical 
value of 5.63 emu K mol~^ expected for short-range antiferromagnetic 
ordering of six 1/2 spins of TNOP and three 5/2 spins of Mn(II) for 
{Mn(hfac) 2 } 3 (TNOP) 2 . As T is lowered, the Xmoi^’ value increases monoton- 
ically in proportion to the increase in the correlation length within the 
network. Together with the lack of a minimum at lower temperature, the room 
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Fig. 15. Crystal structure of the 3 -D metal-radical complex {Mn(hfac)2}3 ■ (TNOP)2. The CF3 
and (CH3)C groups are omitted for clarity, a, b and c denote the orthorhombic crystal axes. 
The Mn(l) and Mn( 2 ) ions are represented by closed circles 



temperature Xmoi^’ value points to the operation of strong (more negative than 
-300 K) antiferromagnetic coupling between the Mn(II) ion and the nitroxide 
radical of TNOP, in which the onset of the intramolecular ferromagnetic 
coupling is meaningful. 

The temperature dependence of the magnetization M for a single crystal 
sample of [{Mn(hfac) 2 } 3 (TNOP) 2 ] was investigated at 5 Oe. When the sample 
was cooled within the field of 5 Oe, field-cooled magnetization displayed an 
abrupt rise at Tc = 46 K (Fig. 17). 

The field dependence of magnetization at 5 K showed two important 
features. First, the magnetization rose sharply at low field reached a value 
of ca. 9 /Tb (50,000 emu G moP^) at 220 Oe, and became saturated. The 
saturation value is in good agreement with the theoretical value of 9 ps 
(5/2 X 3 - 3/2 X 2 = 9/2) expected for the antiferromagnetic coupling be- 
tween the Mn(II) ion and the S - 3/2 triradical TNOP. Secondly, a 
conspicuous magnetocrystalline anisotropy was found in which the easy axis 
of magnetization lies along the c axis of the crystal lattice and the hard axis lies 
perpendicular to it (Fig. 18) [62, 63]. 




Fig. 16. Schematic drawing of a crystal of the 3 -D metal-radical complex 
{Mn(hfac)2}3 ■ (TNOP)2 




Temperature /K 



Fig. 17. Temperature dependence of magnetization of a single crystal of 
{Mn(hfac)2}3 ■ (TNOP)2. (a) Zero field cooled magnetization (ZFC, O) along the c-axis. 
(b) Zero field cooled magnetization (ZFC, A) along the fc-axis 

5.3 

Conclusion 

It is concluded that the perfect 3-D ferro/ferrimagnet with a heterospin 
ferromagnetic (Ji > 0)-antiferromagnetic {J 2 < 0) network is realized in 
{Mn(hfac) 2 } 3 (TNOP 2 ). The saturation value is indicative for the antiferromag- 
netic coupling between the Mn(II) ion and the S = 3/2 triradical TNOP. 
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Fig. 18. Magnetization curves of {Mn(hfac) 2}3 ■ (TNOP )2 at 1.8 and 25 K along the three 
principal crystallographic axes 



6 

Conclusions 

Supramolecular approaches have been successfully applied to the construction 
of high-dimensional network structures which are otherwise difficult to 
construct. These techniques utilize intermolecular interaction considerably 
weaker than the conventional covalent bond, e.g., van der Waals forces, 
hydrogen bonding, hydrophobic interaction, etc. The products are formed 
under thermodynamic, rather than kinetic control. Such network structures 
are intended for use as molecular switches, host molecular cages like “organic 
zeolites”, and enzymatic pockets, electric conductors, non-linear optical 
materials, etc. As far as the magnetic properties are concerned, the exchange 
couplings through these weaker bonds are rather weak. For the coupling to be 
strong, unpaired electrons have to be bound to each another through one or 
two (7-bonds or several 7r-bonds. Thus, the supramolecular approaches have 
been limited to metal coordination compounds. Design of appropriate radicals 
is key to our strategy of employing high-spin oligo-nitroxide radicals as 
bridging ligands. 
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A series of thiazyl radicals related to the trithiadiazolylium radical cation, S3NJ* are 
described. In many instances the compounds exist as spin-paired singlets which are 
consequently diamagnetic. However, when the dimerisation process can be inhibited, these 
open shell molecules exhibit very strong exchange interactions, characterised by Weiss 
constants, 6, up to 10^ K. Magneto-structural correlations show that, in the majority of 
instances, the magnetic exchange interactions are propagated via close intermolecular S- • N 
and S- • S contacts (typically in the region 3. 1-3.7 A). Examples of thiazyl radicals exhibiting 
magnetic ordering temperatures in excess of 50 K are described. In addition, the 
phenomenon of bistability (in which both open-shell monomeric and closed-shell dimeric 
forms are stable over the same temperature range) is discussed and examples of thiazyl 
radicals exhibiting bistability up to room temperature described. 

Keywords: Thiazyl radical. Organic magnet. Molecular magnet, Ferromagnet, Antiferro- 
magnet. Weak ferromagnet. Bistability 
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1 

Introduction 

In the last decade we have witnessed remarkable developments in the field of 
7i-electron magnetism, specifically in the field of organic molecular magnets 
[1]. Fifteen or twenty years ago, the thought that a purely organic compound 
containing no metal ions could exhibit the magnetic behaviour which we 
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typically associate with metals, their oxides, and ceramics was barely tenable. 
However, two reports [2, 3] in 1991 revolutionised the area of organic 
magnetism - suddenly the possibility that an organic material could exhibit 
long-range magnetic order became a reality. Work by Kinoshita and co- 
workers in Japan showed [2] that the open-shell radical p-nitrophenylnitronyl 
nitroxide, p-NPNN (Table 1) underwent a phase transition to a magnetically 
ordered state close to absolute zero (0.6 K), whilst Wudl’s group [3] found that 
the fullerene charge-transfer salt, Ceo ' TDAE (Table 1) ordered as a ferro- 
magnet at 16 K. Subsequently, there has been a continued flow of reports of 
organic radicals, especially nitroxide and nitronyl nitroxide radicals (Table 1) 
[4], that undergo magnetic ordering (ferromagnetism, antiferromagnetism, or 
weak ferromagnetism) at finite temperatures, although almost invariably 
below liquid helium temperature (4.2 K). In nearly all cases, the rather modest 
ordering temperatures can be ascribed to the weakness of the magnetic 
exchange interaction between neighbouring molecules. This arises because of 
the considerable steric protection required to inhibit dimerisation of the often 
reactive radical centre. In this article we review developments in thiazyi 
(sulfur-nitrogen) radical chemistry. In particular, we concentrate on a group of 
71 radicals that are sterically unhindered and which provide close approach of 
radical centres and the potential for strong magnetic exchange interactions. 



Table 1. Selected organic magnets and their transition temperatures (FM = ferromagnet, 
WFM = weak ferromagnet) 



Compound Ref. Compound Ref. 




[4] 




[4] 




[4] 
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2 

Trithiadiazolylium Cations and Related Radicals 

The trithiadiazolylium cation, SsN^*, has been known for well over 
100 years; it was first reported [5] by Demarcay in 1880, but it was not 
until the 1970s that its structures (both molecular and electronic) were 
elucidated [6]. A number of structures have been determined in the solid 
state [7] and all have been found to be dimeric with the two radicals 
associated via a long S- • S contact (ca. 2. 9-3.0 A) between rings. This S- • S 
contact is significantly longer than a conventional S — S bond (2.08 A) [8], 
although somewhat shorter than the sum of the van der Waals radii 
perpendicular to the heterocyclic ring (4.06 A) [9]. EPR studies [6] combined 
with theoretical calculations [6] showed that the unpaired electron resides in 
a n* orbital of fl 2 symmetry. The dimerisation process can therefore be 
considered as a four centre, two electron n*-n* bonding interaction between 
S atoms. A qualitative MO diagram of this interaction, together with the 
structure of [S 6 N 4 ]Cl 2 is shown in Fig. 1. 

A number of related radicals derived from the S 3 NJ* cation can be prepared 
by replacing by the isoelectronic groups R — C or N (see Scheme 1). These 
radicals will be described in more detail in later sections. More recently, there 
has been considerable success in developing selenium analogues of these 
radicals [10, 11]. In this article we will focus on the family of S/N radicals 1-5 
depicted in Scheme 1. 




Fig. 1. Qualitative MO diagram for the spin-paired association of two S3N+* radicals, and 
structure of [S3N2l2Cl2 
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Scheme 1. A family of heterocyclic radicals derived from S3N+* 



2.1 

Electronic Structure 

The electronic structures of these n radicals have been carefully probed by a 
combination of experimental (EPR) and theoretical calculations by a number 
of research groups and the results are summarised here. 



2.1.1 

Theoretical Calculations 

Theoretical calculations on SsN^* and its derivatives have been made at a 
variety of levels, from Hiickel theory [12] to ab initio methods [13] to Density 
Functional Methods [14]. The n manifolds of these thiazyi radicals are closely 
related to those of the cyclopentadienyl anion, Cp~ [12a]. However, the 
degeneracies of the n system observed in Cp~ are lifted due to the lower 
symmetry and mismatch in energies of the constituent atoms (C/N/S). In all 
cases the symmetry of the singly occupied molecular orbitals correspond 
closely to either 714 or of Cp~. The n manifold of Cp~ is shown in Fig. 2, 
along with the SOMOs of 1-4. Both 1 and 2 have SOMOs of 02 symmetry under 
the C 2 V point group corresponding to 714 of Cp~. Radicals 3 and 4 have orbitals 
corresponding to of Cp“. Notably in derivatives of 1-4 there is little 
delocalisation of unpaired spin density away from the heterocyclic ring. This 
arises because of the difference in the energies of the 71 -manifolds of the 
heterocyclic ring and the aryl substituent. Calculations on the open shell 
radicals also indicate significant polarity of the S — N bond, in the sense 
The localisation of the spin density and the polarity of the 
heterocyclic ring have important consequences for the magnetic behaviour of 
thiazyi radicals (see Sect. 2.2.4). 
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7T molecular orbitals of 
the cyclopentadienyl 
anion, Cp’ 
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Fig. 2. The 7c-orbitals of Cp and the singly occupied molecular orbitals of 1-4, corre- 
sponding to either 714 or 715 of the Cp^ ring 
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2.1.2 

EPR Spectra 

EPR studies on 1-5 all show extensive delocalisation around the heterocyclic 
ring. Their EPR spectra are characterised by isotropic g-values typically in the 
region 2.005-2.010 with the deviation from the free electron value caused in 
large part by the larger spin-orbit coupling associated with S. Indeed, where 
selenium analogues are known, this gives rise to even greater shifts in the g- 
factor [10b, 15]. In addition, frozen glass spectra indicate considerable 
anisotropy in the g-value. These radicals exhibit a rhombic EPR spectrum with 
one set of large hyperfine coupling constants to N in the high-field component. 
This is entirely consistent with a 71 -based radical in which the 71 -type nature of 
the orbital means that significant hyperfine coupling is seen in only the 
direction perpendicular to the heterocyclic ring plane. This is clearly 
evidenced in the EPR spectrum of PhCNSSN (Pig. 3). Some typical EPR data 
for derivatives of 1-5 are given in Table 2. There is typically a good 
correspondence between the calculated unpaired spin density distributions 
(both s and p contributions) from EPR [10b, 15] and those predicted from 
theoretical calculations [13, 14]. 
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Fig. 3. Frozen glass X-band EPR spectrum of PhCNSSN at 77 K. Note the rhombic pattern 
and large hyperfine coupling to N (marked *) to just the gi component, characteristic of 
these 7c-based radicals 



2.2 

Association of Thiazyi Radicals 



2.2.1 

Structural Studies 

A number of structural studies on 1 have shown a strong tendency for these 
7i-radicals to associate through an out-of-plane n*-n* interaction (described 
previously and shown in Fig. 1). In derivatives of 1, the S- • S separation 
between monomers is typically 2. 9-3.0 A [7], somewhat less than the sum of 
the van der Waals radii perpendicular to the ring. 

Whilst only one mode of association (trans-antarafacial) has been 
observed for 1 , the modes of association of the dithiadiazolyl radical 2 are 
considerably more diverse. The mode of association has been clearly linked 
to the bulkiness of the R substituent. When R is not sterically demanding 
then the common mode [16] of association is a ds-cofacial geometry (Fig. 4). 
For more sterically demanding groups (e.g., R = Me, CF 3 , Me 2 N, adamantyl) 
[13a, 13c, 17] then a twisted geometry is obtained (Fig. 4b). More recently, 
phenyl substituents containing a strong dipole across the molecule have also 
been found to induce twisting to maximise dipolar interactions [18]. In one 
reported instance a frans-antarafacial geometry has been observed [16b] 
(Fig. 4c). A fourth trans-cofacial mode of association has also been reported 
(Fig. 4d) [19]. With the exception of this latter mode, all these associative 
conformations exhibit one or more intradimer S- • S contacts, typically in the 
range 2.9-3. 1 A [16-18]. In the case of the trans-cofacial mode depicted in 
Fig. 4d, the S- • N contacts are in the range 3.12-3.35 A, similar to the S- • N 
contacts in the twisted configuration (3.18-3.24 A) [13a]. In a few instances 
monomeric or intermediate structures have been observed, and these will be 
described later. 

In all these modes of association, a bonding interaction between SOMOs can 
be clearly identified. Theoretical calculations by Mews and co-workers [13a] 
indicated that the calculated difference in energy between the modes of 




Table 2. EPR data for derivatives of 1-5 [hyperfine coupling constants are quoted in G] 
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Harrison S, Pilkington, RS, Sutcliffe LH (1984) / Chem Soc Faraday Trans I 80: 669. 
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Fig. 4. Modes of association of dithiadiazolyl radicals: (a) cis-cofacial, (b) twisted, 
(c) frans-antarafacial, (d) frans-cofacial 



association was small ('^5 kj mol”^) and the overall geometry is likely to be 
dictated, in a large part, by other packing forces. 

Only one structure of a 1,3,2,4-dithiadiazolyl radical, 3, has been reported; 
the diradical P-SNSNCC 6 H 4 CNSNS also associates in the solid state in a trans- 
cofacial fashion, held together by S- • S (3.21 A) and S- • N (3.35 A) contacts 
[20]. Attempts to prepare other derivatives have been thwarted by the 
propensity of derivatives of 3 to rearrange to 2 either in solution [21] or in the 
solid state [22]. More recently, Passmore prepared crystals of PhCNSNS, 
although these crystals decomposed in the X-ray beam [23]. 

A number of structures of the 1,3,2-dithiazolyl ring, 4, are known and two 
common modes of association have been observed, either cfs-cofacial [24] or 
trans-cofacial [25] (Fig. 5). In both cases the dominant interaction would 
appear to be intradimer S- • S contacts (3.14-3.18 A). However, for many 
derivatives of 4, a far greater diversity of structure is observed, with many 
compounds exhibiting stacked structures with intermolecular separations 
significantly greater than 3.5 A. These will be discussed later. 

Only one derivative of 5 has been structurally characterised (5, = CeFs, 

= Cl) [26]. It too associates, in a twisted configuration (to avoid steric 
hindrance between perfluoroaryl groups) with an intradimer S- • S separation 
of 3.30 A. 




Fig. 5. Crystal structures of (a) [F3CCSNSCCF3]* and (b) C6H4S2N' 
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2.2.2 

EPR Investigations of the Association Process 

In solution, these dimers typically dissociate into monomers but a monomer- 
dimer equilibrium exists (see Eq. 1). The equilibria can be monitored by UV- 
visible spectroscopy [27] or by measuring the intensity of the EPR spectrum 
[6c, 21b, c, 28]. Careful EPR studies have allowed the thermodynamics of the 
equilibrium to be determined in some cases. Some data are given in Table 3. 
Its importance is that it provides an estimate of the dimerisation energy; if 
these thiazyl radicals are to be used in the design of magnetic materials, then 
this dimerisation energy must be overcome. As can be seen from Table 3, the 
dimerisation energy for 1 is ca. 45 kj mol“^ and invariably the solid state 
structures of derivatives of 1 are dimeric. In comparison, replacement of one S 
by C, leads to a reduction in dimerisation energy such that AHjim for 2 and 3 
are ca. 35 and 20 kJ moP^ respectively. Whilst the majority of derivatives of 2 
are also dimeric in the solid state, a few are monomeric and will be discussed 
in more detail later. For 4, replacement of N by C leads to a further reduction 
in AEf^im ( 0 kJ mol and a large number of these radicals are found to be 
monomeric in the solid state. 



2.2.3 

Theoretical Investigations into the Dimerisation Process 

A number of theoretical studies have been made on the dimerisation of 1, 2, 
3, and 4. All conclude that there is a n*-n* interaction between singly 
occupied molecular orbitals (SOMOs) which leads to a net bonding 
interaction and a closed shell singlet ground state [13a, 25, 29]. A variety 
of modes of association is possible that give rise to efficient overlap of n 
orbitals, and all require the rings to be essentially eclipsed, or rotated by 90° 



Table 3. 



12 3 4 




“ Ref. [6c]; 

’’ Ref [21bj; 
‘ Ref [21c]; 
Ref [28], 
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or 180° with respect to one another and these are all observed structurally. A 
more quantitative analysis of the association has proved problematic, and we 
highlight this by a recent, detailed ab initio study [13d] of the dimerisation 
process by Oakley and co-workers at the University of Guelph. They 
examined how the theoretical model of the dimerisation process changed as a 
function of the basis set employed, paying particular attention to the refined 
geometry and energy of dimerisation. Initial findings showed that valence 
minimal basis sets produced significant distortions of the dimeric structures 
and appeared to maximise N- • N rather than S- • S inter-ring bonding 
interactions. Additional d-functions (3-210’*^) were successfully employed and 
produced much better agreement with the observed crystallographic para- 
meters. However, the calculated energies of dimerisation proved to be 
particularly sensitive to the theoretical method used. At low orders of 
perturbation theory (MP2) the dimerisation energy was found to be 
176 kj moU^ but at third (MP3) and fourth (MP4SDQ) order these energies 
reduced successively to 27 and then 17 kJ moU^ (cf. experimental values of 
ca. 35 kJ mol“^ in Table 3). This variation in enthalpy of dimerisation with 
basis set provides a substantial theoretical problem and calculated enthalpies 
of association should not be taken quantitatively. It should, of course, be 
remembered that these calculations are gas-phase and the energies considered 
are small. In such circumstances, both solvent and certainly lattice effects 
may provide a significant perturbation. Nevertheless, a set of comprehensive 
studies using the same basis set on derivatives of 1-5 might, however, provide 
a useful correlation (e.g., between calculated and observed values of AHdim) 
for future work. 



2.2.4 

Implications of the Dimerisation Process on the Design of Thiazyi-Based Magnets 

These thiazyi radicals exhibit a tendency to associate both in solution and in 
the solid state through a spin-paired dimerisation process, generating a closed 
shell singlet ground state. If they are to be utilised in the design of molecular 
magnets, then this dimerisation process must be overcome. The smaller the 
dimerisation enthalpy is, the more favourable is the possibility of forming 
paramagnetic structures in the solid state. This is borne out by a simple 
correlation of the structural data available on the known derivatives of 1-5 and 
the estimates of the dimerisation energy for derivatives of 1-5. Evidently, 
severe steric and/or electronic demands must be placed on radicals such as 1-3 
if the dimerisation process is to be inhibited. In contrast, derivatives of 4 are 
likely to be monomeric in many instances. 

There are several further important factors which should be considered 
here. Firstly, the unpaired spin density in these thiazyi radicals is, to a large 
degree, localised on the heterocyclic ring. The consequence is that we can 
vary the substituents to modify the molecular structure, without signifi- 
cantly changing the electronic structure. Secondly, the thiazyi bonds are 
extremely polar and there is a great tendency for these radicals to pack in 
such a way as to maximise the intermolecular interactions. This 
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leads to close approach of regions of unpaired spin density. Provided that 
the out-of-plane SOMO-SOMO dimerisation process can be inhibited, the 
close approach of radical centres should lead to strong exchange interac- 
tions. Since the heterocyclic N and S atoms both bear considerable 
unpaired spin density the S- • N interactions are likely to propagate 
antiferromagnetic exchange interactions (on the basis of McConnell theory 
[30]). In the sections which follow we examine the magnetic properties of 
derivatives of 2 and 4. 



3 

Magnetic Properties of 1,2,3,5-Dithiadiazolyl Radicals 

In this section we examine how the dimerisation process (ca. 35 kj moP^) can 
be overcome through a combination of electronic and steric factors. Since the 
spin density can be assumed to be localised on the heterocyclic ring, then 
particular attention will be paid to the intermolecular interactions between 
dithiadiazolyl rings since these are likely to propagate the magnetic exchange 
pathway through the solid. Correlations between magnetic behaviour and solid 
state structure will be made. 



(a) 



(b) 




Fig. 6. Crystal structures of (a) 2a and (b) 2d 
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3.1 

Monomer or Dimer? 

The majority of dithiadiazolyl radicals associate as spin-paired dimers in the 
solid state in which the intermolecular separation is in the region 2.9-3. 1 A. 
However, one example has been reported in which the intermolecular 
separation is larger. The bis-radical 2a also has a long intermolecular stacking 
distance with the shortest contact along the stacking direction (see Fig. 6a) at 
3.48 A [31]. The long intermolecular separation is still within the sum of the 
van der Waals radii of S (4.0 A). Radical 2a exhibits some paramagnetism at 
room temperature (equivalent to ca. 30% unpaired spins), but its susceptibility 
drops rapidly such that it is rendered diamagnetic by 200 K. A second 
radical 2b was reported to exhibit a stacking motif with a regular lattice 
spacing of 3.54 A, coincident with the crystallographic c-axis. [18a] However, 
2b appeared to be diamagnetic throughout the temperature range studied, 
despite the longer intermolecular separation! A recent re-investigation [18b] of 
the crystal structure of 2b showed that a super-cell was present in which there 
is regular packing of dimers along the 7i-stacking direction. As with other 
dimeric dithiadiazolyl radicals, the short intra-dimer contacts render the 
material diamagnetic. 




2c 



3.2 

Radical Cations 



F 





A number of derivatives of 2 have been reported in which the radical bears a 
cationic substituent. Under these circumstances the ionic packing forces may 
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dominate the packing and overcome the dimerisation energy. In 1993, 
Banister and co-workers reported [32] the radical cation 2c. Whilst there was 
no crystallographic data, the magnetic moment at room temperature was 
1.8 /Tb and it exhibited Curie-Weiss behaviour with 6 = -65 K. Such large 
values for the Weiss constant will appear as a common theme in this area of 
chemistry and are consistent with strong communication between radical 
centres. More recently, the radical cation 2d was isolated [33]. Its structure 
is shown in Fig. 6. This compound too has a room temperature moment 
consistent with one unpaired electron (ca. 1.3 Pb)- On cooling the 
susceptibility passes through a broad maximum around 275 K consistent 
with the onset of short-range antiferromagnetic order. It then decreases 
slowly to ca. 50 K before undergoing an abrupt transition to a diamagnetic 
state at 15 K. This has been ascribed to a structural phase transition to a 
diamagnetic (dimeric) state [33]. 



3.3 

Fluorinated Radicals 

In recent years our own research has developed a series of perfluoroaromatic 
derivatives of 2. In these derivatives the electrostatic repulsions between 
the orfho-fluorine and the heterocyclic ring N leads to a large twist angle 
between the two rings (30-60°), thereby making the substituent more 
sterically demanding. In addition, fine control of the electronic effects of the 
substituent have been achieved. For a series of derivatives 2e-k, a correlation 
is observed (Fig. 7) between both the electron-withdrawing/electron-releasing 
nature of the para-substituent (Hammett parameter) and the ease of 
reduction, i.e., radical stability. In addition the structures which appear to 
bear strongly electron-withdrawing groups have been found to be mono- 
meric. 
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Fig. 7. Correlation between the Hammett parameter, Cp of the para-substituent and the 
reduction potential of the corresponding [p-XC6F4CNSSN]^ ion. Solid state structures have 
been determined for most derivatives P-XC6F4CNSSN. Those marked (® ) are dimeric, those 
marked (O) are monomeric 



3.3.1 

Magnetic Properties of p-BrC6p4CNS5N 

The electronic structure of p-BrC 6 F 4 CNSSN, 21, has been investigated [34] by 
EPR studies (X-band solution EPR and ENDOR) and DET calculations. 
Solution EPR spectra clearly resolve the nitrogen hyperfine interaction, and the 
much smaller F and Br couplings can be detected in the ENDOR experiments. 
An analysis of the spin density distribution is in excellent agreement with DET 
calculations. The resultant spin density distribution is shown in Eig. 8 and 
clearly shows that essentially all the unpaired spin density is localised on the 
N 2 S 2 fragment of the heterocyclic ring. The dominant magnetic exchange is 
therefore likely to be via close intermolecular contacts between heterocyclic 
rings. 

The structure of 21 is shown in Eig. 9. In 21 the intermolecular contacts 
along the stacking direction (crystallographic a axis) are in the range 3.67- 
4.00 A, although the twist between the rings precludes efficient orbital overlap 
[35]. Instead molecules link together via electrostatic • •N‘^~ interactions 
(S- • N = 3.18 A) along the crystallographic c axis. These chains are isolated 
from each other by ca. 10.2 A, so that the magnetic exchange pathway can, at 
best, be considered to be two-dimensional. The room temperature magnetic 
moment (1.5 is close to that expected for an S = 1/2 ion. Above 60 K the 
compound exhibits Curie-Weiss behaviour (0 = -27 ± 1 K). Using a mean 
field approximation, the compound would be expected to undergo a 
paramagnetic-antiferromagnetic phase transition at 27 K. However, no evi- 
dence of long-range order is observed down to 1.8 K, consistent with the low- 
dimensional nature of the structure. Attempts to model the behaviour using 
the Bonner-Eischer model for a Heisenberg chain of S = 1/2 ions with a 
constant exchange term (/) proved unsuccessful since the observed data did 
not exhibit a maximum in the susceptibility. An excellent fit was obtained 
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Fig. 8. (a) Solution EPR spectrum of 21 and (b) theoretical spin density distributions. 
Estimated values from experimental measurements are given in parentheses 

(Fig. 9b), however, if / was allowed to be temperature-dependent, such that the 
magnetic exchange between radical centres decreased on cooling. Further 
experiments will be needed to determine whether there is any physical basis 
for this phenomenon, e.g., changes in intermolecular contact, etc. 



3.3.2 

Magnetic Properties of P-NCC6F4CNSSN 



3.3.2.1 

Polymorphism 

A recurrent theme in thiazyl chemistry and indeed in the area of organic 
molecular magnetism is the propensity for compounds to crystallise in two 
or more space groups, i.e., polymorphism. Since the solid state structure 
determines the strength and nature of the magnetic exchange interaction, then 
the magnetic behaviour of different polymorphs cannot be expected to be the 
same. This phenomenological problem has been reviewed recently [36]. 

In the area of dithiadiazolyl and diselenadiazolyl chemistry, at least four 
compounds have been reported to be polymorphic [16b, 37-39] in recent 
years. The dithiadiazolyl radical P-NCC 6 F 4 CNSSN, 2m, provides a classic 
example of this phenomenon. Radical 2m crystallises in either of two space 
groups [38, 40], depending on the sublimation conditions: sublimation under 
high vacuum at low temperatures leads to the formation of the kinetic 
phase, a-2m (triclinic, P-1) [38] whereas slow sublimation at more elevated 
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Fig. 9. (a) View of 21, viewed down the crystallographic a axis, illustrating the close in-plane 
S- • N interactions (3.18 A) and (b) Plot of molar susceptibility of 21 as a function of 
temperature. The fit to a Bonner-Fischer model with temperature-dependent / 
[/ = -(0.31T -I- 8.2)] is shown as a solid line in the Figure. The dashed line represents the 
best curve-fit to the high temperature data with a fixed / term 



temperatures leads to the thermodynamic phase, p-2m (orthorhombic, Fdd2) 
[40]. The structures of the a- and ;S-phases are shown in Fig. 10. There are 
some marked similarities in the two phases: both phases are monomeric with 
molecules linked together into chains via electrostatic CN‘^“- • interactions. 
Studies on thin films of 2m indicate [41] that this CN- • S interaction plays a 
role in directing crystal growth. Indeed in P-2m, the crystals grow as long 
needles with the CN- • S interaction (along the crystallographic c axis) 
coincident with the needle axis of the crystal [42]. However in the (x-phase, 
these chains of radicals crystallise in an antiparallel fashion with inter-chain 
interactions between the S and the cyano-N. In the ;S-phase, the interchain 
interactions are between the S and the heterocyclic N and generate coparallel 
alignment of the chains. The lack of an inversion centre in ^-2m has important 
ramifications for its magnetic behaviour. 



3.3.2.2 

Electronic Structure 

The electronic structure of 2m has been probed by solution EPR studies 
(X-band and ENDOR) and also through DFT calculations [34]. The electronic 
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Fig. 10. Crystal structure of (a) a-p-NCC6F4CNSSN and (b) /!-p-NCC6F4CNSSN 



Structure of 2m is almost identical to that of 21 and there is good agreement 
between theoretical and experimental spin density distributions. 



3.3.2.3 

Magnetic Behaviour of (i-p-NCC^ 4 CNSSN 

The a-phase of 2m exhibits a room temperature susceptibility of ca. 1.6 Pb 
consistent with an S = 1/2 ion [38]. It obeys Curie-Weiss behaviour above 
175 K with 6 = -25 K. On cooling further, deviation from Curie-Weiss 
paramagnetism is observed, with the sample passing through a minimum in 
at about 8 K; a little lower than that expected from the mean field theory 
approximation. The problematic synthesis of the (x-phase (the thermodynam- 
ically stable jS-phase grows preferentially) has precluded further study of its 
low temperature behaviour. These studies should confirm whether the 
compound undergoes long-range order at 8 K or whether it is a short range 
phenomenon. An analysis of the magnetic exchange pathways [42] indicates a 
number of intermolecular S- • S contacts (3.60-3.67 A) which should propa- 
gate the magnetic exchange in at least two dimensions. 



3.3.2A 

Magnetic Behaviour of P-P-NCC 6 F 4 CN 5 SN 

The jS-phase of 2m also exhibits a room temperature susceptibility of 1.6 It 
also obeys Curie-Weiss behaviour above 120 K, but with a much larger Weiss 
constant, 9 - -102 K indicative of a strongly antiferromagnetically coupled 
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regime [40]. As the temperature decrease further, x passes through a broad 
maximum at 60 K, indicative of the onset of short-range order. A fit of the high 
temperature data to a Bonner-Fischer expression for an S = 1/2 Heisenberg 
linear chain yields an intrachain exchange term of -83 K, but also a large 
interchain term [z/in,er//intra 0.1]. The large magnitude of the interchain 
interactions means that P-2m is not a good one-dimensional system [40]. An 
analysis of the intermolecular contacts shows that each molecule forms four 
close heterocyclic-N- • S contacts to other molecules in a distorted tetrahedral 
geometry. Propagation of these magnetic exchange interactions through the 
lattice lead to a distorted diamond-lattice for the magnetic exchange 
interactions [43]. Inclusion of the related radical P-NCC 6 H 4 CNSSN, which 
has a very similar spin density distribution, to 2m into the host-lattice of PHTP 
has confirmed that there is negligible exchange via the CN- • S interaction [43]. 

On cooling below 36 K, the susceptibility of p-2m becomes field dependent 
and an out-of- phase component in the ac susceptibility is observed (Fig. 11), 
consistent with the onset of long range order [44]. The field-dependent 
magnetic behaviour of P-2m below 36 K can be fitted to Eq. 1, 

M = yH-FMs (1) 

The small residual spontaneous magnetisation observed in zero-field can be 
explained in terms of spin-canting of the antiferromagnetic lattice leading to a 
small spontaneous magnetic moment. Hysteresis can be observed (Fig. 11) at 
low temperature with a small coercive field and remnant magnetisation. 

Canted antiferromagnetism (weak ferromagnetism) is favoured in low- 
dimensional systems and results from antiferromagnetic alignment of spins on 
two magnetic sublattices which are equivalent in number but not exactly 
antiparallel [45]. The sublattices are canted in the absence of a magnetic field 
provided that the total symmetry is unchanged and that this leads to a net 
magnetisation. Spin-canting cannot be observed in compounds whose spins 




Fig. 11. (a) ac susceptibility of /?-p-NCC 6 F 4 CNSSN below 80 K and (b) hysteresis of 
JS-P-NCC 6 F 4 CNSSN recorded at 1.8 K 
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(a) Pressure (kbar) 



(b) Pressure (kbar) 



Fig. 12. (a) The variation of Tc of /?-p-NCC6F4CNSSN with applied pressure, and (b) the 
variation of the spontaneous magnetisation, Ms, with applied pressure 



(on different sublattices) are related through an inversion centre (e.g., a-2m) 
but can be achieved in non-centrosymmetric space groups, such as ^-2m. 
Canting arises through the Dzyaloshinskii-Moriya interaction which 
minimises the exchange energy when the spins are orthogonal. The magnitude 
of the canting angle is proportional to the anisotropy in the g-tensor (Eq. 2) 
and, as a consequence, can be estimated from EPR studies. The canting 
angles determined from the saturation magnetisation extrapolated to 
0 K (1.5 X 10 ^ Pb)> and predicted from EPR studies are in good agreement, 
cf 0.085° and 0.117° [46]. 

6 = 1/2 arctan(|d//|) where \d/J\ ~ (g - gj/g (2) 

The onset of long-range magnetic order has been probed using specific heat 
capacity studies [44], powder [44] and single crystal [47] neutron diffraction, 
and muon-spin relaxation [48]. Single crystal symmetry precludes c-axis spin- 
alignment for a canted structure and the resultant magnetic moment must 
align parallel with the crystallographic a or b axes [44]. A careful analysis of 
the single crystal, powder, and muon data is in agreement that the moments 
are oriented parallel to b, with a small canted moment along the crystallo- 
graphic a axis. 

The magnetic behaviour of P-2m is truly exceptional. The strength of the 
magnetic exchange interaction, propagated through the close S- • N contacts, 
gives rise to a Tc of 36 K - the first example of an open shell radical to exhibit 
a Tc above liquid helium temperature. Recent studies have shown [49] that the 
application of pressure leads to a further increase in Tq. This presumably 
arises out of a shortening of the intermolecular S- • N contacts thereby 
enhancing the magnetic exchange term, /. The increase in Tq appears linear 
with the applied pressure and fits the empirical formula, 

Tc(P) = Tc(Po) + flPTc(Po) (3) 

It should be noted that the ordering temperature under a pressure P is 
dependent not only on the parameter a but also on Tc(P = 0). Since Tc(P = 0) 
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is large, then the enhancement in Tc is marked, reaching 50 K at 7 kbar 
(Fig. 12). Notably, however, the increase in the magnitude of the exchange 
term / also leads to a reduction in \d/J\. Necessarily (according to Eq. 2), this 
also leads to a reduction in canting angle and the value of the spontaneous 
magnetisation also decreases (Fig. 12) [49]. 

4 

Magnetic Properties of 1,3,2-Dithiazolyl Radicals 

In comparison to the dithiadiazolyl radicals, 2, described above, the dithiazolyl 
radicals, 4, exhibit dimerisation energies close to 0 kj mol~\ As a conse- 
quence, many are found to be monomeric in the solid state. We will separate 
the magnetic behaviour of this class of radical into two groups based on some 
structural motifs: those which contain 7i-stacked structures in the solid state 
and those which do not. We will concentrate in the first instance on those in 
which TT-stacking is not observed. 




4a 



4b 



4.1 

1,3,2-Dithiazolyl Radicals without ;i-Stacking 



4.1.1 

Methylbenzo- 1,3,2-Dithiazolyl, MBDTA 

The MBDTA radical was first reported by Wolmershauser et al. in 1984 [50]. 
They reported that both the benzo-derivative BDTA, 4a and the tolyl 
derivative, MBDTA, 4b, were paramagnetic at room temperature. Subsequent 
studies by Passmore showed that 4a was in fact a dimer [25]. However, our 
own structural studies on MBDTA indicate that it is indeed monomeric [51]. 



4.1. 1.1 

Electronic Structure 

EPR studies on 4b exhibit a triplet pattern (g = 2.003, a^ = 11.4G) consistent 
with localisation of the unpaired spin density on the heterocyclic ring [51]. The 
electronic effect of the methyl group is likely to be minimal and the spin 
density distribution can be expected to closely parallel that of 4a calculated by 
Passmore and co-workers with over 75% of the unpaired spin density residing 
on the C 2 S 2 N ring. Whilst there is evidently more 71 -delocalisation of the 
unpaired spin density in derivatives of 4, (cf. 2) it is by no means extensive. 
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Thus, heterocyclic contacts are likely to continue to play a dominant role in 
determining the magnetic properties of these radicals. (The SOMO of 4b is 
shown in Fig. 13). 



4.1. 1.2 
Structure 

Radical 4b crystallises in the orthorhombic space group Pbca. The asymmetric 
unit of 4b contains a single molecule of unexceptional geometry, with the fused 
ring essentially planar. Each molecule has six nearest neighbours (Fig. 13) with 
intermolecular contacts in the range 3.71-3.82 A. This set of contacts lie in the 
crystallographic hc-plane generating a two-dimensional network. Propagation 
along the crystallographic a direction relies on a close N- • H interaction at 
3.7 A to an aromatic-H atom, and the degree of spin delocalisation onto the 
ring [51]. 



4.1. 1.3 

Magnetic Behaviour 

The room temperature moment of 4b (like the dithiadiazolyl derivatives, 2, 
described above) is 1.3 pe close to, but a little less than that expected for an 




Fig. 13. Structure of 4b with (inset) the symmetry of the SOMO of the closely related 
radical, 4a 
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S = 1/2 paramagnet. This is entirely consistent with the presence of strong 
antiferromagnetic exchange interactions between radical centres, a common 
theme of thiazyi radical chemistry. The susceptibility passes through a broad 
maximum at 140 K indicative of the onset of low-dimensional antiferro- 
magnetic order. The exceptionally high temperature for y^ax is indicative of 
the strength of the antiferromagnetic exchange interactions between radicals. 
Below 25 K the susceptibility curve increases rather sharply most likely due 
to the contribution from a small mole fraction of non correlated 
paramagnetic centres in the sample (<1% radical defects). The low 
dimensionality of the high temperature magnetic behaviour is consistent 
with the two-dimensional nature of the thiazyi interactions observed in the 
solid state structure. An analysis of the high temperature susceptibility as a 
two-dimensional square-planar Heisenberg model of S = 1/2 centres pro- 
vided a good agreement between calculated and experimental data with 
/ = -72 K. The adequacy of the fit is reflected in both the position and 
height of the maximum in y. In Fig. 14 the magnetic susceptibility of 4b, 
corrected for diamagnetic and uncorrelated paramagnetic contributions is 
represented in reduced units and compared with representative examples of 
the S = 1/2 square planar Heisenberg model as taken from the literature [51]. 



0.050 




0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 

kj / |J| S(S+1) 

Fig. 14. Experimental temperature dependence of the susceptibility, in reduced units, of 4b 
(•) compared to representative examples of the S = 1/2 square planar Heisenberg model; 
[Cu(C5H5NO)6(BF4)2] (A), CuF2 ■ 2H2O (□), (2-NH2-5-Cl-C5H5N)2CuBr4 (O), (2-NH2-5- 
CH3-C5H5N)2CuBr4 (■), [Cu(pz)2-(N03)][PF6] (♦), [Cu(pz)2][C104] (T). Data are scaled to 
the high temperature series prediction (solid line) 
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4.1.2 

Naphthalene- 1,3,2-Dithiazolyl, NOT A 



4.1. 2.1 

Crystal Structure 

The crystal structure of NDTA, 4c, [triclinic, P-1] was reported [52] by Oakley 
and contains two molecules in the asymmetric unit. The molecules are packed 
in a close-packed herringbone arrangement similar to that observed for 4b and 
the closest intermolecular contacts are in the range 3.60-3.87 A, a little shorter 
than those observed in 4b. 



4.1. 2.2 

Electronic and Magnetic Behaviour 

The EPR spectra of 4c (g = 2.0067, aN = 11.06G) indicates significant 
localisation of the unpaired spin density on the heterocyclic ring [53]. An 
analysis of the susceptibility of 4c indicates that at room temperature it is 
paramagnetic with a magnetic moment of ca. 1.5 analogous to 4b. On 
cooling the susceptibility passes [52] through a broad maximum around 200 K 
before undergoing a phase transition to a more strongly antiferromagnetically 
coupled state. It passes through another maximum in its susceptibility around 
25 K with eventual loss of paramagnetism [53]. 



4.2 

1,3,2-Dithiazolyl Radicals with 7i-Stacking 

Considerable attention has been paid in recent years to the development of 
neutral radicals, such as 2 and 4, in the design of conducting materials [54]. 
Overlap of the n* SOMOs of these neutral radicals along the stacking direction 
has been proposed to lead to a conduction band [55]. The extent of overlap 
evidently depends on the closeness of the intermolecular separation along the 
stacking direction. If the overlap is small then the unpaired spins can be 
considered localised, whereas if the overlap is good a band forms and the 
magnetic properties more closely match those of a low-dimensional conductor 
[56]. In many instances a charge-density wave driven (CDW), Peierls 
distortion of the lattice occurs and the stack becomes irregular with alternating 
long and short contacts. This leads to a band gap and the material becomes 
insulating. From a magnetic viewpoint, this can be envisaged as a simple spin- 
paired dimerisation process and the resultant species is diamagnetic. This 
Peierls distortion is a phenomenological problem for many low-dimensional 
systems and a number of methods have been employed to inhibit this process 
(metal-insulator transition) [54]. Doping of the structure with halogens, for 
example, removes some electrons from the conduction band and this inhibits 
dimerisation, although more complex CDW-driven periodic distortions of the 
lattice can also occur to open up a band gap at the Fermi level [57]. 
Alternatively, maximising inter-stack interactions [39] can also prove to be 
beneficial since these provide an energy barrier to inhibit the distortion. 
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As a consequence, the magnetic properties of 7i-stacked radicals can often 
be complex, with magnetic moments significantly below that anticipated for 
S = 1/2 spins, the exact value depending on the extent of the development of 
the band structure. 





4.2.1 

Quinoxaline- 1,3,2-Dithiazolyl, QDTA 

This radical, 4d, was first reported by Wolmershauser in 1990 [58], a detailed 
EPR investigation of this radical was then reported in 1992 [59], but it was not 
until 1997 that its crystal structure was reported [52]. 



4.2.1. 1 

Electronic Structure 

EPR studies of 4d have, in contrast to other dithiazolyl radicals described 
above, shown considerable hyperfine coupling to substituents on the 
dithiazolyl ring; not only can N-coupling to the quinoxaline substituent be 
resolved, but also coupling to the four H atoms of the quinoxaline ring [59]. 
The second derivative EPR spectrum of 4d at 222 K is shown in Pig. 15 and, 
unlike the majority of thiazyi radicals of this class, shows considerable fine 
structure due to hyperfine interacfions with the benzo-fused ring. 



4.2.1. 2 

Crystal Structure 

The crystal structure of 4d (monoclinic, P2i) [52] exhibits a slipped 7i-stack 
motif (Pig. 15b). The molecules 7i-stack along the crystallographic n-axis with 
a separation (at 293 K) corresponding to the length of the crystallographic 
fl-axis, 3.7105(8) A. In addition to this n-n interaction, there is only one other 
contact within 4.0 A; an interstack S- • S contact at 3.84 A. 
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(a) 



(b) 

Fig. 15. (a) EPR spectrum of 4d illustrating hyperfine interactions to the N and H atoms of 
the fused ring and (b) its solid state structure 





4.2.1. 3 

Magnetic Properties of QDTA 

The regular 7i-stacked structure gives rise to some degree of band formation, 
despite the relatively long intra-stack separation which is beginning to 
approach the sum of the van der Waals radii of S perpendicular to the ring 
plane, 4.06 A [9]. At room temperature the magnetic moment corresponds to 
just 0.9 or 30% unpaired spins. On cooling the number of unpaired spins 
decreases steadily such that at 120 K all paramagnetism is quenched [53]. 



4.2.2 

The Benzo-1,2:4,5-Bis(Dithiazolyl) Biradicai, BBDTA 

This compound, 4e, was initially reported independently by Wudl [60] and 
Wolmershauser in 1986 [61]. Further work by Wudl provided a detailed EPR 
study [62]. This was followed in 1998 by further work by Wolmershauser who 
had partially oxidised the diradical to the radical cation and had isolated the 
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tetrachloroferrate salt [63]. More recently, Oakley and co-workers have 
published the structure of 4e and provided a detailed characterisation of the 
electronic structure [64]. 



4.2.2.1 

Electronic Structure 

The electronic structure of 4e has been the subject of some considerable 
debate. Theoretical calculations on the benzo-fused derivative BDTA indicate 
strong mixing of the bi SOMO of the dithiazolyl ring with the 7i-system of 
benzene [64]. This mixing also occurs in 4e, although the extent and nature of 
the intramolecular exchange between radicals becomes important. Theoretical 
calculations (cep-31 -t-g”*^* basis set) indicate that the singlet state lies slightly 
higher (175 cm“^) in energy than the triplet [64]. The EPR spectrum of 4e is 
both sample and solvent dependent [62, 64] with the signal intensity being 
very weak in relation to the sample concentration, suggesting the major 
species in solution is closed shell. Recent analysis by Oakley [64] suggests that 
considerable association of molecules occurs in solution, leading to either 
closed-shell dimers, or open-shell dimers in which the radical centres are well 
separated leading to a pair of equivalent, non-interacting doublets. Dissolution 
in toluene disfavours the association and a spectrum more closely resembling 
a pentet structure (superimposed with the triplet spectrum arising from the 
non-interacting spin doublets) is observed. The pentet structure is consistent 
with a diradical in which the intramolecular exchange coupling (/ex) between 
radical centres is greater than the hyperfine coupling, a^. 



4.2.2.2 

Crystal Structure 

The crystal structure of 4e is shown in Fig. 16. Diradical 4e crystallises in the 
monoclinic space group PIJc [64]. The structure consists of discrete 
molecules of 4e which form strongly slipped stacks (i.e., the ring plane is 
not orthogonal to the stacking direction) along the crystallographic n-axis. The 
molecular plane is tilted at 32.4° to this stacking direction. The mean 
interplanar separation along the stacking direction is 3.49 A, shorter than that 
observed for 4d. Additional inter-stack interactions between heterocyclic rings 
range from 3.12-3.74 A. 



4.2.2.3 

Electronic Properties 

The compound is essentially diamagnetic at room temperature, although a 
diamagnetic correction reveals some unpaired spin density, equivalent to ca. 
1% of Curie spins [64]. This is consistent with the much closer intra-stack 
separation 3.49 A (cf. 3.71 A for 4d). Although there is a considerable tilt of 
the molecular plane to the stacking axis, the better overlap evidently gives rise 
to a more well-developed band structure. An analysis of the conductivity 
indicates that it is a semiconductor with a band gap of 0.22 eV [64]. 
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Fig. 16. Crystal structure of 4e, viewed down the crystallographic b-axis 



4.2.3 

1,2,5-Thiadiazolo-Dithiazolyl-Pyrazine, TDP-DTA 



4.2.3. 1 

Electronic Structure 

The EPR spectrum of 4f comprises a triplet of pentets due to coupling to 
the dithiazolyl N and the two pyrazine N atoms (g = 2.0071, a^ = 9.50G, 
aN = 2.07G) [53]. The assignment of the smaller hyperfine coupling was made 
on the basis of spin density calculations at the MNDO level. These calculations 
indicate significant delocalisation of the SOMO onto the pyrazine substituent 
and considerably more than that observed for the benzo-derivative, 4a [53]. 



4.2.3.2 

Crystal Structure 

At room temperature 4f crystallises in the triclinic space group P-1 with two 
molecules in the unit cell [53]. The radicals form a slipped 7i-stack structure 
along the crystallographic c axis, but with significant slippage, such that 
equivalent atoms are separated by some 4.45 A, considerably beyond the sum 
of the van der Waals radii of S. 

On cooling to 150 K, a phase transition is observed. The low temperature 
phase has the same space group but doubles in size, with now four molecules in 
the unit cell [53] . In the low temperature phase the structure is built up of dimers 
which are linked through three S- • S contacts in the range 3.40-3.48 A. The 
change in structure upon cooling has profound effects on the magnetic 
behaviour. The phase change is reversible and on warming the high temperature 
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phase is reformed. However, the temperature at which the phase change occurs 
on cooling (broad transition between 150 K and 50 K, based on magnetic 
response) is different from the transition on warming (abrupt transition at 
175 K). The compound thus offers a region of bistability between 50 and 150 K 
in which the magnetic behaviour is temperature dependent. The crystal 
structures of both high and low temperature phases are shown in Fig. 17. 



4.2.33 

Magnetic Properties of TDP-DTA 

The room temperature moment for 4f is ca. 1.3 a little lower than that 
expected for an S = 1/2 paramagnet, but consistent with previous results on 
thiazyi radicals in which there is strong antiferromagnetic exchange [53]. On 
cooling below room temperature, the susceptibility passes through a broad 
maximum around 200 K consistent with the onset of short range antiferro- 
magnetic order. Around 150 K, y begins to drop markedly reaching a value of 
'^400 emu/mol at 50 K, equivalent to ca. 10% unpaired spins. Although the 
structural information clearly indicates a dimerisation process with quenching 
of the unpaired spins, the compound never becomes diamagnetic. Instead 
there is essentially a constant 10% of unpaired spins below 50 K which follow 
Curie-Weiss behaviour with a Weiss constant of -34 K. Oakley and co- 
workers indicate an antiferromagnetic state, although is unspecified [53]. 
Evidently further studies are warranted to examine the detailed magnetic 
behaviour of this radical. 

On warming the susceptibility remains essentially constant up to 170 K, but 
then begins to rise steeply. At 200 K the magnetic data coincides with the high 
temperature measurements once again [53]. 

This is, to our knowledge, the first reported observation of bistability in 
a thiazyi radical. The variation in susceptibility between the two phases is 




Fig. 17. High temperature (a) and low temperature (b) phases of 4f 
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most marked around 150 K at which point the susceptibility for the high 
temperature phase is approximately twice that of the low temperature 
phase. 



4.2.4 

Trithiatriazapentalenyl Radical, TTTA 

The radical 4g was first reported in 1989 by Wolmershauser [65]. Little had 
been reported since, apart from an EPR study [66], until 1999 when Fujita and 
Awaga reported that this radical exhibited bistability at room temperature 
[67]. Other studies have since been reported [68]. 



4.2.4.1 

Electronic Structure 

EPR studies by us [68] and other workers [66], coupled with DFT calculations 
[68], have provided a clear understanding of the electronic structure of 4g. The 
resolution of hyperfine coupling to both types of N atom in the dithiazolyl and 
thiadiazolyl ring has permitted an accurate analysis of the unpaired spin 
density [68]. The anisotropy in the g-values coupled with large hyperfine 
coupling to the high field gi component in the frozen glass spectra clearly 
indicate a 7i-based radical. An analysis of the spin density distribution 
provides a good correlation with theoretical calculations. These indicate that 
the majority of the unpaired spin density resides on the dithiazolyl ring, with 
approximately 22% of the unpaired spin density localised on the NSN unit 
(Fig. 18). 



4.2.4.2 

Crystal Structure 

This radical exhibits bistability and its structure has been investigated by 
several groups. The structure of the high temperature phase was originally 
reported by Wolmershauser [65]. Its structure at room temperature was 
reported by Fujita [67] and variable temperature studies down to the phase 
transition temperature at 225 K have been made by Rawson and co-workers 
[68]. The low temperature phase has been reported at room temperature by 
Fujita and at 150 K by Rawson [68]. 

The high temperature phase (monoclinic Pljc) exhibits a regular 
7i-stacking arrangement along the crystallographic fi-axis. The intermolecular 
separation coincides with the crystallographic fi-axis (3.71 A) and the radicals 




Fig. 18. Singly-occupied molecular orbital of 4g 
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are inclined to the stacking direction by 21°. In addition to the vertical 
stacking, there are a number of lateral S- • S and S- • N interactions in the 
range 3.05-3.64 A [68]. These are illustrated in Fig. 19. 

In the low temperature phase (triclinic, P-1), the radicals dimerise with 
intermolecular S- • S contacts in the range 3.24-3.32 A, with slightly closer 
contacts between S atoms at the dithiazolyl ring which bears the larger spin 
density. The disruption of the regular 7i-stack motif is accompanied by a 
slipping of the dimer pairs away from a perfect stacking arrangement. A set of 
lateral interactions in the range 2.92-3.65 A are observed. These are generally 
shorter than those observed in the high temperature phase [68]. The 
differences in the 7i-stacking arrangement of the high temperature and low 
temperature phases are shown in Fig. 20. 

X-ray diffraction studies indicate that the high temperature phase under- 
goes a phase transition to the low temperature phase on cooling below 225 K. 
However, the reverse process does not occur until 317 K. Thus, at room 
temperature the compound can exist in either of two interconvertible forms, 
i.e., it is bistable. 



4.2.43 

Magnetic Properties 

The magnetic properties of 4g are remarkable. The room temperature 
behaviour of 4g is dependent on sample history. If it has been annealed at 
340 K, then its magnetic behaviour is that of the high temperature phase [67]. 




Fig. 19. Lateral contacts in the high temperature phase of 4g, from ref. [68] 
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Fig. 20. 7i-stacking in (a) the high temperature and (b) low temperature phases of 4g, from 
ref. [68] 



Its room temperature moment is ca. 1.0 ^b> indicative of a very strong 
antiferromagnetically coupled regime. A fit of the data (350 to 230 K) to a one- 
dimensional antiferromagnetic chain yielded an intrachain exchange term, 
/intra> of “320 K, but required a large inter-chain interaction, / = -60 K [67]. 
The magnitude of / to /intra indicates that the one-dimensional model is not a 
good one, and the behaviour should be considered at least two-dimensionally. 
In conjunction with the low temperature crystallographic phase change, the 
susceptibility rapidly decreases below 230 K, and the compound becomes 
entirely diamagnetic by 150 K. The thermal hysteresis of 4g (measured by EPR 
spectroscopy) is shown in Fig. 21. 




Temperature / K 



Fig. 21. Thermal hysteresis observed in 4f, measured by EPR spectroscopy from ref. [68] 
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5 

Comment 



5.1 

Is It Possible to Make Room Temperature Organic Magnets? 

The magnetic behaviour of thiazyi radicals is characterised by strong 
antiferromagnetic exchange interactions between radical centres. The ex- 
change interactions can be so strong that short-range order is often observed 
in the region 10^-10^ K. Using a mean field theory approximation, the 
preparation of organic radicals that order at room temperature seems an 
achievable goal. In the majority of cases it would appear likely that they will he 
antiferromagnets since they exhibit strong antiferromagnetic exchange inter- 
actions. However, the lower symmetry of many molecular structures may 
favour the formation of spin-canted antiferromagnets, as has already been 
observed for P-NCC6F4CNSSN, 2 m. In these circumstances a small spontane- 
ous magnetisation may be observed. 

However, for this goal to he achieved, the dimerisation energy associated 
with many of these radicals must he overcome and the radicals must pack in 
such a way that magnetic exchange interactions can propagate throughout the 
solid. The control of solid state structure is a complex and arduous task, hut 
one which provides a goal worthy of achieving in the new millennium. 



5.2 

Bistability in Organic Radicals 

Two dithiazolyl radicals (4e and 4f) show thermal hysteresis; 4e exhibits 
thermal hysteresis between a paramagnetic state and an antiferromagnetic 
phase [53], whereas 4f exhibits thermal hysteresis between paramagnetic and 
diamagnetic states [67, 68]. A number of other dithiazolyl radicals have shown 
evidence that they undergo phase transitions, although the reversibility of 
their behaviour has not been reported. 

The bistability in these systems presumably arises through competition 
between the spin-density wave driven dimerisation process and lattice 
packing forces. In the case of derivatives of 4, the dimerisation energy is 
approximately zero and both open-shell (two doublet states) and the closed- 
shell configurations must have similar stabilities. The energy to interconvert 
between structures depends on the energy barrier for structural reorgani- 
sation. If the energy barrier is too high then two non-interconvertible 
polymorphs can be expected. Polymorphism has been observed in a number 
of thiazyi and selenazyl radicals, including 2f. In these cases the crystal 
packing in the different phases is substantially different and interconversion 
is not possible in the solid state. For smaller energy barriers, i.e., requiring 
minimum lattice reorganisation, then interconversion can be expected. In 
the case of 71 -stacked radicals the reorganisation energy can be expected to 
be low since it requires only minor displacements along the stacking 
direction. Contraction or expansion of the crystal lattice on cooling or 
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warming can lead to a variation in the relative stabilities of the two 
polymorphs. 

In the case of the dithiadiazolyl radicals, 2, the enthalpy for dimerisation 
clearly favours dimeric structures over monomeric ones and (in the majority 
of instances) the dimeric phase is only observed; the asymmetry in the 
potential energy wells inhibits bistability. Interestingly, a number of regular 
7i-stacked derivatives of 2 have recently been prepared [69] and it will be of 
interest to see whether, with appropriate control, bistability can be induced in 
dithiadiazolyl radicals. 

The observation of bistability at room temperature in 4f is evidently of 
significance in the development of organic data recording devices. The 
development of new thiazyl radicals with different ranges of bistability will be 
of significant interest. Moreover, the methods which might be employed to 
switch the bistability in such devices will lead to exciting new avenues in 
materials chemistry. 
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The discovery of Cgo and other members of the fullerene family opened new horizons in the 
design and synthesis of novel materials. Particularly exciting have been the results of the 
synthesis of Ceo-based materials, which exhibit transitions to ferro- and antiferro-magnetic 
ground states. In this chapter, we describe structural, electronic, and magnetic properties of 
fullerene salts in which magnetic transitions have been reported, i.e., the TDAE-Cso 
ferromagnet, the (NH 3 )K 3 C 6 o antiferromagnet, and the ACeo and Na 2 AC 6 o (A = K, Rb, Cs) 
polymers. The effects of different thermal treatments, changes of the interfullerene distance, 
and the role of orientational order/disorder of CJJq" ions on the magnetic properties are 
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LUMO lowest unoccupied molecular orbital 

/rSR muon spin relaxation 

SC superconducting 

SDW spin density wave 

ZF zero field 

1 

Introduction 

Organic materials are molecular solids in which intermolecular distances are 
commonly much longer than the interatomic distances found in inorganic 
metallic and ferromagnetic solids. The molecules are normally diamagnetic, 
closed-shell species and as such they do not usually display interesting 
electronic and magnetic properties. Fortunately, in the past decade, there has 
been enormous progress in the development of molecular-based magnetic 
materials [1] with large molecular units as building blocks. In these, the 
magnetic moments are derived exclusively from unpaired p-electrons, as 
opposed to the more commonly encountered d- and /-electron systems. 
Basically, one can divide organic molecules which are open-shell entities into 
two groups containing either neutral, localised unpaired electron-containing 
functional groups (free electrons) or charged species (radical ions). The 
nitroxyls are examples of members of the former category. Radical ions result 
by electron transfer from the 7i-orbitals of an electron-rich molecule or 
transfer of an electron to the n molecular orbital of an electron-deficient 
molecule. In this chapter, we review new magnetic materials of the donor- 
acceptor type, based on the buckminsterfullerene, Ceo, molecule. 

A new chapter in organic chemistry and condensed matter physics was 
opened recently with the discovery of Ceo (Fig. la) and other fullerenes [2]. In 
Ceo, 60 carbon atoms are arranged in the form of a truncated icosahedron 
(soccer ball). This molecule has the high symmetry of an icosahedron (point 
group 4) with six five-fold, ten three-fold axes, fifteen two-fold symmetry axes, 
and inversion symmetry. From an organic chemist’s perspective, Cso repre- 
sents an entirely new type of organic molecule and opens new horizons in the 
design and synthesis of novel compounds. On the other hand to a physicist. 




Fig. 1. Molecular structures of (a) C^o and (b) tetrakis-dimethyl-aminoethylene (TDAE) 
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the Ceo molecule represents a highly symmetric, stable molecular analogue of a 
very large atom, which in the solid state can adopt a crystal structure with 
relatively large interstitial sites that can be filled with atomic or molecular host 
units (Fig. 2). 

Even before its production in macroscopic quantities, the electronic 
structure of molecular Ceo had attracted considerable theoretical interest. 
Due to the small overlap between the electron clouds of neighbouring 
molecules, the electronic structure of solid Ceo can be described as a relatively 
small perturbation of its molecular orbital levels scheme [3]. Pristine solid Ceo 
is an insulator and there is a =1.5 eV gap between the (filled) highest occupied 
molecular orbital (HOMO) and the lowest unoccupied molecular orbital 
(LUMO) (Fig. 3) [4]. The LUMO of tiu symmetry is triply degenerate and can 
accept up to six electrons per Ceo molecule. Intercalated Ceo salts should thus 
become metallic if fewer than 6 electrons are injected into the tiu-derived 
conduction band. This idea stimulated a great amount of work on intercalated 
Ceo compounds, resulting in some spectacular results. Certainly one of the 
most important achievements has been the discovery of superconductivity 
with relatively high Tq in alkali fullerides with stoichiometry A^C^o, where 
(A = Na, K, Rb, Cs) [5]. 

The importance of Ceo as a new building block for the design of novel 
magnetic materials lies with its great flexibility: 

(a) there is vast number of ways of doping pristine Ceo that can result in 
materials with exciting electronic and magnetic properties; 




Fig. 2. The cryolite structural type adopted by the K3C60 fulleride. The shaded and white 
spheres represent potassium ions residing in the octahedral and tetrahedral sites, 
respectively 
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Fig. 3. Molecular orbital level scheme of Ceo (adapted from ref. [4]) 

(b) the Ceo-Ceo distances can be easily controlled through the use of 
appropriate ions with different sizes (changing K"^ to Cs"^, for instance), 
of co-ordination complexes (co-intercalation of alkali ions and neutral 
NH 3 ) and of large organic donors; 

(c) different thermal treatments of the samples can lead to stabilisation of 
competing phases characterised by the formation of interfullerene bridging 
C — C bonds; and 

(d) the electronic and magnetic properties of the materials are sensitive to the 
presence of orientational order/disorder of Ceo units. 

Most attempts to synthesise intercalation compounds of Ceo have employed 
the introduction of dopants into the empty interstitial space between the Ceo 
units (Fig. 2). In some cases - the building blocks are shown in Fig. 1 - such 



The ionic radii of the alkali ions are: 0.69 A (Li"^), 0.95 A (Na^), 1.33 A (K^), 1.48 A (Rb^), 
and 1.69 A (Cs"^). These should be compared with the size of the tetrahedral (1.12 A) and 
octahedral (2.06 A) interstices. 
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doping led to very interesting magnetic behaviour. In this review, we focus our 
attention only to those fulleride systems in which magnetic transitions have 
been encountered, i.e., the TDAE-Ceo ferromagnet, the (NH 3 )K 3 C 6 o antiferro- 
magnet, and the ACeo and Na 2 AC 6 o (A = K, Rh, Cs) polymers. 



2 

Ferromagnetism in TDAE-Cgo 



The purpose of this section is to review the magnetic properties of one of the 
most interesting new magnetic organic compounds of the past decade, namely 
the ferromagnetic fulleride tetrakis-dimethylamino-ethylene-Ceo (TDAE-Ceo) 
(Eig. 4), first discovered in 1991 by Wudl and co-workers at the University of 
California, Santa Barbara [6]. At the time of discovery, the Curie temperature, 
Tc = 16 K of TDAE-Cso was one order of magnitude higher than the existing 
record [7] and brought the research field of p-electron ferromagnetism from 
the realms of the esoteric into the mainstream. 

When the Santa Barbara group measured the magnetic properties of 
TDAE-Ceo in powder form, a steep rise in the real part of the ac susceptibility 
(Fig. 5) was evident below 16 K, implying the occurrence of a ferromagnetic- 
type transition at a surprisingly high transition temperature, Tq- 16 K for a 
compound comprising only elements of the first row (C, N, H) of the periodic 
table. Since then, numerous experiments have been performed on both 
powders and single crystalline samples of TDAE-Ceo, leading to a compre- 
hensive understanding of its intriguing behaviour. 

It was soon realised that the ferromagnetic transition in TDAE-Ceo was not 
of a conventional type, which could be simply described within some standard 
theory. Some of the puzzling properties of TDAE-Ceo were described in detail 
in the original report on its synthesis and characterisation [6]: 

(i) despite the sharp increase of the magnetisation below Tc, its temperature 
dependence did not follow conventional mean field theory behaviour; 

(ii) within experimental error, no hysteresis was observed between measure- 
ments under cooling and heating protocols; 

(iii) similarly, the coercive field and the remnant magnetisation were found to 
be zero; and 

(iv) the spontaneous magnetisation was very small, leading to a low- 
temperature saturation moment of 0.11 /rg/Cgo- 
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Fig. 4. Schematic diagram of the TDAE reduction 
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Fig. 5. Temperature dependence of the real part of the ac susceptibility (v = 3 kHz) of a 
single crystal of TDAE-Ceo; showing the onset of the ferromagnetic transition at Tc = 16 K 

In order to explain these unusual magnetic properties, a variety of models 
appeared in the early literature, describing TDAE-Ceo as an itinerant 
ferromagnet [6], a superparamagnet [8], a spin-glass [9], or a weak 
ferromagnet [10]. 

When TDAE (a liquid at room temperature) is added to C^o in solution, the 
resulting salt crystallises as small particles, which can be removed by filtration. 
TDAE+CgQ is a charge-transfer salt with a monoclinic structure, different from 
the high-symmetry cubic structures of pristine Cso and its alkali intercalated 
derivatives [11]. The crystal structure of TDAE-Ceo was first determined by 
Stephens et al. [12] to be monoclinic (space group C2lm) with two formula 
units per unit cell. However, a subsequent X-ray single-crystal structural 
determination [10] identified the existence of a superlattice, arising from the 
presence of symmetry-inequivalent TDAE"^ units in the unit cell. Two subcells 
are present, with the TDAE"^ ions shifting along the b axis by '^0.02 A in 
opposite directions and resulting in doubling of the unit cell along the c axis 
(Eig. 6). The structure at room temperature is then monoclinic (space group 
C2/c) with cell dimensions a = 15.858(2) A, = 12.998(2) A, c = 19.987(2) A, 
P = 93.37(3)° and four formula units per unit cell. The TDAE units reside in 
the (1/2, 1/2 + 5, 3/4), (1/2, 1/2-3, 1/4), (0, <5, 3/4), and (0, -d, 1/4) sites with 
3 = 0.002 and the C=C double bond parallel to the c axis, whereas the Ceo 
molecules are centred at the (0, 1/2, 0), (0, 1/2, 1/2), (1/2, 0, 0), and (1/2, 0, 1/2) 
sites. The resulting Ceo-Ceo centre-to-centre distances are shortest along the 
crystal c axis (9.99 A at room temperature as compared to 10.24 A in solid Ceo 
[11]). At 80 K, the centre-to-centre distances along the c-axis contract further 
to 9.87 A. The unusually short distances between adjacent Ceo units are crucial 
in leading to coupling between spins localised on neighbouring Cgg ions and 
the formation of the ferromagnetically correlated spin state. 
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Fig. 6. Schematic diagram of the monoclinic unit cell of TDAE-Cgo (space group C2/c) 
consisting of two subcells stacked along the c axis. The arrows indicate the small shifts of the 
TDAE units along the h direction 



Since charge transfer from TDAE to C^o leads to a partially full molecular 
orbital (Fig. 3) and because of the close contacts of the Ceo units, it might be 
expected that TDAE-Ceo is a metal with a highly anisotropic electronic 
structure, displaying itinerant ferromagnetism [6]. However, powder infrared 
absorption spectroscopy revealed that TDAE-Ceo displays no Drude tail and 
exhibits no absorption at low frequencies [13], implying that the material is 
most probably an insulator. Microwave absorption measurements [14] also 
supported the insulating character of TDAE-Ceo- The issue was finally settled 
unambiguously by dc and ac conductivity measurements on single crystals of 
TDAE-Ceo as a function of temperature [15]. These experiments ruled out the 
possibility of an itinerant ferromagnetic ground state below Tc = 16 K. 

The first unambiguous evidence for long-range magnetic order came from 
zero-field muon spin relaxation (ZF-p'^SR) experiments performed on 
TDAE-Ceo powder [16]. p'^SR spectroscopy is an extremely powerful technique 
in cases of small-moment magnetism and in all instances where magnetic 
order is of a random, very short-range, spatially inhomogeneous or incom- 
mensurate nature [17]. The are implanted into the solid sample and after 
they come to rest at an interstitial site, they act as highly sensitive microscopic 
local magnetic probes. In the presence of local magnetic fields, they will 
precess with a frequency = {y^/2n){B^), where = 13.55 kHz/G. In the 
absence of an applied external field (ZF), the appearance of a precession 
signals the onset of an ordering (ferromagnetic or antiferromagnetic) 
transition. Moreover, application of a magnetic field parallel to the initial 
muon spin polarisation (LF) allows the decoupling of the spin from the 
static internal fields. 

In well-annealed powder TDAE-Ceo, a heavily damped oscillating signal 
[v^ = 0.92(2) MHz at 3.2 K] was observed in zero field (Fig. 7), providing 
definite proof of the existence of long-range magnetic order below 16.1 K. The 
strong muon spin relaxation - whose quasi-static nature was confirmed by LF 




136 



D. Arcon • K. Prassides 




Time (psec) 



Fig. 7. Evolution of the zero-field (ZF) spin polarisation, P^{t) between 3.2 and 16.8 K for 
TDAE-Ceo 



measurements - found in ZF implies the presence of substantial spatial 
disorder and inhomogeneity effects. The local field distribution, 
= 48(2) G is only smaller than (B^) = 68(1) G by a factor of 1.4. In 
addition, the temperature dependence of (Fig. 8) which also mirrors that of 
the magnetisation measured at 100 G behaves on approaching Tc in a different 
manner from that expected by the mean-field treatment of a 3D Heisenberg 
exchange model. It follows well Bloch’s T^^^-law and thus it appears that 
spin wave (magnon) excitations dominate the magnetic behaviour, even at 
temperatures close to Tq. Alternatively, the magnon-like behaviour may reflect 
a more complicated physical picture of the system, where the intrinsic 
orientational disorder associated with the fullerene molecules gives rise to 
a broad distribution of exchange constants (and ordering temperatures) 
[17, 18]. 




Fig. 8. Temperature dependence of the pA precession frequency, (solid circles) in 
TDAE-Ceo powder. The solid line is a fit to Bloch’s 7^^^-law 






Magnetism in Fulierene Derivatives 



137 



thus provided clear evidence for the long-range ordering below 
Tq = 16 K, but many details were still missing. Additional information has 
been obtained from electron spin resonance (ESR) experiments. The ESR line 
width (AHpp) in the paramagnetic phase gives information about the spin 
dynamics in the system under investigation. Eor TDAE-Ceo, a dramatic 
decrease in AHpp was detected below Tq = 160 K (Fig. 9). Above Tq, the line 
width is determined by the dipolar interactions between the neighbouring 
spins, while below Tq the exchange interactions become important and rapid 
fluctuations of the electronic spins on the ESR timescale (10“^° s) cause the 
narrowing of the ESR line (so-called exchange narrowing). The sharp decrease 
of the ESR line width is related to the slowing down of the rotational motion of 
the Ceo units, as also observed by ^^C-NMR spectroscopy [19]. Therefore Tq is 
the orientational ordering temperature of the Cgg ions. In this sense, TDAE-Ceo 
is a highly interesting material as the orientational ordering of Ceo units 
directly influences the magnetic interactions [20, 21] through the different 
relative orientations of the molecular orbitals and different overlap of the 
neighbouring electronic wave functions. This relationship between the 
orientational and spin degrees of freedom has been subsequently confirmed 
by the isolation and characterisation of two structural modifications of TDAE- 
Ceo single crystals, designated as a- (ferromagnetic) and a'- (non-magnetic) 
phases [22, 23]. The two structural forms essentially differ only in the degree of 
the orientational order of the Cgg ions [24]. Annealing of the a'-form at or 
above room temperature leads to a transformation to the ferromagnetic 
a-modification. 

Orientational disorder effects, associated with the quasi-spherical nature of 
the fulierene molecules are abundant in Ceo-based compounds. Orientational 
order in pristine Ceo can be understood in terms of the combination of van der 
Waals and electrostatic interactions which result in the most electron-poor 
regions (the pentagonal faces) of the molecules facing the most electron-rich 
regions (the higher bond-order 6:6 bonds) of adjacent molecules [25]. 




Fig. 9. Temperature dependence of the ESR (9.6 GHz) line width in TDAE-Cso single crystal. 
The solid line is a guide to the eye 
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However, disorder is still present, as a second “contact motif’ is also present 
involving hexagons facing 6:6 bonds of neighbouring molecules [103]. Although 
the orientational potential in TDAE-Ceo is much more complicated, we expect 
that similar structural motifs should be found and for optimal contact between 
adjacent Cgg ions the ferromagnetic interactions are favoured. In fact in a 
recent study, Kambe et al. [26] studied in detail the effects of annealing on the 
magnetic and structural properties. It was found that well-annealed samples 
undergo a structural phase transition at around 180 K - most probably 
associated with the Cgg orientational ordering previously detected by ^^C-NMR 
- and that in this case the saturated magnetisation increases to 0.9(1) 
ion. On the other hand, the as-grown crystals undergo no structural phase 
transition at least down to 30 K and do not order ferromagnetically. This 
underlines the importance of orientational order/disorder effects on the 
magnetic properties of TDAE-Ceo compound. A theoretical model has been 
proposed in which the effective ferromagnetic exchange coupling strength, Jeff, 
shows pronounced variation with varying relative orientation of the Cgg ions 
[21]. This model naturally explains the controversy raised by a number of 
contradictory reports on the magnetic properties. 

Magnetic resonance experiments on magnetically ordered phases differ 
from ordinary ESR on paramagnetic samples in which only one electronic spin 
interacting with weak local fields is excited. The presence of large exchange 
fields in magnetic systems, typically orders of magnitude larger that the 
applied field, causes a coherent precession of all the electronic spins. The 
entire magnetisation of the sample is thus excited and the precession 
frequency depends on the magnitudes of the exchange, anisotropy, and 
demagnetisation fields. The relationship between the resonance frequency and 
the resonance field becomes strongly non-linear and depends on the nature of 
the magnetic ordering [27]. While in simple ferromagnets only one resonance 
mode is predicted theoretically, in antiferromagnets and weak ferromagnets 
two resonance modes should be found. Eor instance, for both uniaxial and 
weak ferromagnets, when the external magnetic field is perpendicular to the 
easy axis, a resonance mode with a dip in the resonance field-resonance 
frequency relation at a resonance field equal to the anisotropy field is found. 
However, in weak ferromagnets, an additional high-frequency antiferromag- 
netic type mode is also present. 

On the other hand, the magnetic resonance in spin-glasses is very 
complicated and strongly dependent on thermal history [28]. Eor instance, 
due to hysteresis effects, the resonance line position under field- cooling 
conditions depends on whether the measured field is parallel or antiparallel to 
the cooling field. Superparamagnetic materials [29] behave in a similar fashion 
to normal paramagnets and the resonance position of the observed lines 
depends linearly on the resonance frequency. As a result, the (anti) ferromag- 
netic resonance technique provides an extremely sensitive tool for the 
determination of the ground state and the microscopic parameters in a variety 
of magnetic materials. The theoretical curves depicting the resonance 
frequency-resonance field relationship in ferromagnets, antiferromagnets 
and weak ferromagnets are summarised in Eig. 10. 
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Fig. 10. Resonance frequency-resonance field relationship for (a) ferromagnets, (b) anti- 
ferromagnets, and (c) weak ferromagnets. Here Hg is the exchange field, Hk is the 
anisotropy field, and Hdm is the Dzyloshinskii-Moriya field 



Clear evidence for the presence of a ferromagnetic resonance in TDAE-Ceo 
was found in the radiofrequency region [30]. A single resonant line was 
observed in high-frequency experiments performed between 110 MHz and 
220 GHz [31]. Below 110 MHz, a new line emerges at zero field which shifts to 
higher resonant fields with decreasing resonant frequency. The two resonant 
lines, high and low field, merge together below 50 MHz. The dependence of the 
resonant frequency on the resonant field is strongly non-linear and shows a 
characteristic dip at H = 30 G (Fig. 11), as expected for both uniaxial and 
weak ferromagnets (Fig. 10). Since no other resonant modes are observed at 
higher frequencies, the possibility of a weak ferromagnetic ground state in 
TDAF-Ceo can be ruled out. The strongly non-linear behaviour (Fig. 11a) also 
eliminates the possibility of superparamagnetism and leads to the description 
of the magnetic ground state in TDAF-Cgo single crystals as that of a normal 
Heisenberg ferromagnet. Agreement between theory and experiment becomes 
quantitative, if the demagnetising field effects are also taken into account. This 
leads to values of the anisotropy field. Hr = 29 G and the demagnetising field, 
Hdem - ~39 G. The non-linear dependence of the resonance field on the 
resonance frequency disappears above the transition temperature in the 
paramagnetic phase (Fig. 11b). 

In conclusion, we stress that TDAF-Cso is thus far the only authenticated 
ferromagnetic material based on Geo- Despite considerable efforts, doping with 
other organic donors invariably leads only to paramagnetic materials. It 
appears that the fortuitous combination of a number of factors, including 
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Fig. 11 . Dependence of the resonance frequency on the resonance field (fl||H) in TDAE-Cs, 
single crystal: (a) at T = 5 K (<Tc) and (b) at T = 20 K (>Tc) 



orientational ordering of Cgg ions, short Ceo-Ceo distances along the c-axis, and 
the role of TDAE"^ ions results in the formation of a ferromagnet with such an 
unusually high Curie temperature, unique among currently known non- 
polymeric, purely organic systems. 



3 

Antiferromagnetic Ordering in (NH3)K3Ceo 

Intercalation of solid Ceo with alkali metals affords a series of salts with 
composition (A = alkali metal, 0 < x < 6). Compounds with stoichio- 

metry AjCso are invariably metallic, becoming superconducting with critical 
temperatures, Tq as high as 33 K at ambient pressure [32] with Tq increasing 
monotonically as the interfullerene separation, d increases [5, 33]. This can be 
rationalised in terms of increasing density-of-states at the Fermi level, N(sp) 
with increasing d, resulting from the decrease in the overlap between the 
molecules that leads to band narrowing. As a consequence, in order to obtain 
high-Tc fullerides, large interfullerene spacings are needed. Ammoniation has 
proven an excellent method to achieve large expansions of the unit cells of 
fulleride salts, as neutral NH3 molecules co-ordinate to the alkali ions, leading 
to large effective radii for the resulting (NHs)^^ species. Successful ammo- 
niation of superconducting Na2CsC6o (Tc = 12 K) has afforded the cubic 
(NH3)4Na2CsC6o salt with Tq of 29.6 K [34]. On the other hand, reaction of 
K3C60 (Tc = 19 K) with NH3 led to the isolation of the orthorhombic 
(NH3)K3C6 o phase which is non-superconducting [35]. However, application 
of a pressure >1 GPa leads to recovery of superconductivity in (NH3)K3C6 o 
with Tc = 28 K [36]. 
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Structural characterisation [35, 37] of (NH 3 )K 3 C 6 o at room temperature 
revealed an anisotropic expansion of the fulleride array leading to a face- 
centred orthorhombic cell with lattice constants a = 14.917(10) A, 
b = 14.971(10) A, c = 13.692(4) A (space group Fmmm). The observed 
increase in the interball separation and the reduction in crystal symmetry 
should have important consequences for electron localisation and the loss 
of superconductivity. Below 150 K, a phase transition occurs to another 
orthorhombic structure (space group Fddd) in which the unit cell dimensions 
double along all three crystallographic directions [37]. The origin of the 
observed superstructure was originally identified with the antiferroelectric 
ordering of the K'^-NH 3 pairs, which reside in the octahedral sites [37]. Its 
signature was also evident in neutron inelastic scattering measurements in 
which a librational mode at '^6.5 meV, implying freezing of the ammonia 
rotation, appears below 150 K [38]. However, recent neutron diffraction 
experiments have revealed additional subtleties in the adopted structure 
arising from orientational ordering of the Ceo units, which adopt two different 
orientations related by a 90° rotation about the c orthorhombic axis and order 
along the a axis (antiferrorotative ordering. Fig. 12) [39]. Structural work at 
elevated pressures also showed that above '^1 GPa the system is isostructural 
with the ambient-pressure low- temperature phase [39]. With increasing 
pressure, (NH 3 )K 3 C 6 o becomes increasingly more anisotropic with an in- 
creased size difference between the lattice constants a and b. 

The position of the ^^C-NMR line (195 ± 2 ppm) in (NH 3 )K 3 C 6 o [40, 41] 
falls within the range of other Cgg^ ions and is only slightly larger than that of 




Fig. 12. Low-temperature orthorhombic crystal structure of (NH3)K3C6o (space group Fddd), 
projected on the [ 110 ] plane. The observed superstructure arises from the antiferrorotative 
order of the C^o units and the antiferroelectric order of the K'^-NH3 pairs 
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K 3 C 60 (187 ppm). The narrow ^^C-NMR line at room temperature implies 
quasi-isotropic fast reorientation on the NMR timescale of the Cgg ions. 
The ^^C-NMR spectra broaden below 150 K due to the freezing out of the 
molecular motions. molecular dynamics has been studied in detail by 
the low-energy inelastic neutron scattering measurements as a function of the 
scattering vector Q and in the temperature range between 310 K and 30 K [42]. 
The CgQ ions were found to execute small-amplitude librations about their 
equilibrium positions, giving rise to well-defined peaks near 3.1 meV at low 
temperature (Fig. 13). The librational energy is smaller than in the parent 
KsCgo salt (4.04 meV at 12 K), and the peaks are much broader, reflecting a 
weaker and more anisotropic orientational potential than that encountered in 
the parent KaCgo [42]. 

The electronic and magnetic properties of (NH 3 )K 3 Cgo were measured 
with the NMR, ESR, and pSR techniques. Measurements of spin-lattice 
relaxation rate, 1/Ti, showed that (NH 3 )K 3 Cgo is metallic above 40 K and 
the Korringa law is followed between 40 and 100 K. However below 40 K, 
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Fig. 13. Representative NIS spectra of (ND 3 )K 3 Cgo in the temperature range 30-310 K at 
constant Q = 5.5 
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Fig. 14. Temperature dependence of the ESR intensity in NH3K3C60 



NH3K3CgQ powder 



• • 



• • 



1/Ti suddenly drops corroborating the occurrence of a phase transition 
[40, 41], 

More information about the low temperature transition at 40 K can be 
deduced from the ESR results which showed that (NH 3 )K 3 C 6 q is a narrow 
band metal [40, 41, 43, 44]. The room temperature spin susceptibility, xs 
determined by double integration of the integrated derivative ESR spectra, is 
5.0 X lO”'^ emu/mol. The value of xs is characteristic for the metallic state of 
fulleride salts and is comparable with that measured in K 3 C 60 [45]. On 
cooling, the ESR spin susceptibility first slightly increases down to 40 K 
(Eig. 14). The increase in xs is rather exceptional for the A 3 C 60 systems, 
where the susceptibility normally decreases due to the lattice contraction and 
thus the reduction of the density of states at the Eermi level [46]. The 
variation of xs with temperature is difficult to explain simply in terms of 
metallic or localised spin behaviour. The increase of the susceptibility is too 
weak to originate from Curie- Weiss dependence [47]. However, the same 
type of behaviour was found in polymeric RbCeo and has been ascribed to 
the strong antiferromagnetic correlations in a narrow band metal system. 
Below 40 K, the ESR susceptibility suddenly drops by more than an order of 
magnitude. Such a drop in xs could be a result of antiferromagnetic ordering, 
transition to a spin-density wave (SDW) or even to non-magnetic charge- 
density wave (CDW). It is of interest to note that the ESR line width, which 
should increase in the vicinity of the transition due to the increased spectral 
density of the local field fluctuations, does not show any anomaly. If the 
ground state is characterised by long range ordering of electronic spins, 
antiferromagnetic or SDW type, then antiferromagnetic resonance should be 
detected in high-frequency ESR experiments. Magnetic resonance experi- 
ments at frequencies between 9 and 225 GHz [44, 47] indeed identified an 
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additional resonance line, which has been attributed to the antiferromagnetic 
resonance. Below 40 K, the resonance line broadens and shifts as the 
magnetic order develops. For the antiferromagnetic resonance in powdered 
samples, the line width, AHpp should decrease with increasing resonance 
field. Ho {Ho » Hsf, where ffsp is the spin-flop field) according to the 
relationship: AHpp oc Hjp/Ho. However, while AHpp between 9 GHz 

(Ho = 0.3 T) and 75 GHz (Hq = 2.7 T) narrows, for frequencies between 
75 GHz and 225 GHz, the line width depends only weakly on external field 
and does not follow the 1 /Hq field dependence. One of the reasons could be 
that the line width narrowing is compensated by the g-factor anisotropy 
broadening effects [47]. Experiments on single crystals of NH 3 K 3 C 60 should 
certainly help to clarify the nature of the observed line. 

Even though the early ESR and NMR measurements of the electronic and 
magnetic properties have shown that (NH 3 )K 3 C 6 o is a narrow band metal 
which exhibits a transition to an insulating ground state at about 40 K, the 
nature of the low temperature state remained controversial as the experimental 
results were unable to distinguish between a magnetic or a non-magnetic 
ground state. The problem of the existence of non-zero internal local magnetic 
fields at low temperature in (NH 3 )K 3 C 6 o and its deuterated analogue, 
(ND 3 )K 3 C 6 o was addressed by using 100% spin polarised positive muons 
in the absence of external fields [48, 49]. At temperatures higher than 
40 K, only a non-oscillating signal is evident. In this temperature range, the 
spin relaxation is determined by the fluctuations of the nuclear moments, 
which appear frozen into a disordered spin configuration. In addition, the 
spin is relaxed by the rapidly fluctuating electron spins, which give rise to the 
Lorentzian relaxation component with 2 = 0.030(2) /is~^ at 100 K. Cooling- 




Fig. 15. Evolution of the zero-field (ZF) spin polarisation, P^(t) at 14 and 40 K for 
(ND3)K3C60- The solid lines through the data are fits to the functions of Eqs. 1 and 2, 
respectively 
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Fig. 16. Temperature dependence of the ZF muon precession frequency, (open circles) in 
(ND3)K3C6 o; the solid line is a fit of the data to the critical law, Eq. 2. The inset shows the 
temperature dependence of the volume fraction (open circles) of the magnetically ordered 
component in (ND 3 )K 3 C 6 o; the broken line represents a guide to the eye 

down leads to a slowing-down of the electron spin dynamics within the 
paramagnetic domains with X reaching a value of 0.046(1) /ts~^ at 40 K [49]. 

Below 40 K (Fig. 15), the shape of the time-dependent /t'^SR spectra 
changes, as a short-lived oscillating component, whose depolarisation grad- 
ually increases with decreasing temperature, appears. The observation of a 
precession signal in zero external field is clear evidence of the existence of 
coherent ordering of the electronic spins and indicates unambiguously the 
onset of long-range antiferromagnetic order. The precession frequency, is 
0.564(4) MHz at 14 K, corresponding to a static local field at the muon site, 
(B^) = 41.6(3) G. In addition, the depolarisation rate, X has a value of 
2.14(4) at 14 K, implying a distribution of local fields with a width 
= 25.1(4) G, only smaller than (B^) by a factor of 1.7. The muons 
thus experience a local field with large spatial inhomogeneities, which may be 
due to a number of physical reasons, including orientational disorder effects of 
the fulierene molecules. The temperature evolution of the muon frequency, v^, 
is described by the equation: = Vo[l — (T/Tn)]^, with Vq = 0.64(1) MHz, 

Tn = 36.5(3) K, and jS - 0.32(3) (Fig. 16). The critical exponent, /I, is close to 
the value expected for a conventional 3D Heisenberg antiferromagnet 
(/I = 0.36). As the y*' site is not known for a powder sample, a value of the 
magnetic moment cannot be extracted directly. However, we note that Vq is 
1.4 times smaller than the frequency observed in (TDAE)C6o, in which the 
magnetic moment is '^1 /iB/molecule [26]. Assuming similar stopping sites for 
the two samples, we find /r(T = 0 K) is in the order of 0.7 pe/molecule. 

Data collected at 14 K in a 200 G longitudinal field reveal a complete 
recovery of the asymmetry. As the effect of applied longitudinal fields is to 
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allow the depolarisation due to dynamic or fluctuating moments to be 
decoupled from that due to static components, we conclude that the origin of 
the observed relaxation in zero field is of quasi-static nature. 

Recent ^^C-, ^H-, and ^^K-NMR measurements have also confirmed the 
antiferromagnetic nature of the low-temperature ground state in (NH 3 )K 3 C 6 o 
[50]. However, in disagreement with earlier ESR and NMR results, the authors 
conclude that that even above 40 K, the system can be described as an S = 1/2 
localised spin system, and not as a highly-correlated metal. 

In conclusion, (NH 3 )K 3 C 6 o and its deuterated analogue show a transition 
below '^37 K to a long-range ordered antiferromagnetic state, characterised by 
considerable spatially inhomogeneous effects. The suppression of supercon- 
ductivity is thus associated with effects of magnetic origin, providing an 
important analogy with the well-established phenomenology in high-Tc and 
organic superconductors. Thus, the increased interfullerene separation to- 
gether with the symmetry reduction in (NH 3 )K 3 C 6 o have important conse- 
quences for the electronic properties of the material. The overlap between the 
molecules decreases substantially, leading to a reduced bandwidth, W, and an 
increased (UIW) ratio, which for fixed band filling drive the system to an 
antiferromagnetic Mott insulating state [51]. Then the recovery of supercon- 
ductivity at pressures higher than 1 GPa (Tc = 28 K) [36, 39] should 
presumably be associated with the suppression of the magnetic transition. 
Fig. 17 displays a schematic electronic phase diagram for CgQ compounds as a 
function of the volume per CgQ^ anion [or equivalently increasing (UIW) ratio]. 




Fig. 17. Schematic electronic phase diagram of Cj]] compounds, showing the approximate 
location of the metal (superconductor)-insulator phase boundary. The open symbols are 
literature values of Tc for a variety of superconducting fullerides, while the solid symbols 
mark (ambient pressure) and Tc (>1 GPa) of (NH3)K3Cgo 
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Future studies of the electronic properties of systems lying on either side of the 
metal (superconductor)-insulator phase boundary promise to lead to intrigu- 
ing results. 

4 

Magnetism in Polymeric Fullerides 



4.1 

Magnetic Properties of the ACeo Family 

Stable phases at the x = 1 point in the phase diagram of A^Ceo (A = K, Rb, Cs) 
were not identified until some time after the superconducting x = 3 and 
insulating x = 4 and 6 phases were discovered and characterised [52]. The first 
indications of stable x = 1 phases came from Raman [53], photoemission [54], 
X-ray diffraction [55], and ESR experiments [56]. Soon afterwards it was 
observed that RbCgo can be either conducting or insulating below room 
temperature, depending on the thermal history. Depending on the cooling 
protocols used, a number of distinct phases has now been stabilised (Fig. 18). 

At high temperature (above 400 K), ACeo is fee with alkali ions residing in 
the octahedral sites [55]. Below 400 K, the most stable form was originally 
shown to adopt orthorhombic crystal symmetry [57, 58] with pronounced 
quasi-one-dimensional character and unusually close contacts between the 
fulleride ions; partial covalent bonding between fullerenes along the chains 
and their consequent molecular deformation from icosahedral symmetry have 
been postulated to account for these observations. Structural characterisation 
was originally achieved by synchrotron X-ray powder diffraction [58] studies, 
which established a cell of dimensions, a = 9.11 A, h = 10.10 A, and 
e - 14.21 A (space group Pmnn) for RbCeo at room temperature with fulleride 
linkages along the polymer chains achieved by a [2-1- 2] cycloaddition 
mechanism, in analogy with the photopolymerisation reaction of pristine 
Ceo- The contraction along the polymerisation a axis is '^0.9 A, compared to 
the centre-to-centre interfullerene distances of about 10 A, encountered in 
monomeric fullerenes and their derivatives. However, subsequent structural 
work on KCeo and ACeo (A = Rb, Cs) [59, 60] revealed that they possess 
different chain orientations about their axes, which are described by distinct 
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Fig. 18. Sequences of phases observed in ACgo fullerides as a function of temperature and 
cooling rates 
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space groups Pmnn (orthorhombic) and 72/m (monoclinic), respectively 
(Fig. 19). 

A medium quenching rate to just below room temperature leads to the 
stabilisation of a low-symmetry insulating structure which has been charac- 
terised by synchrotron X-ray diffraction studies to comprise (Cgg)2 dimers, 
bridged by single C — C bonds [61]. The RbCeo dimer phase crystallises in the 
monoclinic space group Plja with lattice parameters a = 17.141(5) A, 
b = 9.929(5) A, c = 19.277(5) A, // = 124.4° at 220 K (Fig. 20) and transforms 
back to the polymer phase on heating. The chain axis is now along b and is 
somewhat longer than that found in the RbCeo polymer, implying a smaller 
contraction of '^0.7 A from the interfullerene separations found in monomeric 
fullerenes. This is consistent with the presence of a single bridging C — C bond, 
instead of the four-membered rings of the polymers. Another metastable 
conducting phase of RbCeo and CsCeo was reported with very rapid quenching 
by immersing the high-temperature sample in liquid nitrogen [62]. In this 
case, the 3D isotropic character of the structure is retained while orientational 
ordering of the Cgg ions leads to a simple cubic structure (space group Pa3) 
with a = 13.9671(3) A at 4.5 K [63]. 

The ACeo (A = K, Rb, Cs) fullerides have attracted particular interest not 
only because of their rich structural phase diagram but also for their 
electronic, conducting, and magnetic properties which have remained the 
subject of controversy for a while [64]. A metal-insulator transition below 50 K 
is accompanied by the stabilisation of a magnetic state in the monoclinic phase 
of RbCeo and CsCeo, whereas orthorhombic KCeo remains metallic to low 




Fig. 19. Linear polymer chains formed by [2 -l- 2] cycloaddition in ACeo salts. Schematic 
drawing of chain orientations for (a) Pmnn space group (KCeo) and (b) 72/m space group 
(RbCeo, CsCeo)- The shaded bars indicate the orientation of covalently bonded polymer 
chains along the a axis by projection on the crystallographic be plane. Ci labels the bridging 
carbon atom, a is the monoclinic angle, and fi is the angle between the cycloaddition planes 
and the c axis (fi = 510, 470, and 460 for KCeo, RbCeo, and CsCeo, respectively) 




Magnetism in Fulierene Derivatives 



149 




Fig. 20. Structure of dimerised phase of RbCso stabilised by medium quenching from high 
temperature to below room temperature 



temperatures [65]. The nature of the magnetic transition has been controver- 
sial, as both quasi-one-dimensional electronic instabilities [57, 64-68] and 
three-dimensional magnetic ordering [69, 70] have been proposed. According 
to density- functional calculations [69], the rehybridisation of carbon orbitals 
at the Ceo — Ceo bonds impedes conduction along the chains and leads to 
strongly three-dimensional electronic properties with a 3D antiferromagnet- 
ically ordered ground state. On the other hand, in the original report on the 
magnetic properties of polymeric RbCeo [57], it was proposed that it behaves 
as a quasi-one dimensional metal. Measurements of the ESR susceptibility find 
a large and approximately Pauli-like spin susceptibility (xs ~ 6 x 10“^ emu/ 
mol. Fig. 21a). The temperature evolution of the static susceptibility, xs> as 
measured by EPR, revealed a sharp decrease to zero in the vicinity of 50 K, 
attributed to a spin density wave (SDW) instability of the conducting Cgg 
linear chains [57]. Antiferromagnetic correlations are also supported by the 
temperature dependence of the ESR line width (Fig. 21b), the g- factor of the 
ESR line (Fig. 21c) [57] as well as the enhanced ^^C-NMR spin-lattice 
relaxation [71]. Zero-field muon spin relaxation (ZF-pSR) studies of RbCeo and 
CsCeo samples [72-74] revealed frozen electronic moments in the low- 
temperature phase, but spontaneous Larmor frequencies which would have 
unambiguously implied long-range magnetic order (LRO), associated with 
antiferromagnetism (AF) or spin density wave (SDW) formation, were absent. 
The muons experienced a local magnetic field that peaked close to zero and 
with large spatial inhomogeneities which may be due to a number of physical 
factors, including a wide distribution of chain (Cgg),^ lengths and orientational 
disorder effects. 

In order to clarify the true nature of the long-range magnetic order in the 
ACeo polymers, antiferromagnetic resonance should be observed. Unfortu- 
nately the lack of single crystalline samples makes it necessary to rely on 
powder data. The first report of the observation of antiferromagnetic 
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Temperature (K) 

Fig. 21. Temperature dependence of the ESR (a) susceptibility, (b) line width, and 
(c) ^-factor measured in polymeric RbCeo 

resonance came from Jannosy et al. [68] who measured the ESR spectra 
of RbCeo at high frequencies (75 GHz, 150 GHz, and 225 GHz). At these 
frequencies, a broad line whose frequency is inversely proportional to the 
resonance field has been observed [64, 68]. The sublattice magnetisation was 
found to be independent of applied magnetic field up to at least 8 T, while 
magnetic fluctuations were present between 35 and 50 K. Comparison with the 
SDW system, (TMTSeF) 2 PFg, implied that the fulleride polymers are also 
quasi- ID SDW systems with 3D ordering at low temperatures. However, as 
discussed by Bennati et al. [64], unequivocal proof for an antiferromagnetic 
ground state cannot be established, as the data could be also interpreted in 
terms of a spin-glass model and final confirmation should await single-crystal 
AFMR studies. If the antiferromagnetic ground state scenario is invoked, a 
Neel temperature, = 25 K, and a spin flop field of about 0.34 T can be 
estimated. 

The low temperature state of the CsCeo polymer seems to be even more 
complicated. Recent and ^^^Cs-NMR measurements [75] have detected the 
appearance of a spin-singlet (non-magnetic) ground state below Ts=13.8 K 
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which coexists with the magnetic order that develops at The application of 
pressure first suppresses the magnetic order and a homogeneous non- 
magnetic ground state is stabilised at 5 kbar. It was proposed that the 
development of the non-magnetic phase might be correlated either with a 
structural change, tentatively ascribed to the occurrence of a spin-Peierls 
transition, or with an electronic instability. Synchrotron X-ray diffraction 
measurements have identified the appearance of a spontaneous thermal 
contraction below 14 K which provides the signature of magnetoelastic 
coupling and could trigger the occurrence of the spin-singlet ground state [60]. 
We note here that a similar situation is also encountered in the prototypical 
quasi-lD systems (TMTTF) 2 Pp 6 [76], (BCPTTF) 2 AsF 6 [77], and (TMTTF) 2 Br 
[78] where the presence of spin-Peierls fluctuations has been recently 
identified. Spin-Peierls fluctuations were detected by X-ray diffraction and 
weak diffusive lines, located midway between the main Bragg reflections, were 
observed [79]. Similar X-ray diffraction experiments on the CsCeo polymer 
should be valuable in order to clarify further the nature of its low-temperature 
spin-singlet state. 

A final comment concerns the spin susceptibilities measured by EPR for the 
ID polymer and the quenched 3D CsCeo phases. The observed Pauli-like 
susceptibilities imply that both phases are conducting [62]. However, while the 
quasi- ID polymer phase shows a metal- insulator instability, the primitive 
cubic phase reveals contrasting electronic behaviour on the local scale, as 
revealed by and ^^^Cs-NMR measurements [80]. A partial spin gap opens 
below T < 50 K, ascribed to the occurrence of localised spin singlets on a small 
fraction of the Cgg ions. ^^^Cs-NMR shows that there are two inequivalent Cs 
sites (Fig. 22), despite the lack of any observed distortions from cubic 
symmetry in structural studies [63]. For one of these sites, the NMR shift and 
(TTi)~^ follow an activated law, confirming the existence of the spin gap. 




Fig. 22. '^^Cs-NMR spectrum of the quenched primitive cubic phase of CsC6o> showing the 
presence of two inequivalent Cs lines 
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4.2 

Magnetic Properties of Na2RbC6o and Related Salts 

Most prominent among alkali metal fullerides are the superconducting salts 
with composition AsCeo- Geometrical considerations dictate that for large 
alkali ions (K"^, Rb"^, and Cs"^) occupying the small tetrahedral interstices (two 
per Geo unit with a radius of 1.12 A), the crystal structure of the A 3 C 60 salts is 
face-centred cubic (fee, space group Fm3m) [81]. Rotational motion of the 
ions is restricted and the adopted structure is understood in terms of the 
strong repulsive interactions between the alkali ions and the fulleride units. 
The latter are forced to expose towards the tetrahedral alkali ions the part of 
their quasi- spherical shape with the largest surface area, namely the hexagonal 
faces, thus leading to maximisation of the A'^-C(C6o) distances. The control of 
the structural properties of A 3 C 60 by the tetrahedral ions is manifested in the 
expanded dimensions of the fee unit cell (a = 14.240 A for K 3 C 60 , a = 14.157 A 
for Geo at room temperature) and the extreme robustness of the structure. The 
superconducting transition temperature is affected in essentially an identical 
way by both physical and chemical pressure in these systems. The N(E^) is 
modulated through the changes of the electronic overlap between the Cgg ions 
leading to the changes in the bandwidth. 

The sensitive control of the properties of intercalated fullerides by the ion 
occupying the tetrahedral interstices is immediately manifested when Na"^ are 
introduced in the tetrahedral holes of a eep array of Cgg ions. As the Na"^ ionic 
radius (0.95 A) is smaller than the size of the hole, there is evidently enough 
space for the Cgjf ions to rotate in such a way as to optimise both the attractive 
Na"^ — CgQ interactions and the Clf — Cgg contacts. The resulting structure for 
Na 2 AC 6 o (A = K, Rb, Cs) is found to be fee, comprising quasi-spherically 
disordered Cgg ions. On cooling, a phase transition occurs in the vicinity of 
room temperature to a primitive cubic structure (pe, space group Pa3) [82, 83]. 
The primitive cubic family of Na 2 (A,A')C 6 o salts displays a much steeper rate 
of change of Tc with interfullerene spacing than that exhibited hyfee fullerides, 
while at the same time, the effects of physical [84] and chemical [85] pressure 
on the superconducting properties are not identical, and chemical pressure 
suppresses Tq much faster than physical pressure does. These peculiarities in 
the electronic properties of Na 2 AC 6 o were originally ascribed to the depen- 
dence of N(Ep) on the orientational order of the Cgo molecules in the primitive 
cubic (pe) structure. However, this conjecture has proven too simplistic, as 
subsequent work [5, 86] revealed a much more complicated structural and 
electronic phase diagram for Na 2 AC 6 o (Fig- 23). Thus, it was found that, if 
special care was taken to cool the sample very slowly down to 180-200 K (with 
cooling rates as slow as 1 K/h), the ground state of Na 2 RbC 6 o below about 
250 K was that of a quasi-one-dimensional polymeric phase with monoclinic 
crystal symmetry and with a short interball centre-to-centre distance of 
'^9.38 A, which is essentially insensitive to cooling [86, 87]. 

The polymerisation of Cgjf ions in Na 2 RbCgo is characterised by a different 
structural motif than that encountered in the extensively studied RbCgg 
polymer phase, involving the formation of a single C — C bridging bond 




Magnetism in Fulierene Derivatives 



153 



pc+polymer phase-slow cooling 
pc-quenching 

0 50 100 150 200 250 300 350 

Temperature ( K ) 

Fig. 23. Phase diagram of Na2RbC6o at ambient pressure 




(Fig. 24) [87, 88]. Partial replacement of Rb by Cs to give Na 2 Rbi_xCS;cC 6 o 
(0 < X < 1) allows the formation of a family of isostructural polymer phases in 
which the amount of polymer at low temperatures strongly depends on the Cs 
content, x, as well as the cooling rate [39, 89-91]. Structural studies have 
confirmed the occurrence of polymerisation for the entire Na 2 Rbi_xCS;cC 6 o 
family at both ambient and elevated pressure. Table 1 summarises the 
dependence of the interfullerene separation, d, in the primitive cubic precursor 
phase just below the onset of polymerisation on the Cs content, x [39]. The 
Ceo-Ceo interfullerene distance is sensitively controlled by x; for instance, it 
increases by '^0.24% on going from Na 2 RbCeo to Na 2 Rbo. 3 Cso. 7 C 60 . This 
expansion is also reflected in the interchain separation of the corresponding 
polymer phases and is found to have dramatic consequences on the electronic 
and magnetic properties of the Na 2 Rbi_xCs,cC 6 o polymer phases. The 
interchain electron hopping is sensitively affected and increased quasi-one 
dimensional behaviour is encountered across the series with increasing x. 




Fig. 24. Crystal structure of the Na2RbC6o polymer. The polymer chains are aligned along the 
c axis of the cell with the fullerenes connected by single C — C bonds. Na"^ and Rb^ ions are 
shown as small and large spheres, respectively 
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Table 1. Dependence of the interfullerene separation on Cs content for the primitive cubic 
series NajEfli.^^-Cs^^Qo (0 < x < 1) at ~250 K 


Compound 


d{A) 


Na 2 CSC 6 o 


9.9754(4) 


Na 2 Rbo. 3 CSo. 7 C 60 


9.9669(2) 


Na 2 Rbo. 5 CSo. 5 C 60 


9.9601(2) 


Na 2 Rbo.sCSo. 2 C 60 


9.9519(3) 


Na 2 RbC 6 o 


9.9431(4) 



The electronic properties of Na 2 Rbi_;cCSxC 6 o samples were investigated 
using the X-band (9.6 GHz) and high-field (110 GHz) ESR techniques [92-94]. 
The ESR spectra at low temperatures when both the cubic monomer and the 
polymer phases co-exist are composed of two lines: a broad line, characteristic 
for the cubic superconducting phase and a narrow component, originating 
from the polymeric phase (Eig. 25). Deconvolution of the X-band ESR spectra 
thus enables one to determine the temperature dependence of the spin 
susceptibility of the polymeric phase. 

The temperature dependence of the intensities of the two components in 
Na 2 RbC 6 o (Fig- 26) clearly reflects the changes in the spectra: the intensity of 
the narrow component gradually increases below '^265 K at the expense of 
the high-temperature broad cubic component. However, the difference in the 
line widths of the two components is not as large as one would expect 
considering the one-dimensional versus three-dimensional structure of the 
polymeric and cubic phases, respectively. This implies that polymerisation in 
Na 2 RbC 6 o results only in a partially reduced dimensionality of the electronic 
structure. It also explains the temperature dependence of the susceptibility 
of the polymeric phase, which decreases with decreasing temperature. 
The Na 2 RbC 6 o polymer remains metallic down to 4 K without any typical 
instability, characteristic of low-dimensional metallic systems. In fact, the 
decrease in the ESR susceptibility with decreasing temperature could well 
result from lattice contraction, suggesting a 3D-like electronic structure [92]. 
It also suggests the absence of any magnetic correlations in this system and 
explains the fact that no low-temperature instability has been found in this 
compound. 

An investigation of the electronic properties of the series of Na 2 Rbi_xCS;cC 6 o 
polymers with the X-band and high-field ESR (HE ESR) techniques allowed a 
systematic determination of the low-temperature magnetic properties. A very 
sensitive way of looking at the electronic dimensionality of these systems is the 
electron spin-lattice relaxation time, Ti [95]. In the 3D metallic regime, the 
electron spin-spin relaxation time, T 2 is expected to become equal to Ti and 
the measured homogeneous line width of the ESR spectrum is given by: 
Fi ~ T 2 = ( 2 h)/{gfigAHii 2 )- It has been shown before for A 3 C 60 cubic systems 
that Ti is principally determined by the spin-orbit coupling of the conduction 
electrons to the lattice [96]. In this case, the spin-lattice relaxation time can be 
expressed by the Elliot expression [97]: 
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H(T) 



Fig. 25. Representative temperature evolution of the high-field (110 GHz) ESR spectra in 
Na 2 Rbi_;tCS;tC 6 o> showing the appearance on cooling of the narrow component, attributed to 
the polymeric phase, superimposed on the broad cubic component. The data shown here are 
for the sample with x = 0.5 



Ti = 2' 

(Ag)^ 



and 



T> 



(la) 



Ti 





T< do 



(lb) 



where Ag is the shift of the g-factor from the free electron value, ge=2.0023, T|| 
is the relaxation time of the electric conductivity, is the Debye temperature, 
and a is constant, which for isotropic metals is on the order of 1 at high 
temperatures. It should be noted that even Elliot in his original contribution 
warned about the crudity of the calculations and that numerical factors 
involved in the equations above should be treated with caution. Tanigaki et al. 
[97] connected Ti with the spin-orbit coupling constant, A, of the intercalated 
alkali ions, introducing the expression, Ag w X/E, where E is the energy 
separation of the considered valence states. This leads to 
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Fig. 26. Temperature dependence of the high-field ESR (110 GHz) intensity of the cubic 
(open squares) and polymeric (closed circles) components in Na 2 RbC 6 o. The inset shows an 
expanded version of the diagram in the temperature range 200-300 K 
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and explains why the ESR line width of the cubic line increases on moving 
from K to Rb to Cs. The reported i values of the alkali metals are: 0.3 cm“^ for 
Li, 17 cm“^ for Na, 58 cm“^ for K, 6820 cm“^ for Rb, and 9810 cm“^ for Cs. 
However, in anisotropic systems, a deviates from unity and should be 
estimated theoretically or experimentally. In a quasi-one-dimensional metallic 
system, one should also consider a cut-off of the ID correlation of electronic 
motion due to the probability of escape from the ID chain. The relaxation 
time, Tx, thus measures the mean lifetime of the ID correlation. Experimental 
work on the prototypical quasi-one-dimensional system, TTF-TCNQ, has 
shown [98] that a could be expressed as: (x w tx/t||. Increasing the degree of 
quasi-one-dimensionality by expanding the chain-chain separation in a 
conducting polymeric chain results in a reduced transverse transfer integral. 
This in turn will effectively increase a and result in line narrowing. 

In Fig. 27 we show the dependence of the high-held ESR (vl =110 GHz) 
homogeneous line width of the polymeric line on the Cs concentration in 
powder Na 2 Rbi_,cCSxC 6 o samples [93]. The homogenous line width hrst 
slightly increases on going from x = 0 to x = 0.2 and then decreases for the 
X = 0.5 and x = 0.7 samples, in striking contrast to the behaviour of the 
homogeneous line width of the cubic phase (Fig. 27 inset). The behaviour of 
the cubic homogeneous line width is in agreement with that reported for other 
A 3 C 60 systems [97] and can be understood using Eq. 2. However, in order to 
understand the behaviour of the polymeric phase with increasing Cs content x. 
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Fig. 27. Variation of the homogeneous line width of the polymeric line in Na2Rbi_j:Cs,tC6o 
samples as a function of the Cs content, x. The inset shows the dependence of the 
homogeneous ESR line width of the cubic non-polymeric phase 
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the variation of the parameter a with x should also be considered. Thus, the 
decrease in the homogeneous line width of the polymeric phase means that 
despite the large increase in the spin-orbit coupling constant, with 
increasing Cs content, a decreases even more substantially. This suggests that 
the electronic structure of Na2Rbo.3Cso.7C60 is much more quasi-one-dimen- 
sional than that of Na2RbC6o- 

We note at this point that in the analysis of the temperature dependence of 
the spin-lattice relaxation times in Na2Rbi_,cCSj.C6o. other relaxation channels, 
which may be active in these systems should also be taken into account. 
Namely, when x increases, the fraction of the polymeric phase decreases and 
consequently the coherence length within which the electron motion is 
restricted decreases with increasing x. Such a restricted type of electronic 
diffusion may additionally affect the electron spin-lattice relaxation. The 
second effect relates to the coupling of the spins in the polymer phase with 
those of the coexisting cubic phase. Determination of the individual 
contributions to the relaxation mechanism necessitates a systematic study of 
the electron spin-lattice relaxation with the pulsed ESR technique. 

It was thus of no surprise to encounter the presence of a magnetic instability 
in Na2Rbo.3Cso.7C60 [ 99 ]. On cooling, the ESR susceptibility, y(T), of 
Na2Rbo.3Cso.7C60 first increases slightly below 160 K, reaching a maximum 
between 45 and 50 K (Eig. 28 ). In a similar fashion to (NH3)K3Ceo. the increase 
of the susceptibility between 160 K and 50 K cannot be described by the Curie- 
Weiss law and its origin is ascribed to antiferromagnetic correlations in a 
narrow band metallic system. Below this temperature, y(T) suddenly drops 
and almost completely vanishes below 10 K, providing evidence of the opening 
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Fig. 28. Temperature dependence of the ESR susceptibility in the Na2Rbo.3Cso.7C60 polymer 
phase. Inset: temperature dependence of the homogeneous width of the polymer line 



of a gap at the Fermi level. The observed changes in x(T) are also accompanied 
by changes in the temperature dependence of the homogeneous line width 
(Fig. 28 inset). Hy 2 changes considerably with temperature, decreasing almost 
linearly with decreasing temperature between 200 and 45 K, where it reaches 
a broad minimum. Below 45 K, the trend is suddenly reversed and Hy 2 starts 
to increase below 25 K. The antiferromagnetic fluctuations at the Larmor 
frequency may be responsible for the increase in the homogeneous line width 
at low temperatures. 

There are at least three possible scenarios, which can lead to the vanishingly 
small ESR susceptibility observed in Na 2 Rb 0 . 3 Cs 0 . 7 C 60 at low temperature: 

(a) 3D spin density wave (SDW) ordering, as originally suggested for the 
RbCeo polymer [57], 

(b) antiferromagnetic (AFM) ordering, as seen in (NH 3 )K 3 C 6 o [49], and 

(c) transition to a non-magnetic charge density wave (CDW) state, as 
theoretically suggested for polymeric Na 2 RbCeo [100]. 

If SDW or AFM ordering occurs, the low-temperature state has a well- 
defined magnetic character and one expects to find an additional antiferro- 
magnetic resonance (AFMR) mode at high frequencies. A closer inspection of 
the low-temperature integrated HF ESR spectra of Na 2 Rbo. 3 Cso. 7 C 60 indeed 
revealed a new line emerging on the low- field side below 15 K. A low field shift 
and a broadening of the new line with decreasing temperature are character- 
istics of both the AFMR and the resonance in spin glasses [28]. Additional 
frequency-dependent experiments are needed to unambiguously assign the 
new line to AFMR. 
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In this section, we have shown that the structural, electronic, and magnetic 
properties of the Na 2 Rbi_^CS;tC 6 o polymeric phases strongly depend on the Cs 
content. The electronic structure becomes progressively quasi-one-dimen- 
sional, as x approaches 1. Consistent with this, Na 2 RbC 6 o remains conducting 
down to 4 K with no low-temperature ground-state instabilities, which may 
have led to the formation of either a SDW or a CDW state. Similar behaviour is 
displayed by the quaternary fulleride polymers with increasing interchain 
separation, despite the increased one-dimensional character, up to and 
including the phase with x = 0.5. However, in the Na 2 Rbo. 3 Cso. 7 C 60 salt, a 
pseudo-gap of magnetic origin opens at the Fermi level below 45 K. The 
presence of an antiferromagnetic resonance in high-held ESR below 
Tn 15 K provides evidence of the onset of three-dimensional magnetic 
ordering. The Na 2 Rbi_^CS;cC 6 o family of polymers thus offers a unique way of 
chemical control of the electronic properties, as the opening of the gap in this 
system of predominantly itinerant electrons is an extremely sensitive function 
of the interchain separation, which is controlled by the Cs content, x. 

5 

Summary and Future Prospects 

A common theme for all the intercalated Ceo systems described in this chapter 
has been the interplay between magnetic instabilities and the metal-insulator 
transition and suppression of superconductivity. As the Coulomb interaction, 
U between two electrons on the same Ceo molecule is by a factor of 2-3 larger 
than the width, W of the partially filled band [101] and the (UIW) ratio 
exceeds the critical value for which the systems should become Mott-Hubbard 
insulators, the fulierene intercalation compounds are highly correlated 
electron systems [5]. Then the electronic properties of the systems described 
here result from the combined effects of low-dimensional structures 
and electron correlation. A recent study of the phase diagram of ammoniated 
alkali fullerides, (NH3)K3_xRbxC6o = 1. 2, 3) [102] showed that the Neel 
temperature, T^, first increases with increasing interfullerene spacing and then 
decreases, showing a maximum at ^76 K for (NH 3 )KRb 2 C 6 o- The results were 
explained within the framework of the generalised phase diagram of the 
Mott-Hubbard transition with an antiferromagnetic ground state. 

Work on solid fullerenes and their numerous derivatives continues unabated 
to produce novel and unexpected results. Ceo-based materials display super- 
conductivity at TqS still only surpassed by the high-Tc cuprate superconductors. 
The recent findings of a variety of magnetic phenomena in systems comprising 
different fulierene structural architectures has generated a surge of interest in 
fundamental aspects of their physics and chemistry as well as in technological 
applications. Although hopes for technological applications are now shifting 
towards carbon nanotubes [103], the simplicity of the Ceo geometry and the 
variety of electronic properties found in its compounds keep this fascinating 
molecule in the mainstream of current research. The experience and the results 
obtained on Ceo-based materials represent invaluable insight, which could be 
rapidly transferred to other nanostructured materials in the future. 
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Spin ladders are low-dimensional magnetic quantum systems that were recently discovered. 
A comprehensive overview of their structural and magnetic characteristics is summarised 
here and the possibilities for supramolecular chemistry in the development of such 
fascinating materials are discussed. By using molecules as spin carrier units the interacting 
spins are located in molecular orhitals with a dominant role of n electrons. On the basis of a 
review of the few existing examples of molecular spin ladders it is clear that typical 
supramolecular and crystal engineering criteria, such as n-n overlap, S- • S and C — H- • S 
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of these magnetic quantum systems with intermediate dimensionalities between ID chains 
and 2D square lattices. 
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1 

Introduction 

Even though the field of magnetic materials has been traditionally confined to 
metals, molecule-based magnets recently appeared on the scene of magnetism. 
The rational design of molecular magnetic materials has permitted the 
modulation of their structural and magnetic dimensionalities and, on going 
from ID to 2D and finally 3D spin-coupled systems, the stabilisation of bulk 
ferromagnetism in many of these compounds has been achieved [1]. 

The main interest of these new materials arises from the fact that the 
interacting spins are located in molecular orbitals (MO) with a large or 
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dominant role of s and p atomic orbitals rather than of d and / orbitals as in 
ordinary metals and magnets. Moreover, the modulation of their crystal 
structures and thereby of their electronic and magnetic properties and 
dimensionalities can be achieved by modifying the type and/or the strength of 
intermolecular interactions present in the crystal. By changing the structure 
and components of the molecular units, modifications of this type can, in 
principle, be accomplished although the methodologies required for the 
control of the crystal packing of molecular materials as well as the factors that 
govern the magnetic interactions between purely organic molecules are still 
not completely understood. Nevertheless, the flurry of activity in this field can 
be traced to the widely held belief that even the most sophisticated properties, 
and combinations of them, can be rationally designed from the bottom up. 
This motivation was further fuelled by the increased synthetic capabilities that 
permit us to synthesise molecules with suitable structures and topologies, and 
also by the spectacular progress of supramolecular chemistry for materials 
development as witnessed in the recent years. As a result of a systematic 
synthetic effort, it has been possible to obtain different kinds of crystalline 
supramolecular organisations, using open-shell molecular building blocks, 
that show a richly diverse array of magnetic properties such as antiferromag- 
netism, ferrimagnetism, kagome antiferomagnetism, weak ferromagnetism, 2D 
and 3D ferromagnetism, and spin-ladder behaviour. The combination of 
magnetic properties with others, such as electrical or optical properties, has 
also been achieved [1]. 

This chapter is devoted to the very few molecular materials showing spin- 
ladder behaviour, describing first the structural and magnetic characteristics 
necessary to achieve this recently discovered, fascinating physical property. 

2 

Characteristics of Spin Ladders 

Spin ladders are low-dimensional magnetic quantum systems that consist of a 
finite number of strongly magnetically coupled chains of spins that are at the 
crossroad between one (ID) and two dimensions (2D) (Fig. 1). In principle, 
one might expect that a smooth crossover from chains to a plane would result 
if one assembled chains to form ladders of increasing width. But, in fact, the 
field of ladder systems emerged strongly in magnetism when Dagotto et al. [2] 
found, by numerical calculations, that the crossover from the quasi-long-range 
order in a chain (ID) of antiferromagnetically coupled S = 1/2 spins to the 
true long-range order that occurs in a plain (2D) is not at all smooth. 

The ladders with S = 1/2 spins that show antiferromagnetic isotropic 
coupling between the nearest-neighbours, and that are generally named as 
S - 1/2 Heisenberg antiferromagnetic spin ladders are especially interesting. 
The Hamiltonian in Eq. (1), defines such magnetic systems, 

L L 

H = — 7|| Si^a ■ S;+I,a — J± S,j • Sj_2 

a=l,2 i=l i=l 



(1) 
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Fig. 1. Schematic illustration of different kinds of low dimensional magnetic Heisenberg 
systems with antiferromagnetic couplings: (a) ID chain; (b) two-leg and three-leg spin 
ladders where J± is the coupling along the rungs and /y is the coupling along the chains; 
(c) 2D square lattice with different exchange coupling constants, / and /'. Black dots denote 
S = 1/2 spin-containing units and continuous and dashed lines represent antiferromagnetic 
interactions of different strength 



where is the spin operator at the site i (i - 1, 2,...,L) on the leg a 

(fl = 1, 2 . . .) of a ladder with L rungs. J±_ and /|| denote the intra- and inter- 
rung exchange couplings, respectively. With this definition, /x and /|| should 
be negative for antiferromagnetic interactions. 

The magnetic properties of such systems are particularly interesting as a 
result of the increased importance of the quantum mechanical effects 
operating within them. Theoretical calculations have predicted that spin 
ladders with an even number of legs have a spin-liquid ground state, so called 
because of their purely short-range spin correlation along the legs. These even- 
leg ladders consist of spin singlet pairs with a spin-spin correlation distance 
along the legs that show an exponential decay produced by the presence of a 
finite spin gap. By contrast, a ladder with an odd number of legs behaves quite 
differently and displays properties similar to those of a ID antiferromagnetic 
Heisenberg chain at a low thermal energy, namely, gapless spin excitations and 
a power-law falloff of the spin-spin correlations which are magnetically 
ordered [3-10]. 

The magnetic properties of two-leg S = 1/2 antiferromagnetic spin ladders 
can easily be described in the simple limit where intra-rung exchange coupling 
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(/x) is much larger than the exchange coupling along the legs (/||); /x ^ J\\. 
The energy of the ground state, Eg^, in this limit is approximately Eg^ = -3/4 
J^N, where N is the number of rungs and 3/4/x is the energy of each rung 
singlet state. In this particular case, the rungs of the ladder interact only 
weakly with each other and the dominant spin configuration in the ground 
state is that with the spins on each rung forming a spin singlet. Therefore, the 
ground state has a total spin S = 0. To produce the lowest excitation in the 
infinite ladder to a configuration with a total spin S = 1, one of the rung 
singlets of the ladder must be promoted to an S = 1 triplet (Fig. 2) with a finite 
energy that must overcome the finite spin gap of the system. Moreover, the 
coupling along the legs creates an energy band of S = 1 magnons with a 
dispersion of about 2/|| leading to an exponential decrease in the magnetic spin 
susceptibility as the temperature is lowered until kT is below the spin gap 
energy (Fig. 3). 




s = o 



S=1 



Fig. 2. Schematic illustration of the lowest excitation to a configuration with a total spin 
S = 1 in a two-leg S = 1/2 Heisenberg antiferromagnetic ladder with /x > J\\ 




k 

4 



2J// 



Fig. 3. Evolution of the energy levels in a spin ladder when the interaction along the legs is 
turned on 
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Another interesting prediction for two-leg ladders is that light hole doping 
can even lead to high Tc superconductivity on basis that effective attraction 
between extra holes may arise from the magnetic interactions between the 
spins of the ladders [2]. Odd-leg ladders, on the other hand, must have 
different behaviour because singlet pair formation on their rungs is not 
possible [5]. In fact, they display properties similar to those of single chains at 
low energies. 

“Ideal” spin ladders are considered to be those in which the exchange 
coupling along the rungs is very similar to the exchange coupling along the 
chains, that is }±/J\\ ~ 1- Therefore, regarding the two-leg spin ladders, such 
systems are intermediate between two limiting cases: 

a) isolated dimers, when J±/}\\ oo, and 

b) isolated ID chains, when }±/J\\ 0. 

Although these limiting cases have been analysed from the theoretical point of 
view, spin ladders of experimental interest are those having both exchange 
coupling constants of the same order of magnitude, that is 0.1 < }±/J\\ < 10. 
Another factor that should be taken into account on considering a system as a 
spin ladder is that ladders should be quite well isolated one from another, 
since appreciable interladder coupling (/' in Fig. 4) can promote a quantum 
phase transition between the spin liquid ground state and a magnetically 
ordered state. Troyer et al., established a critical ratio of the interladder to 
intraladder coupling of /7/x =0.1 that separates the spin liquid from an 
antiferromagnetic ordered state [11]. 

The appealing properties theoretically predicted for spin ladders have only 
recently been found in very few materials. The first examples were connected 
chains of transition metals in inorganic oxides, such as the series of even- and 
odd-leg ladder structures provided by the cuprates Sr„_iCu„+i02„ {n = 3, 5) 
[12] or the two-leg ladders (VO) 2 P 207 [13] and LaCu02.5 [14]. Another system 
that has been described with various models including a spin ladder one is 
Cu 2 (C 5 H 2 N 2 ) 2 Cl 4 with interacting Cu atoms being in the strong coupling limit 
(/x//|l = 5.5) [9]. However, it is important to note that the identification of 




Fig. 4. Scheme of two-leg ladders showing intra-ladder (/x and /||) and inter-ladder (/') 
exchange couplings. Black dots denote S = 1/2 spin-containing units 
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these compounds as spin ladders has been made, primarily, on the basis of 
their macroscopic characterisation (magnetic susceptibility data) and it has 
been shown later on that microscopic probes, such as neutron scattering, 
nuclear magnetic resonance (NMR), and muon spin resonance (pSR), are 
necessary to establish the exact nature of the ground state. In fact, although the 
family of cuprates containing Sr has been confirmed as real spin ladders 
by different experiments [15], NMR and pSR studies on the La-containing 
cuprate [16, 17] and inelastic neutron scattering measurements on the 
(VO) 2 P 207 [18] showed that both compounds are not real spin ladders. 

The search for new compounds showing spin ladder properties is a very 
important task not only for understanding the fascinating quantum ladder 
physics but also for developing new materials, especially superconductors. 
Recently, some cuprates that behave as spin ladders have shown supercon- 
ducting properties upon doping [19]. 

3 

Molecular Spin Ladders 

In the search for new spin ladder compounds, those belonging to the family of 
molecular solids can be of great interest since they have a new molecular 
architecture for interacting spins. In fact, in the field of molecular materials 
science and from the viewpoint of chemists, one of the most exciting topics is 
solid state supramolecular synthesis; that is to combine appropriate molecular 
building blocks in a predesigned way so as to endow the resulting 
supramolecular assembly with desirable physical properties. Because of their 
structural flexibility, molecular solids offer a field of choice to finely tune 
interesting physical properties such as those of the spin ladder materials. It 
should be emphasised that, in materials constructed from open-shell 
molecules, the modulation of their structural, electronic, and magnetic 
dimensionalities has permitted the realisation of molecular compounds 
showing a large variety of physical phenomena ranging from low-dimensional 
metals and superconductors to organic ferromagnets [1, 20]. As the interacting 
spins in these materials are located in molecular orbitals (MO), the molecule- 
based magnetic systems should be commonly understood by using the MO 
concepts and to analyse the magnetic interactions it is not only the 
connectivity between atoms that we should be considering, but also the 
connectivity between orbitals. Thus, the two most important issues to be 
considered for the creation of molecular spin ladders are: 

a) the election of the spin-containing building blocks, or open-shell repeating 

units (molecules), and 

b) the assembling of such molecular building blocks or repeating units leading 

to correct connectivity of orbitals. 

The construction of the simplest molecular spin ladder - a two-leg ladder - 
can be conceptually achieved following two different strategies (Fig. 5). The 
first one consists of the connection of two molecular ID S = 1/2 spin chains, 
one next to the other, and the second one to assemble an infinite number of 
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(a) (b) 



Fig. 5. Illustration of the two strategies that can be used in the construction of molecular 
two-leg ladders: (a) assembling of ID S = 1/2 chains, and (b) connecting of S = 1/2 dimers 



units with two interacting S - 1/2 molecules (dimers). Consequently, the 
election of the appropriate open-shell building blocks, able to form ID chains 
or dimers and that are also able to interact with the neighbouring units giving 
rise to antiferromagnetic interactions, is of primary importance. Although the 
few molecular spin-ladders known so far have been obtained by serendipity, 
the study of these compounds shows that both conceptual strategies have been 
involved in their construction [21-24]. In the next sections the very few 
already known molecular spin ladders constructed following the two men- 
tioned strategies will be described. 



3.1 

Molecular Spin Ladders by the Assembling of Chains 

To follow the first strategy - the assembly of chains - one of the most 
promising building blocks are tetrathiafulvalene (TTF) derivatives (Scheme 1) 
since it is well known that these planar 7i-electron donors form stable radical 
cations, by donating an electron from the HOMO, and have a large tendency to 
stack, thus forming ID chains [20]. 

This packing promotes the overlap of directional 7i-orbitals along the 
stacking direction that results in anisotropic band structures. Furthermore, 
the transfer integrals (t) between neighbouring molecules, which describe the 
intermolecular electronic interactions, are comparable or even smaller than 




Scheme 1. 
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the on-site Coulombic repulsion (U) and thus the electron wave function is 
very often on the borderline between the localised and delocalised regimes. 
When localisation of the electrons is promoted, by lowering the t/U ratio, a ID 
spin chain can be generated that could provide the legs of the ladders if 
interchain magnetic interactions are adequate. In fact, the sulfur atoms of TTF 
derivatives, that have a large contribution to the HOMO of the molecule, 
readily establish intercolumn structural and electronic interactions, thus 
promoting enhanced electronic dimensionality in a great number of metallic 
charge-transfer complexes and radical- cation salts derived from these mainly 
planar molecules [20, 25]. Therefore, S- ■ S interactions could provide the 
rungs of the ladder by connecting the ID chains. Isotropic exchange coupling 
constants, /„ among open-shell molecules in these TTF-based molecular 
compounds are related with the transfer integral, t„ of the corresponding 
intermolecular contact and the on-site Coulombic repulsion, U, of the 
molecule by Eq. [2] having predominantly an antiferromagnetic character: 



Due to the above-mentioned characteristics, structural and electronic ID 
and 2D layered, TTF-based organic conductors are quite common in the 
literature [20], although the achievement of solids with an intermediate 
structural dimensionality - that is a finite number of assembled stacks forming 
ladder-like structures - has only recently been accomplished [21, 22]. Only 
when the proper array of diamagnetic counterions cut a 2D layered TTF 
structure with localised electrons, will a spin ladder molecular material result 
in the donor sublattice. The two existing examples that will be described in the 
following are two-leg spin ladders obtained in such a way since the 2D donor 
layers are cut every two chains by diamagnetic counterions. 

As a first example, the purely organic molecular two-leg spin ladder 
compound prepared in our group by using as a building block the TTF 
derivative dithiophenetetrathiafulvalene (DT-TTF) and the Au complex of the 
maleonitrile dithiolate [Au(mnt) 2 ]~ as a magnetically innocent counterion will 
be described and analysed (Scheme 2) [21]. This planar diamagnetic coun- 
terion has the same size as the donor and also exhibits a great tendency to 
stack, so forming ID chains, and therefore being complementary to DT-TTF. 

DT-TTF and [Au(mnt) 2 ]~ molecules crystallise together to form the mixed 
valence salt (DT-TTF) 2 Au(mnt )2 in which segregated DT-TTF (donor, D) and 
[Au(mnt) 2 ] (acceptor. A) stacks with a herringbone pattern are present (Fig. 6). 




( 2 ) 




DT-TTF 



Au(mnt)2 



Scheme 2. 
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Fig. 6 . (a) Crystal packing of (DT-TTF)2Au(mnt)2 showing the ladder structure formed by 
the organic DT-TTF donor stacks related by a two-fold screw axis; (b) projection of the 
crystal structure along the b axis, where short S- • -S contacts are indicated by dotted lines 
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At room temperature, the DT-TTF stacks are arranged in pairs related by a 
twofold screw axis, that alternate with single stacks of [Au(mnt)2] units along 
the a-c direction. The pairs of organic donor stacks are the legs of a two-leg 
ladder structure since they are strongly linked by three inter-stack S- • S close 
contacts (interaction II in Fig. 6b), providing the rungs for the ladder. In fact, 
the calculated transfer integrals show that the two legs of the ladder are 
interacting strongly since the value of the transfer integral coupling the two 
paired chains of donors (tx = 21 meV) is very similar to that of the transfer 
integral up the chain (f|| = 36 meV), the on-site Coulombic repulsions of 
DT-TTF being of the order of 1 eV. 

The structural arrangement of donors in [(DT-TTF)2Au(mnt)2] seems to be 
governed by strong intermolecular interactions since it is the same as that 
found in the 2D organisation of the neutral donor in the solid state (Fig. 7) [25, 
26]. Therefore, the structure of the [(DT-TTF)2Au(mnt)2] salt can be described 
by separating, in pairs of stacks, the 2D arrangement of the neutral DT-TTF 
with the stacks of anions. Thus, by using a diamagnetic anion of similar size 
and shape to that of the donor molecules, we have succeeded in the realisation 
of an intermediate dimensionality between ID and 2D for a molecular 
compound. 

This structural ladder should have localised interacting electrons in order 
to behave as a two-leg spin ladder; otherwise it could behave as a ID 
metal. Indeed, electrical transport measurements (Fig. 8) show that 
below 220 K the electrons are localised and. X-ray diffuse scattering 
experiments (Fig. 9) confirm that a dimerisation is achieved below 220 K in 
the chain direction b. 

In addition, the finite width of the observed X-ray diffuse lines shows that 
the dimerisation locally breaks the twofold screw axis symmetry that relates 
the two DT-TTF chains that form the ladder. In this picture, each dimer of 
DT-TTF molecules has one localised electron (Fig. 10). The dimers 
[(DT-TTF)^*] constitute the spin carrying units of the spin ladder. 

Susceptibility data on powdered and single crystal samples of the salt [(DT- 
TTF)2Au(mnt)2] show a characteristic behaviour of localised spins with 
strong AF interactions, which displays a spin gap (Fig. 11). The exponential 
decay of susceptibility data at low temperature (from 8 to 45 K) was 
successfully fitted with the Troyer expression [Eq. (3)] [8] for a two-leg spin 
ladder system with the energy gap, A/k, in the spin excitation spectrum found 
to be equal to 78 K. 



Ziadder = «7’ exp(-A//cT) (3) 

The exchange interactions along the legs {}\\/k = -83 K) and the rungs 
(7x/k=— 142K) of the ladder were extracted by fitting of the whole 
susceptibility data (Fig. 11) with Eq. (4), which takes into account the Curie 
contribution due to paramagnetic impurities (defects) in the crystals and with 
Eq. (5). The two-leg ladder model [Eq. (5)] used for such a fitting was that 
developed by Barnes and Riera, where Ci, C 2 , C3, C4, C5, and Ce, depend on /x and 
/|l values [13b[. 
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Fig. 7. Crystal packing of DT-TTF: (a) view along the long axis of the molecules showing the 
organic DT-TTF stacks related by a two-fold screw axis; (b) projection of the crystal 
structure along the b axis. Short S- • -S contacts are indicated by dotted lines 
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A remarkable feature of this two-leg spin ladder is the fact that the 
ratio is 1.7, which is close to that of an “ideal” spin ladder. The spin gap has 
also been calculated from the resulting values of /|| and /x with the theoretical 
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Fig. 8. Electrical resistivity p of [(DT-TTF)2Au(mnt)2] as a function of the reciprocal 
temperature (left) and its first derivative (right) measured along the needle axis b by the 
standard four-probe method 




Fig. 9. Thermal dependence of the HWHM along the b axis of the ll 2 b* diffuse lines of 
[(DT-TTF)2Au(mnt)2]. The inset shows the profile along b* of such a line (y is the number of 
counts in arbitrary units) 

Eq. (6) [8], giving A/k = 83 K, which is in good agreement with the previous 
value. 

a = i;j-|;iiI+/^/2/x 



( 6 ) 
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Fig. 10. Schematic illustration of the two possible two-leg ladders formed by the dimeri- 
sation of DT-TTF stacks in [(DT-TTF)2Au(mnt)2] 
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Fig. 11. (a) Temperature dependence of the paramagnetic susceptibility (•) of the 
(DT-TTF)2Au(mnt)2 salt. The solid line is the fit by Eq. ( 2 ) to a contribution of a spin ladder 

model (-•-•-•) and a Curie tail ( ). (b) Temperature dependence of the paramagnetic 

susceptibility (A) of the (DT-TTF)2Au(mnt)2 salt measured by ESR on a single crystal 



The EPR signal of a single crystal of [(DT-TTF)2] [[Au(mnt)2] has the typical 
parameters (g value and AHpp) of TTF derivatives. The intensity of this FPR 
signal, which is proportional to the paramagnetic spin susceptibility, also 
exhibits the same thermal dependence as the static spin susceptibility 
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(Fig. 11b), confirming that the spin-ladder behaviour occurs on the organic 
[(DT-TTF)^*] stacks. 

The temperature dependence of the magnetic susceptibility clearly inferred 
that the [(DT-TTF) 2 Au(mnt) 2 ] salt is a purely organic molecular system with a 
two-leg spin ladder configuration. Nevertheless, since the calculated inter- 
ladder transfer integral t' - 6 meV, although small, is not negligible and could 
imply a quite sizeable ladder-ladder interaction (/ in Fig. 4) that affects the 
nature of the low temperature magnetic state, independent magnetic 
measurements were performed. Thus, zero-field and longitudinal field 
^■^SR measurements as a function of temperature have confirmed that 
[(DT-TTF) 2 Au(mnt) 2 ] is a molecular material with a real two-leg spin ladder 
configuration since the results obtained corroborate the absence of a magnetic 
ordering, as expected for even-leg spin ladders with the quantum spin-liquid 
state realised [27]. 

There is a second example of a TTF-based organic spin ladder that has been 
reported by Komatsu et al. [22]. The compound, (BEDT-TTF)Zn(SCN) 3 , 
belongs to the fruitful family of bis(ethylenedithio)tetrathiofulvalene 
(BEDT-TTF) radical ion salts (Scheme 3). 

The structure of this compound consists of two parallel uniform stacks of 
completely ionised BEDT-TTF donors along the c-axis that are isolated from 
the neighbouring pairs of stacks by the diamagnetic Zn(SCN) 3 ^ anions, as 
shown in Fig. 12. The (BEDT-TTF)+* stacks form a structural two-leg ladder 
since they are connected by side-by-side S- • S interactions (interaction 3 in 
Fig. 12) between BEDT-TTF molecules as well as by another weaker face-to- 
face interaction (interaction 2 in Fig. 12). 

The low conductivity and the semiconducting behaviour of this salt denote 
that electrons are localised in the BEDT-TTF molecules and, therefore, the 
(BEDT-TTF)+ stacks can be regarded as uniform S = 1/2 chains. Based on the 
calculated overlap integrals (Table 1) and applying Eq. (2), with an assumed 
value of 1.0 eV for U, the authors argue that there are four possible exchange 
coupling constants; one along the chain (/i/k = -966 K), two interchains 
ihlk = —5.2 K and J^/k = —86.4 K), and another, very small {JJk - —0.09 K), 
with an interladder character. Nevertheless, only two of these magnetic 
couplings can be considered as effective: those along the legs (/i = /||) and 
those along the rungs {J3 = J±_) of the ladder. J2, that is 1/17 of J3, can be 
regarded as a very weak perturbation to the ladder interaction (Fig. 13). On the 
other hand, the ladders are well isolated from each other since the inter-ladder 
interactions J4 are smaller that 0.1 K. 




BEDT-TTF 



Scheme 3. 
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Fig. 12. X-ray structure of (BEDT-TTF)Zn(SCN )3 following Komatsu et al. [22a]. (a) View 
along the molecular long axis of BEDT-TTF; (b) view along the c axis. Arrows indicate the 
interactions within (interactions 1-3) and between (interaction 4) the ladder 



The fit of the experimental susceptibility data in the low temperature range 
(Fig. 14) to the Troyer expression [Eq. (3)] gives an energy gap of A/fc = 
340 K. Nevertheless, in this case the energy gap calculated using Eq. (6) from 
the 7|| and /x values, which were estimated from the transfer integrals, gives A/ 
k - 4520 K. This value shows a large discrepancy with the experimental one 
(more than one order of magnitude). Since Eq. (6) is derived from a 
perturbation of the dimer limit /x ^ }\\ [8] and the present case is in the 
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Table 1. The overlap integrals (S) and the exchange energy (/) 



Interaction* 


1 

o 


;rVK 


1 


-20.4 


-966 


2 


1.5 


-5.2 


3 


-6.1 


-86.4 


4 


0.2 


0.09 



Interactions 1-4 correspond to those in Fig. 12. 




Fig. 13. Schematic representation of the two-leg spin ladder formed by the two chains of 
(BEDT-TTF)+* in (BEDT-TTF)Zn(SCN )3 (reprinted with permission [22a]) 




7/K 



Fig. 14. Plot of xT vs. T for the (BEDT-TTF)Zn(SCN )3 salt (open circles). The solid line is the 
best fitting curve to the spin ladder model [Eq. (3)]. Temperature dependence of ESR line- 
width (closed circles) of an oriented single crystal is also shown (reprinted with permission 
[22a]) 

opposite limit of weakly interacting antiferromagnetic chains, another 
approximation to the spin gap, that takes into account this characteristic, 
was used. Such an estimation is given by a Monte Carlo calculation [28] as 
A = 0.41 /x, which yields an energy gap A/k = 35 K; a value that is one order 
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of magnitude smaller than the experimental one. The authors stated that this 
result could probably be accounted for by the frustrating effect of the 
antiferromagnetic interaction J 2 (Fig. 13). It must be pointed out that the }±/}\\ 
ratio of 0.1, estimated from the transfer integrals, place this two-leg spin 
ladder far from an “ideal” spin ladder, more precisely at the limit in which it 
could also be described as pairs of weakly coupled ID chains. 



3.2 

Molecular Spin Ladders by the Assembling of Dimers 

The second strategy to the two-leg ladder systems - the assembly of an infinite 
number of dimers - has been also successful in the achievement of molecular 
spin ladders. By contrast with the previously described examples, in all known 
cases the spin carrying units are not purely organic but metal dithiolene 
complexes. Nevertheless, the magnetic interactions present in such complexes 
are mediated by overlaps between the p or n SOMO of the neighbouring spin 
carrying units. 

Fourmigue et al. [23] have found spin ladder properties in some members 
of the family of charge-transfer complexes formed by the acceptor tetra- 
fluorotetracyanoquinodimethane (TCNQF4) and the flexible organomefallic 
cyclopentadienyl/dithiolene (dithiolene = dmit^“, dmid^“, dsit^“) Mo and W 
complexes as the donors (Scheme 4). These organometallic complexes are 
good donors forming stable radical cations that have the SOMO orbital 
delocalised in both the CP 2 M fragments and the dithiolene ligands. In every 
complex, the MY 2 C 2 (M = Mo, W; Y = S, Se) metallacycle is folded along the 
Y- ■ Y axis (see Fig. 15) and the folding angle varies strongly with the nature of 
the metal and the dithiolene ligand. Depending on the folding angle, the 
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S' 





dmit^' 



dmid^' 



dsit^ 




F F 
TCNQF4 




Cp2M(dithiolene) 



Scheme 4. 
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Fig. 15. Side view and top view (relative to the dmid^ ligand) of the [{Cp2Mo(dmid)+*}2] 
dimer in the X-ray structure of [Cp2Mo(dmid)] [TCNQF4] after Fourmigue et al. [ 23 ] 

contribution of each fragment orbital to the SOMO varies and as a 
consequence so do the electronic characteristics of the building blocks. 

The analysis of the X-ray structures of different compounds formed by 
oxidation of these donors reveal that they have a remarkable tendency to form 
paramagnetic [{Cp2M(dithiolene)+*}2] radical cation dimer motifs which 
interact to give different structural organisations [ 23 , 29 ]. Therefore, these 
complexes are good building blocks to form spin ladders following the second 
strategy proposed if the proper interaction between them is achieved. 

The strong electron acceptor TCNQF4 oxidises the organometallic 
[Cp2M(dithiolene)j donors that are organised in the resulting solid as head- 
to-tail centrosymmetrical dimers with parallel dithiolene planes. The reported 
[Cp2M(dmid)][TCNQF4] (M = Mo, W) complexes are isostructural being the 
[Cp2M(dmid)]+* radical cations dimerised by dmid/dmid n-n overlap 
interactions (Fig. 15 ). Such dimers constitute the rungs of a ladder since they 
are strongly interacting due to short S- • S contacts, giving rise to the two-leg 
ladder type structures that are isolated from each other by the chains of 
TCNQF4 (Fig. 16 ). Since TCNQF4^* radical anions are strongly dimerised 
into diamagnetic [TCNQF4]2^ species, only the [Cp2M(dithiolene)]+* cation 
radicals are responsible for the magnetic properties shown by these solids. 

In accord with the described structural characteristics, the exponential 
decay of susceptibility data below the maximum can be fitted to the Troyer 
equation [Eq. ( 1 )] giving different spin gaps for both compounds (Fig. 17 ); 
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Fig. 16. (a) X-ray structure of [Cp 2 Mo(dmid)] [TCNQF4] after Fourmigue et al. [23]. 
(b) View of the ladder structure formed by the S = 1/2 cations in the 
[Cp 2 Mo(dmid)][TCNQF 4 ] salt (reprinted with permission [23]) 



A/fc = 74 K and A/k = 13 K for the Mo and W salts, respectively. Thus, the 
replacement of W for Mo, while not altering the structural characteristics of 
the supramolecular structures, modifies the electronic characteristics of the 
[{Cp 2 M(dmid)+*} 2 l building block and as a consequence the electronic 
structure of the molecular spin ladder. 

The authors extracted the values of /|| and combining the Eqs. (6) and (7). 

The obtained exchange coupling constants for the Mo salt were /x = 107 K 
and 7|| = 41 K and for the W salt J± - 23 K and /|| = 16 K. The relative values 
of these constants are in accordance with the smaller dithiolene contribution 
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Fig. 17. Plot of magnetic susceptibility vs. T of the [Cp 2 M(dmid)] [TCNQF4] (M = W, Mo). 
Insets show the fits of the low temperature data to [Eq. (3)] [23] 

to the SOMO for the W complex and are also in very good agreement with the 
calculated interaction energies within the dimers. Remarkable also are the 
7x//|| ratios of 2.6 and 1.5 for Mo and W salts, respectively, since they are close 
to that of an “ideal” two-leg spin ladder. 

r(ymax)/|/±l = 0.625 - 0.018357|j//x -f 0.2532(7||/7x)" (7) 

The structure of the complex [Cp 2 W(dmit)] [TCNQF4] also presents 
interacting dimers although with another type of dimer association, due to 
S- ■ S interactions through their p orbitals instead of n-n overlaps (Fig. 18). 
The interaction of the dimers gives rise to a structural ladder motif (Fig. 19). 

The susceptibility data (Fig. 20) were analysed as a two-leg spin ladder 
system, as in the other complexes, giving a spin gap value of A/k = 40 K [23]. 

The last known example of molecular spin ladder is provided by another 
dithiolene complex, the anionic 71 -radical [Ni(dmit) 2 ]^* (Scheme 5) [24]. This 
planar S = 1/2 ion also forms dimeric units that interact in a ladder-like 
structure in the salt formed with a paramagnetic counterion, the p-N- 
ethylpyridinium (X-nitronyl nitroxide (p-EPYNN]). 

In the crystal structure of [p-EPYNN] [Ni(dmit) 2 ] salt, the open-shell 
p-EPYNN cations form chains along c axis that isolate the chains of dimers 
formed also along c axis by the dithiolate radical-anions (Pig. 21). The dimers 
of [Ni(dmit) 2 ]^* are formed by the ordinary plane-to-plane n-n overlap and 
the connection between the dimers are provided by very short (3.31 A) S- • S 
contacts that take place along the c axis. 

In contrast to the previous [Cp 2 M(dithiolene)] [TCNQP 4 ] compounds, the 
crystals of [p-EPYNN] [Ni(dmit) 2 ] contain two different magnetic subsystems: 
the p-EPYNN cation radical chains and the one dimensional ladder-chain 
of [Ni(dmit) 2 ]^* anion radicals. The chains of p-EPYNN cation radicals 
are arranged in very similar way to that observed in the salt 
[p-EPYNN] [Au(dmit) 2 ] in which only the organic radicals are responsible 
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Fig. 18. Side view and top view (relative to the dmit^ ligand) of the [{Cp2W(dmit)+*}2] 
dimer in the X-ray structure of [Cp2W(dmit)][TCNQF4] after Fourmigue et al. [ 23 ] 




Fig. 19. View of the ladder-like structure formed by the S = 1/2 cations in the 

]Cp 2 W(dmit)] [TCNQF 4 ] salt after Fourmigue et al. [23]. Intra-dimer interactions: , 

inter-dimer interactions: 



for the magnetic properties [24b]. For this compound, the paramagnetic 
susceptibility can be fitted to the expression for a ID ferromagnetic chain [Id]. 
The authors analysed the susceptibility data of the [p-EPYNN][Ni(dmit) 2 ] 
compound as originating from two independent magnetic subsystems 
(Fig. 22). Thus, they considered the observed two distinct temperature 
dependencies, one above 150 K and the other below 40 K, to originate from 
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Fig. 20. Plot of magnetic susceptibility vs. T of the [Cp 2 W(dmit)] [TCNQF4]. Insert show the 
fit of the low temperature data to [Eq. (3)] [23] 
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these two subsystems. The low temperature region, where an increase of XpT 
with decreasing temperature is observed, has been ascribed to the p-EPYNN 
cation radical chains showing a ID ferromagnetic chain behaviour, since the 
behaviour is very similar to that found in the [p-EPYNN] [Au(dmit) 2 ] 
compound (Pig. 23) which shows the same structural arrangement of radicals. 
The high temperature region, where there is an decrease of ypT with 
decreasing the temperature is ascribed to a spin ladder behaviour for the chain 
of dimers formed by the [Ni(dmit) 2 ]^* anion radicals. The fitting of the 
susceptibility data shown in Pig. 22 to the sum of expressions for a spin ladder 
[Eq. (3)] and for a ID chain with ferrromagnetic coupling [Id] gave a spin gap 
of A/k = 940 K for the Ni(dmit )2 ladder and a ferromagnetic exchange 
interaction of Jlk - 0.16 K for the [p-EPYNN] chains. 

Extended Huckel calculations of the overlap integrals between adjacent 
molecules along the legs and in the rungs gave a small value (0.13 x 10“^) for 
the overlap along the [Ni(dmit) 2 ] chains and a significantly larger value 
(-17.0 X 10“^) for the plane-to-plane overlapped molecules in the rungs. The 
evaluation, from these values, of the exchange energies Jx and /|| of this system 
gives a }±/}\\ ratio of 10^ that is far away from the limit of “ideal” spin ladders 
and is more in accordance with a non-interacting dimer description. 
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Fig. 21. Crystal structure of [p-EPYNN][Ni(dmit) 2 ] salt after Imai et al. [24b]. (a) Both 
chains of cations and anions are shown. Short S- • -S interdimer contacts are denoted by 
doted lines, (b) Molecular arrangement in the p-EPYNN chain 
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Fig. 22. Plot of susceptibility of [p-EPYNN] [Ni(dmit) 2 ] vs. T [24] . The solid curve is the fit to 
the sum of ladder and ferromagnetic chain equations and the black and grey regions 
represents the contribution from [p-EPYNN] and [Ni(dmit) 2 ], respectively. The broken 
curve is the best fit using a simple dimer model 
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Fig. 23. Plot of susceptibility of [p-EPYNN][Au(dmit) 2 ] vs. T [24]. The solid curve is the fit 
to the ferromagnetic chain equation 



Nevertheless, the best fitting of the susceptibility data to a simple dimer model 
(broken line in Fig. 22) with a A/k = 870 K, manifests significant deviation 
from the experimental plots and the observed relatively large conductivity for 
a half-filled insulator, however, suggests that the interaction along the legs may 
possibly be larger than that estimated. 

4 

Perspectives 

It is important to emphasise that although spin ladder physics has emerged 
only very recently as a fertile subfield of condensed matter studies, molecular 
spin ladders have appeared at the initation of this area of research. This fact 
reveals both the high level reached by the field of molecule-based magnets and 
the spectacular progress of supramolecular chemistry for materials develop- 
ment witnessed in the recent years. In fact the examples presented here, show 
that supramolecular chemistry is a very powerful tool to construct solid 
compounds with tailored physical properties such as spin ladders. Typical 
supramolecular and crystal engineering criteria like n-n overlap, S- ■ S 
interactions and complementary nature of size and shape, used with success 
in the tuning of dimensionality and properties of molecular conductors and 
magnets have also proven to be useful in the construction of intermediate 
magnetic dimensionalities which are the spin ladders. Although it is not an 
easy task, application of these criteria will help in increasing the number of 
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examples of molecular spin ladders that will aid in the understanding of low- 
dimensional spin systems exhibiting spin gaps. Another interesting aspect of 
the molecular spin ladders is that the interacting electrons are mainly n- 
electrons since they are located in molecular orbitals with a large or dominant 
role of s and p atomic orbitals. Important challenges in this field that deserve 
to be reached with molecular compounds are the construction of spin-ladders 
with different spin gap values able to be measured experimentally. Also 
important is the preparation of spin ladders with different degrees of ladder- 
ladder interactions and molecular spin ladders with more than two legs. 
Molecular compounds with such characteristics will help in the corroboration 
of theoretical findings about the limits of ladder-ladder exchange coupling 
interaction values that can destroy the ladder characteristics, as well as to 
verify the dependence of the bulk magnetic properties with the odd/even 
number of legs in the ladders. 
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